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1.  FORWARD 

Advanced  Countermeasure  Systems  (ACM  Systems)  has  completed  the  Phase  1  SB1R  awarded 
on  August  6,  1990  under  the  U.S.  Army  Research  Office  contract  number  DA AL03-90-C-0025. 
This  document  is  the  final  report  of  this  Phase  I  study  effort. 

The  original  RFP  solicitation,  A90-099  Adaptive  Antennas  and  Processing,  stated  the  objective 
to  be  "Perform  research  to  provide  technology  for  adaptive  antenna  arrays  and  processing  for 
Army  tactical  ground  radio  communications."  The  solicitation  went  on  to  describe  the  general 
requirements  to  be  "Army  tactical  ground  radio  systems  operate  in  an  environment  of  strong 
and  complex  interference.  Adaptive  antennas  have  the  potential  to  provide  20-30  db  gain. 
However  the  systems  must  have  the  ability  to  converge  rapidly  in  the  presence  of  worst  case 
countermeasures  and  must  be  small  since  Army  communication  terminals  need  to  be  mobile. 
Research  is  also  needed  on  the  operation  of  such  antennas  in  wideband  systems  such  as  spread 
spectrum,  frequency  hipping  and  direct  sequence.  Techniques  must  also  be  found  for  the 
utilization  and  integration  of  adaptive  antennas  in  Army  networks."  The  specific  Phase  I  goal 
was  stated  to  be  "  ...  demonstrate  the  feasibility  of  small  adaptive  antennas  under  conditions 
of  high  interference  (jamming)." 

These  objectives  and  goals  of  Phase  I  have  been  met  and  exceeded  by  ACM  Systems,  as  is 
evident  in  the  final  report  which  follows.  A  receive  adaptation  algorithm  which  is  tailored  for 
spread  spectrum  communications,  frequency  hopping  and  direct  sequence,  has  been  developed. 
A  receive  architecture  which  is  compatible  with  existing  Army  radios,  yet  flexible  enough  to 
accomodate  a  variety  of  algorithms,  using  commercially  available  and  proven  components,  and 
with  minimum  complexity  to  allow  size  and  power  reductions  in  future  designs  to  result  in 
compatibility  with  the  Army’s  mobile  applications,  has  been  developed.  The  performance  of 
the  algorithm  and  limitations  of  the  architecture  have  been  simulated,  and  the  expected 
performance  is  evaluated. 

In  addition,  as  part  of  our  Phase  I  effort,  ACM  Systems  explored  the  possibility  of  using  the 
information  acquired  during  receive  adaptation  tosteer  transmit  beams  in  either  retrodirective 
or  directive  fashion.  Two  applications  in  particular  are  addressed:  Low  Probability  of^ntercept 
(LPI)  and  Anti-Jam  (A/J)  communications.  The  concepts  for  application  of  an  Adaptive  Array 
Processor  (AAP)  to  LPI  and  A/J  communications  is  presented,  and  the  expected  performance 
is  simulated  and  evaluated. 

The  goal  for  Phase  II  of  this  SBIR,  as  stated  in  the  original  solicitation,  is  "...  demonstrate 
hardware  implementation  of  small  adaptive  antennas."  ACM  Systems’  response  to  this  Phase  II 
goal  is  the  recently  submitted  proposal,  dated  January  7,  1991,  and  entitled  A  Multi-Purpose 
Adaptive  Array  Processor  Hardware  Development  for  Frequency  Hopped.  Direct  Sequence 
Spread  Spectrum  Radios.  This  proposal  addresses  and  exceeds  the  Government’s  expectations 
for  Phase  II  of  this  SBIR  program. 
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4.  BODY  OF  REPORT 

4.1  STATEMENT  OF  PROBLEM  STUDIED 

The  Phase  I  technical  objectives  were  stated  in  the  proposal  as  follows: 

The  goal  o f  Phase  /  is  to  determine  the  ability  of  property  restoral  algorithms  to  be  used  for  blind 
signal  extraction  and  d  irective /retrodirective  transmission,  and  to  develop  low-cost  architectures  for 
performing  this  processing.  This  goal  is  accomplished  using  the  following  procedure. 

o  Identify  realistic  test  cases  for  developing  and  evaluating  the  blind  directive  / 
retrodirective  communication  system. 

o  Design  a  candidate  property  restoral  algorithm  that  is  best  able  to  blindly  adapt  the 
receive-path  array  weights  in  the  environment  identified  above,  and  develop  methods  for  using  the 
receive-path  weights  to  compute  the  transmit-path  array  weights  that  maximizes  the  transmitted 
signal  energy  in  the  direction  of  the  received  signal  communicator. 

o  Develop  a  realistic  architecture  for  performing  the  signal  reception  and  directive  / 
retrodirective  transmission  using  the  algorithm  designed  above,  and  determine  the  cost  / 
performance  tradeoffs  for  each  of  these  structures  under  realistic  dynamic  range  /  noise  conditions. 

o  Evaluate  the  performance  of  the  processing  algorithm  and  reception  /  transmission  system 
developed  above,  using  theoretical  analysis  and  computer  simulation. 

Two  primary  approaches  are  to  be  considered  for  performing  the  blind  emitter  reception:  LSCMA 
adaptation,  with  initialization  /  recapture  via  transience  exploitation;  and  spectral  self -coherence 
restoral. 

Rather  than  proceeding  blindly  towards  the  stated  objectives  of  the  proposal,  ACM  Systems 
contacted  several  interested  Government  technical  representatives  early  in  the  Phase  I  program, 
to  ascertain  the  areas  of  greatest  interest  for  investigation  and  development  under  this  contract. 
The  stated  objectives  of  the  proposal  were  widely  accepted  and  encouraged,  with  the  following 
points  of  emphasis  suggested  as  well: 

o  Any  receive  adaptation  algorithm  would  need  to  be  optimized  for  application  to 
spread  spectrum  communication  systems,  both  frequency  hopped  and  direct  sequence. 

o  Previous  Government  studies  and  independent  research  efforts  have  yielded  many 
receive  adaptation  algorithms.  An  effort  which  addresses  architectural  implementations  of 
these  various  algorithms,  as  well  future  algorithms,  was  needed. 

o  The  architectures  developed  as  part  of  this  study  must  consider  the  intended 
application:  existing  Army  tactical  ground  based  communication  systems.  This  implies  low  cost, 
small  size,  low  power  dissipation,  and  of  course  frequency  hopping,  direct  sequence  spread 
spectrum  signal  compatibility. 

The  objectives  as  stated  in  the  proposal,  modified  to  emphasize  the  areas  suggested  by 
concerned  technical  representatives  from  the  Government,  were  accomplished  via  a  Phase  I 
Work  Plan  which  included  the  following  tasks:  I)  Algorithm  Development,  2)  Architecture 
Development,  3)  Performance  Evaluation  and  4)  Documentation.  Detailed  descriptions  of  the 
results  of  this  study  are  organized  according  to  these  tasks  within  section  4.3. 
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4.2  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 

The  most  important  results  obtained  in  the  course  of  this  Phase  I  study  are  stated  briefly  below. 
Detailed  descriptions  and  explanations  are  provided  in  section  4.3. 

o  A  realistic  computer  model  of  jammers,  signals  of  interest,  AAP  front  end  processing, 
dynamic  range  constraints,  the  antenna  array,  directions  of  arrival,  etc.  has  been  established. 
Many  key  parameters  are  selectable  variables,  thereby  allowing  flexibility  to  study  a  variety 
of  factors  associated  with  each  of  these  topics. 

o  A  new  receive  adaptation  algorithm  has  been  developed  which  is  tailored  for 
frequency  hopped  applications,  the  Dominant  Mode  Block  Update  algorithm.  A  derivation  is 
provided,  along  with  formal  analysis  of  the  algorithm’s  convergence,  stability,  and  steady  state 
performance.  A  proof  indicating  this  algorithm  to  be  a  maximum  likelihood  estimator  was 
begun  and  will  be  submitted  for  publication  at  a  later  date.  Acknowledgement  to  this  project 
for  its  contribution  to  the  paper  will  be  made  at  that  time. 

o  The  algorithm  was  modeled  on  the  computer  and  its  performance  was  simulated  for 
various  signal,  antenna  array,  and  timing  conditions.  The  simulated  performance  results  are 
encouraging: 


Frequency  hopped  performance  is  verified,  with  convergence  the  maximum 
attainable  SINR  within  a  fraction  of  a  timeslot  (EPLRS  timing  assumed). 

t>  The  antenna  patterns  were  generated,  revealing  simultaneous  nulls  in  the 
direction  of  jammers  and  a  beam  in  the  direction  of  the  desired  signal. 

t>  The  number  of  host  radio  data  samples  required  for  integration  by  the  AAP 
is  investigated  by  simulation;  16  is  sufficient. 

>  The  impact  of  a  desired  signal  not  being  present  during  a  timeslot  in  which  it 
is  expected  was  explored;  there  is  little  degraded  performance,  and  is  easily 
adjusted  for  in  the  architecture.  This  is  likely  in  network  communications. 

t>  The  algorithm  converges  to  the  maximum  SINR  in  response  to  changing  desired 
signal  directions  of  arrival  within  a  fraction  of  a  timeslot. 

>  Several  architectural  limitations  were  identified  and  simulated;  most  proved 
not  to  impact  AAP  performance,  while  one  result  caused  redesign  of  the  off  time 
data  collection  scheme. 

>  Performance  for  varying  jammer  conditions  were  simulated,  including 
directions  of  arrival,  number  of  simultaneous  jammers,  jammer  power  levels,  and 
types  of  jammer  waveforms.  The  AAP  perf ormance  was  successf ul  and  consistent 
for  all  conditions  studied. 

»  Antenna  array  parameters  were  investigated  for  their  impact  on  AAP 
performance,  including  the  number  of  elements,  the  geometry  of  the  array,  and 
the  spacing  of  the  elements. 
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4.?  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS,  (continued) 

o  A  receive  AAP  architecture  was  developed  which  is  compatible  with  existing  Army 
tactical  radios  and  is  capable  of  supporting  the  Dominant  Mode  Block  Update  algorithm  as  well 
as  many  others.  The  features  of  this  architecture  include: 

t>  Complex  weights  can  be  implemented  at  RF  if  required  by  the  host  Army 
radio,  but  additional  flexibility  is  available  if  IF  may  be  used.  Either 
implementation  is  supported  by  the  architecture. 

t>  The  statistical  measurements  required  for  adaptation  are  made  directly  at  IF 
for  on  time  processing,  and  at  baseband  for  off  time  data  collection.  This 
combination  is  optimal  for  fast  adaptation  while  allowing  a  variety  of  algorithms 
to  be  implemented. 

>  The  weights  are  calculated  digitally,  allowing  maximum  flexibility  in 
algorithm  selection.  This  also  allows  for  storage  and  retrieval  of  weights,  thereby 
greatly  enhancing  performance  in  frequency  hopped  systems. 

>  Techniques  for  correcting  or  calibrating  all  significant  sources  of  error  are 
identified. 

»  Options  for  all  circuits  required  for  hardware  development  of  the  AAP  are 
identified  and  discussed. 

o  Methods  for  generating  antenna  beams  for  transmission  were  developed.  Both 
directive  and  retrodirective  techniques  were  investigated.  Of  particular  interest  is  applications 
to  Low  Probability  of  Intercept  (LPI)  and  Anti-Jam  (A/J)  communications.  Techniques  for 
smoothing  aperature  estimates,  proven  by  simulation  to  be  necessary,  are  suggested. 

o  Both  LPI  and  A/J  scenarios  are  simulated  and  evaluated,  showing  significant  system 
performance  improvements.  Antenna  patterns  and  areas  exceeding  critical  power  levels  in  the 
vicinity  of  the  transmitter  are  generated  to  support  the  evaluations. 


4.3  WORK  CARRIED  OUT  /  RESULTS  OBTAINED 
4.3.1  INTRODUCTION 

The  details  of  the  work  carried  out  for  this  Phase  I SBIR  and  the  results  obtained  are  presented 
in  this  section.  The  work  is  divided  into  four  primary  efforts:  Algorithm  Development, 
Architecture  Development,  Computer  Modeling  and  Simulations,  and  Documentation.  The  work 
carried  out  and  the  results  obtained  for  each  of  these  efforts  are  included  in  the  following 
sections. 
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4.3.2  ALGORITHM  DEVELOPMENT 
4.3.2. 1  INTRODUCTION 

Receive  adaptation  algorithms  have  been  customized  for  a  variety  of  applications.  An  Army 
mobile  tactical  ground  based  communication  system  requires  that  the  algorithm  include: 

o  Compatibility  with  spread  spectrum  techniques,  both  frequency  hopping  (FH)  and 
direct  sequence  (DSSS).  Most  Army  radios  employ  FH,  some  use  DSSS  techniques;  these  are 
definitely  the  trends  for  future  radios. 

o  Fast  adaptation  times  are  desired  for  compatibility  with  FH  systems,  preferably 
providing  convergence  in  a  fraction  of  the  transmit  time  period  of  the  network  radios.  If  this 
is  not  achieved,  the  adaptation  must  occur  over  several  timeslots.  In  a  network  communication 
system,  the  radio  attempts  to  receive  signals  from  several  other  radios  in  the  network,  many 
having  different  directions  of  arrival  (DOAs).  In  this  case,  steering  a  beam  as  well  as  nulls  over 
several  timeslots  is  difficult  since  the  DOA  of  the  desired  signal  changes  during  adaptation. 

o  Steering  nulls  in  the  directions  of  interferes  is  required.  Reduced  interference  output 
power  will  generally  improve  the  SINR  provided  to  the  host  radio,  but  not  necessarily.  The 
desired  signal  may  have  a  DOA  which  is  in  the  direction  of  an  inadvertant  null  formed  by 
minimizing  the  interference  output  power.  Simultaneously  steering  a  beam  in  the  desired 
signal’s  DOA  assures  the  SINR  will  be  improved,  and  is  therefore  highly  desirable. 

o  Ease  of  implementation,  thereby  resulting  in  a  design  which  minimizes  the  amount 
of  hardware  required.  This  requirement  reflects  the  need  for  the  AAP  to  be  small  and  mobile. 

o  Tolerance  to  the  severe  multi-path  conditions  associated  with  ground  communications. 

o  Frequency  compensation  or  memory  of  previous,  frequency  dependent  weights, 
necessary  to  support  the  wide  operating  bandwidths  of  the  FH  systems. 

o  Capability  to  support  the  wide  cancellation  bandwidths  required  for  DSSS  signals. 

o  Capable  of  performing  blind  adaptation  -  without  special  purpose  antennas,  special 
sensor  geometries,  matched  receiver  to  array  calibrations,  or  special  training  signals. 

o  Optimized  for  operation  in  severe  jamming  environments. 

A  brief  summary  of  the  receive  adaptation  algorithms  considered  as  part  of  this  study  is 
provided  in  this  section.  A  baseline  algorithm  is  selected  for  focus  in  the  architecture 
development  and  computer  modeling  and  performance  simulations.  This  selection  is  justified 
and  further  developed  to  meet  the  specific  requirements  of  this  application.  An  overview  of 
the  selected  algorithm  is  provided,  followed  by  its  derivation,  an  analysis  of  its  convergence 
and  stability  properties,  and  an  analysis  of  its  expected  performance. 

The  concepts  of  directive  and  retrodirective  communications  are  addressed  as  well.  The 
algorithms  for  generating  the  transmit  weights  based  on  receive  adaptation  information  is 
described,  with  emphasis  on  two  applications:  Low  Probability  of  Intercept  (LPI) 
communications  and  purely  anti-jam  (A/J)  capabilities.  Figures  of  merit  for  evaluating  the 
directive  communication  system’s  performance  are  established,  and  the  optimal  transmit 
weights  for  both  cases  of  interest  are  calculated. 
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4. 3.2.2  RECEIVE  ALGORITHM  INVESTIGATION 

Several  approaches  were  considered  in  this  Study  for  adapting  the  receive-site  antenna  array, 
based  on  the  possible  attributes  of  the  modulation  format  to  be  employed  by  the  communication 
link.  These  algorithms  can  be  grouped  into  five  basic  classes  of  techniques: 

o  nonblind  techniques,  which  adapt  the  processor  using  a  nonblind  steepest-descent  or  least- 
squares  acquisition  algorithm  [1],  followed  by  a  post-acquisition  tracking  algorithm  such 
as  a  decision-direct  [2]  or  pilot-directed  [3,4]  algorithm; 

°  coherence-exploiting  techniques,  such  as  the  Griffiths  P-vector  and  Frost  algorithms  [5,6], 
tne  self-coherence  restoral  (SCORE)  algorithms  [7,8],  and  higher-order  cumulant  and 
spectral  analysis  techniques  [9,10],  which  adapt  the  processor  to  exploit  or  restore  known 
spatial,  spectral,  conjugate,  or  higher-order  self-coherence  of  the  transmitted  signal 
wavef  orms; 

°  modulus-restoral  (MORE)  techniques,  such  as  the  constant-modulus  algorithm  (CMA) 
[1 1,12],  known  modulus  algorithm  [13],  multiple-modulus  algorithm  [14],  and  adaptive 
modulus  algorithm  [7],  which  adapt  the  processor  to  restore  known  or  estimable  modulus 
properties  of  the  transmitted  signal  waveform  to  the  processor  output  signal; 

°  regenerative  techniques,  such  as  the  regenerative  DSSS  acquisition  algorithms  developed 
in  [15,16]  and  the  more  general  reference-directed  adaptation  techniques  developed  in 
[17,18,19],  all  of  which  exploit  the  known  spectral,  temporal  or  statistical  structure  of 
the  transmitted  signal  to  generate  a  reference  signal  for  us  in  a  conventional  steepest- 
descent  or  least-squares  adaptation  algorithm;  and 

°  Maximin  techniques,  such  as  the  steepest-asccnt  Maximin  algorithms  developed  in 
[20,21,22]  and  the  dominant-mode  Maximin  algorithm  developed  in  [23,24]  (referred  to 
as  the  burst  acquisition  algorithm  in  those  reports),  and  the  adaptive  detection  and  DSSS 
acquisition  algorithms  developed  in  [25,26]  (which  can  be  interpreted  as  rank-1  Maximin 
algorithms),  all  of  which  adapt  the  processor  to  restore  known  or  estimable  temporal, 
spectral,  or  spreading  support  of  the  transmitted  signal  to  the  processor  output  signal. 

However,  it  is  possible  to  quickly  narrow  these  techniques  to  a  few  obvious  choices  after 
considering  the  stringent  communication  system  requirements  described  in  the  previous  section. 
The  requirement  that  the  communication  link  be  secure  and  LPI  immediately  precludes 
nonblind  techniques,  which  require  transmission  of  a  know  (and  therefore  easily  detectable) 
preamble  signal  in  order  to  acquire  the  communication  signal  at  the  beginning  of  the 
transmission.  The  requirement  that  the  techniques,  since  the  preamble  signal  must  then  be 
retransmitted  at  the  beginning  of  every  frequency  hop  to  track  the  frequency-dependent 
variation  of  the  dehopped  transmission  channel. 

These  requirements  also  preclude  a  number  of  the  blind  adaptation  algorithms.  The  SCORE 
algorithms  are  precluded  for  LPI  modulation  formats,  since  these  formats  must  exhibit  a  very 
low  degree  of  spectral  or  conjugate  self  -coherence  to  prevent  detection  by  conventional  f  eature 
detection  technique.  In  particular,  the  convergence  time  of  the  SCORE  algorithms  would 
exceed  the  hop  duration  of  most  practical  frequency-agile  communication  waveforms  at  the 
receiver  SINRs  likely  to  be  encountered  by  such  systems.  Higher  order  spectral  analysis 
techniques  are  inapplicable  to  this  communication  waveform  for  similar  reasons.  The 
requirement  that  the  processor  operate  without  knowledge  of  the  transmit  signal  DOA,  special 
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4. 3.2.2  RECEIVE  ALGORITHM  INVESTIGATION  (continued) 

array  geometries,  or  array  calibration  data  precludes  the  remaining  coherence  exploiting 
techniques,  ali  of  which  rely  on  this  information. 

The  additional  requirement  that  the  AAP  be  small  and  low-cost  mitigates  against  the  majority 
of  modulus  restoral  techniques.  MORE  techniques  have  been  developed  that  can  adapt  quickly 
enough  to  track  frequency  agile  communication  signals  [27,28],  and  to  operate  in  the  presence 
of  multiple  interferers  [23,29,30].  However,  this  capability  comes  at  the  cost  of  increased 
computational  complexity.  In  particular,  the  multitarget  MORE  rlgorithm  described  in 
[23,29,30]  requires  (and  achieves)  simultaneous  capture  and  separation  of  all  of  the  signals 
received  at  the  transmit  signal  frequency.  This  capability  is  excessive  in  the  radios  of  interest 
here1,  which  only  require  extraction  of  the  intended  transmit  signal  form  the  receive 
environment. 

The  regenerative  and  Maximin  techniques  can  provide  a  simpler  means  for  satisfying  these 
requirements  if  the  communication  signal  has  a  DSSS  and/or  frequency  agile  modulation 
format.  All  of  these  techniques  a.e  extremely  robust,  requiring  only  knowledge  of  the 
spreading  code  or  hopping  sequence  to  train  the  receiver  processor2.  In  addition,  the 
techniques  are  inherently  selective  providing  only  the  intended  transmit  signal  at  the  processor 
output. 

Of  these  techniques,  the  Maximin  approach  is  superior  from  robustness,  performance  and 
implementation  viewpoints,  particularly  if  the  communication  signal  has  a  frequency  agile 
communication  format.  The  Maximin  approach  assumes  no  knowledge  of  the  baseband  signal 
structure,  and  can  operate  effectively  if  the  baseband  signal  has  no  discernable  (or  exploitable) 
structure,  e.g.,  if  the  baseband  signal  is  a  stationary  Gaussian-disuibuted  random  process. 
Consequently,  the  approach  is  inherently  less  sensitive  to  error  in  the  baseband  signal  model, 
for  example,  if  the  transmit  signal  is  subjected  to  over  modulation,  bandlimiting  or  clipping 
operations  over  the  communication  channel.  Moreover,  the  Maximin  approach  is  insensitive 
to  many  commonly  occurring  errors  in  the  spreading  format  used  to  modulate  the  baseband 
signal,  such  as  sub-hop  timing  error  or  frequency  offset  in  FH  spreading  formats  [23,24]. 

In  fact,  the  Maximin  approach  can  usually  outperform  regenerative  techniques  under  even 
ideal  conditions  where  the  transmitted  signal  waveform  is  not  subject  to  modelling  error.  The 
ability  of  Maximin-adapted  antenna  arrays  to  capture  frequency  agile  signals  with  newly  rhw 
maximum-attainable  signal-to-noise  ration  (SINR)  of  the  array  has  been  demonstrated  via  both 
theoretical  analysis  [23,24]  and  computer  simulation  [20]-[24],  and  for  exceptionally  sever 
interference  powers  and  reception  geometries.  No  similarly  powerful  result  has  been  reported 
in  the  literature  for  any  of  the  regenerative  techniques  developed  to  date. 


’However,  this  capability  may  be  of  interest  in  some  specialized  applications,  e.g,  if  the 
radio  is  a  transponder  or  central  node  in  a  combat  net  radio  system. 

2The  regenerative  techniques  reported  in  [  1 5]-[  1 9]  also  employ  nonlinearities  that  all  (at 
least  implicitly)  assume  additional  information  about  transmitted  signal  baseband,  for  example, 
that  it  has  a  low  modulus  variation  or  a  known  modulation  format. 
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4. 3.2.2  RECEIVE  ALGORITHM  INVESTIGATION  (continued) 

Moreover,  it  can  be  shown  [31]  that  the  Maximin  approach  can  be  derived  from  the  maximum- 
itkelihood  (ML)  estimate  of  the  baseband  (dehopped)  data  signal,  given  simple  regularizing 
conditions  on  the  transmitted  and  received  signal  waveforms.  Specifically,  if  the  received 
signal  data  signal  x(n)  can  be  modelled  by 

x(n)  =  as(n )  +  i(n) 

over  discrete  time  interval  { 1  ,...,N),  where  i(n)  is  ^  white  circularly-symmetric  complex  Gaussian 
(CSCG)  interference  waveform  with  mean  zero  and  unknown  autocorrelation  matrix  R^,  and 
where  s(n)  has  limited  statistical  support,  such  that  a  linear  thinning  operator  Pg  exists  that 
satisfies 


Ps{s(n))  =  s(ri) 


over  { 1  ,...,N},  then  the  joint  ML  estimate  of  s(n),  a,  and  Ry  is  given  by 


full")  =  wuLxpaJn) 

auL  =  /  iM. 

=  Rxx  -  6 UL^M^l^S, 

XpaJn)  =  P,(x(n)) 


*4tcsxnj*ct  ML 


MV*  ML 


X reject^")  =  *(»)  ~  *n o»(«> 


where  (  )H  denotes  the  conjugate-transpose  (Hermitian)  operation  and  <>N  denotes  time 
averaging  over  the  reception  interval,  and  where  wML  is  the  maximal  point  of  the  Maximin 
objective  function 


FulW  = 


*  t&ct*  Ht*a 


w 


The  thinning  operator  P,  can  be  interpreted  as  a  generalized  bandpass  operation  that  passes  the 

region  of  statistical  support  containing  s(n)  and  rejects  the  remaining  data  in  the  environment. 

Example  thinning  operators  include: 

°  bandiimiting  operators ,  applicable  to  environments  where  s(n)  has  limited  spectral 

support,  i.c.,  where  s(n)  exists  over  a  narrow  frequency  band  within  the  overall  received 
data  passband  (for  example,  if  s(n)  is  a  CW  tone); 

°  time-limiting  operators,  applicable  to  environments  where  s(n)  has  limited  temporal 

support,  i.e.,  where  s(n)  exists  over  a  narrow  subset  of  the  overall  reception  interval  (for 
example,  if  s<n)  is  a  burst,  pulsed,  or  time-hopped  waveform);  and 
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4.3.2. 2  RECEIVE  ALGORITHM  INVESTIGATION  (continued) 

o  dehop-rchop  operators,  applicable  to  environments  where  s(n)  has  limited  and  time- 
varying  spectral  support,  i.e,  where  s(n)  exists  over  a  narrow  "instantaneous”  frequency 
band  that  hops  over  the  receive  data  passband  during  the  reception  interval. 

Other  thinning  operators  that  can  be  postulated  for  different  signal  modulation  formats 
include:  despread-respread  operators ,  applicable  to  DSSS  signals;  and  generalized  despread- 
respread  operators,  applicable  to  hybrid  FH/DS  modulation  formats. 

Given  this  derivation,  the  following  four-step  procedure  can  be  used  to  determine  the  ML 
estimate  of  the  transmitted  signal  waveform: 

1.  pass  each  element  of  x(n)  through  the  generalized  bandpass  and  complementary 
bandreject  filter  that  is  matched  to  s(n),  to  form  thinned  signa-plus-interference  data 
signal  xpM,(n)  and  signal-free  data  signal  xreject(n); 

2.  compute  the  sample  autocorrelation  of  xpM,(n)  and  xreject(n)  over  the  reception  interval; 

3.  compute  the  maximal  value  of  the  Maximin  objective  function  FML(w)  defined  above; 

and  C., 

4.  compute  the  transmitted  signal  estimate  form  the  thinned  data  signal  xpMB(n),  using  the 
formula  given  above. 

From  this  procedure,  it  can  be  seen  that  the  Maximin  algorithms  reported  in  [20,21,22]  provide 
the  ML  estimate  of  an  FH  signal  using  a  stecpest-desccnt  algorithm3.  Similarly,  it  can  be  seen 
that  the  burst-acquisition  algorithm  reported  in  [23,24]  provides  the  estimate  of  a  burst  signal4 
using  the  exact  solution  to  Fj^w),  given  by  the  dominant  mode  (eigenvector  with  maximum 
cgenvalue)  of  the  Maximin  eigenequation 

z  w  -  Rx x  w 

With  some  additional  effort,  the  adaptive  detection  algorithm  reported  in  [25]  and  the  DSSS 
acquisition  algorithm  reported  in  [26]  can  also  be  interpreted  as  rank-1  approximations  to  the 
ML  estimator  of  a  burst  or  DSSS  signal,  where  the  statistical  support  for  xpaM(n)  is  so  sparse 
that  the  autocorrelation  of  xpBIi(n)  has  rank  1.  In  this  environment,  Kelly  and  Dlugos  have 
independently  derived  statistical  models  for  s^n)  and  FML(wMI/)  based  on  Wishart 

distribution  of  the  ACMs  of  xpMI(n)  and  x^^n). 


3More  exactly,  the  Maximin  algorithms  provide  a  close  approximation  to  the  ML  estimate 
of  the  FH  signal,  since  the  algorithm  computes  x^^n)  using  narrowband  and  causal  processor 
operations,  and  updates  w  using  a  stochastic  steepest-ascent  algorithm. 

4More  exactly,  the  burst  acquisition  algorithm  provides  a  close  approximation  to  the  ML 
estimate  of  the  FH  signal,  since  the  algorithm  estimates  xpMI(n)  and  x^j^n)  using  the  optimal 
value  of  Fj^fw)  as  an  ON-time  detection  statistic,  and  almost  certainly  estimates  wML  using  an 
estimator  with  significant  timing  error. 
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4.3. 2. 2  RECEIVE  ALGORITHM  INVESTIGATION  (continued) 

This  result  has  an  important  bearing  on  the  performance  of  the  Maximin  approach.  The  strong 
link  between  this  approach  and  the  ML  estimator  indicates  that  the  Maximin  approach  should 
provide  the  best  possible  estimate  of  the  transmitted  signal  waveform  over  a  given  reception 
interval5.  Consequents,  the  Maximin  approach  should  provide  faster  convergence  than  the 
regenerative  approaches  in  environments  where  the  structure  of  the  transmitted  signal  is 
completely  characterized  by  its  statistical  support6.  This  prediction  is  borne  out  by  the 
computer  simulations  conducted  to  date  for  all  of  these  techniques. 

The  analytical  tractability  of  the  Maximin  approach  (relative  to  the  regenerative  approaches) 
also  has  an  important  bearing  on  the  reliability  of  these  techniques.  The  ability  to  analyze  the 
convergent  behavior  of  the  adaptation  algorithm  prior  to  actual  hardware,  by  allowing  the 
system  engineers  to  predict  system  performance  for  a  variety  of  errors  scenarios,  and  to  design 
reliable  safeguards  against  any  failure  modes  detected  during  this  analysis.  This  design 
philosophy  is  demonstrated  in  Section  4. 3. 2.6,  where  performance  predictions  are  made  for  a 
number  of  error  scenarios  in  addition  to  the  nominal  reception  scenario. 

This  analytical  tractability  also  extends  to  the  convergence  performance  of  certain  of  the 
Maximin  algorithms.  In  particular,  the  convergence  speed  of  the  dominant-mode  Maximin 
algorithm  described  in  [23,24]  and  developed  on  this  Study  can  be  expressed  in  terms  of  the 
eigenstructure  of  the  Maximin  eigenequation,  and  related  to  key  environment  parameters  such 
as  the  maximum-attainable  SINR  of  the  array.  This  analysis  is  demonstrated  in  Section  4. 3.2.5- 
4. 3.2. 6,  where  it  provides  a  powerful  estimate  of  the  convergence  speed  of  the  processor,  as  well 
as  important  information  for  further  improving  the  performance  of  the  processor.  In 
particular,  the  modified  dominant-mode  Maximin  algorithm  described  in  Section4.3.2.5  is  a 
direct  result  of  this  analysis. 

In  addition  to  the  performance  advantages  described  above,  the  Maximin  algorithm  possesses 
important  implementation  advantages  over  the  regenerative  techniques  developed  to  date. 
Almost  all  of  the  regenerative  techniques  developed  to  date  impose  a  nonlinearity  of  some  sort 
on  the  beamformer  output  signal  during  the  adaptation  process7.  This  nonlinearity  complicates 


5Of  course  a  better  estimator  can  be  devised  if  additional  information  about  the  transmitted 
signal  is  incorporated  into  the  received  data  model.  For  example,  in  [31]  it  is  shown  that  the 
ML  estimator  reduces  to  a  hybrid  Maximin/MORE  algorithm  if  additional  knowledge  of  the 
transmitted  signal  modulus  is  incorporated  into  the  received  data  model.  In  fact,  this  processor 
can  bear  a  close  resemblance  to  the  regenerative  processors  for  some  transmitted  signal  types. 

6However,  the  regenerative  techniques  may  have  some  performance  advantages(under  ideal 
conditions)  if  additional  information  about  the  transmitted  signal  is  incorporated  into  the 
received  data  model.  Even  then,  the  robustness  and  implementation  advantages  of  the  Maximin 
approach  will  mitigate  against  the  regenerative  techniques  in  the  application  of  interest  here. 

7  A  notable  exception  is  [18]  which  expressly  removes  the  limiting  operation  from  the 
regenerative  processor  in  order  to  remove  the  undesirable  effects  imposed  by  this  limiter. 
However,  the  resultant  processor  structure  is  still  more  complex  than  the  equivalent  Maximin 
processor  structure,  and  does  not  converge  to  the  maximum-attainable  SINR  solution  in  most 
environments. 
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4. 3. 2. 2  RECEIVE  ALGORITHM  INVESTIGATION  (continued) 

the  implementation  of  the  resultant  processor,  and  can  result  in  hysteresis  effects  such  as 
oscillatory  weight  instability  in  the  processor  weights  [17,32]. 

In  contrast,  the  Maximin  algorithms  do  not  impose  a  nonlinearity  on  the  beamformer  output 
signal,  and  are  consequently  not  subject  to  the  hysteresis  effects  experienced  by  the 
regenerative  techniques.  In  addition,  the  computational  complexity  of  these  techniques  can  be 
much  lower  than  the  complexity  of  any  of  the  other  candidate  algorithms  considered  on  this 
Study(including  the  regenerative  techniques). 

Most  notably,  the  computational  complexity  of  the  dominant  mode  Maximin  algorithm  is  lower 
than  any  of  the  "rapidly-converging"  algorithms  considered  on  this  Study,  including  the  cross- 
SCORE  and  least-squares  MORE  algorithms,  the  least-square  extensions  of  the  regenerative 
techniques  developed  to  date,  and  the  nonblind  least-squares  algorithm.  This  complexity  savings 
results  from  two  sources: 

°  the  dominant  mode  Maximin  algorithm  does  not  require  cross-correlation  data  to 
operate,  only  autocorrelation  data  (which  is  being  collected  by  other  algorithms  as  well 
as  the  Maximin  algorithm);  and 

o  the  signal-free  autocorrelation  matrix  is  only  inverted  once  per  hop)  rather  than  every 
block  as  is  required  in  all  of  the  other  algorithms). 

Together  with  the  performance  and  robustness  features  of  the  Maximin  algorithms,  these 
implementation  features  provide  the  Maximin  algorithms  (and  particularly  the  dominant  mode 
Maximin  algorithm)  with  formidable  advantages  over  the  competing  approaches. 


4. 3. 2. 3  BASELINE:  DOMINANT  MODE  MAXIMIN  ALGORITHM 

The  history  of  the  Dominant  Mode  (Block  Update)  Maximin  algorithm  was  presented  in  the 
previous  section,  and  its  choice  as  the  baseline  algorithm  for  this  study  was  justified.  The  basic 
operation  of  the  algorithm  is  explained  in  this  section  using  three  steps:  1)  an  introduction  to 
dominant  mode  eigenvectors,  2)  a  description  of  the  power  method  of  determining  eigenvectors, 
and  3)  a  discussion  of  how  the  eigenvector  estimates  are  updated  using  blocks  of  stochastic 
data. 

Dominant  Mode  Eigenvectors  - 

The  Dominant  Mode  Block  Update  algorithm  attempts  to  maximize  the  SINR  at  the  output  of 
the  AAP.  The  relevant  equation  is  F(w)  =  SINR  +  1  ■  whRjw/whR0w.  Obviously,  maximizing 
this  function  also  maximizes  the  SINR.  The  weight  vector  which  gives  the  maximum  SINR  is 
the  solution  (A.max  w„,.,)  to  the  generalized  eigenequation  A.R0w  =  Rtw.  It  is  proved  in  section 

4. 3.2.4  that  the  eigenvector  solution  to  this  eigenequation  with  the  largest  eigenvalue  (the 
eigenequation’s  dominant  mode)  maximizes  the  function  F(w).  This  vector  gives  the 
complex  weights  which  should  be  applied  to  the  input  signals  from  the  antenna  array  elements 
to  provide  the  maximum  output  SINR.  Note:  Complete  descriptions  of  these  variables  and  the 
derivations  are  presented  in  section  4. 3.2.4. 
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4.3. 2.3  BASELINE:  DOMINANT  MODE  BLOCK  UPDATE  (continued) 

Power  Method  - 

There  are  many  ways  to  determine  the  eigenvalues  of  a  matrix.  One  method  that  is  particularly 
easy  to  implement  and  will  converge  quickly  if  there  is^one  dominant  eigenvalue  is  the  power 
method  (33],  It  starts  with  an  initial  guess  w0  and  successively  calculates  better  estimates  by 
multiplying  the  last  estimate  by  the  same  original  matrix  A:  wk+1  =  Awk.  After  k  iterations,  the 
estimate  is  wk  =  A^Wq.  Assuming  A  has  a  full  set  of  eigenvectors  xk  the  vector  wk  is  given  by 
the  formula  for  a  difference  equation: 

wk  =  c1X1kXi  +  ...  +  cnAnkxn. 

If  the  eigenvalues  are  numbered  according  to  increasing  magnitude,  and  the  maximum 
eigenvalue  is  much  larger  than  any  of  the  others,  then 

wkAnk  =  Cl(  V*n)k*l  +  -  +  Cn-l(*n-lAn)kxn-l  +  CnXn. 

As  long  as  the  initial  guess  w0  has  some  component  of  xB,  so  that  cn  is  not  0,  the  cnxn  term  will 
dominate.  Since  A.n  is  larger  than  any  other  Ak,  the  terms  (Aj/A.n)  are  less  than  1;  as  the  number 
of  iterations  increase,  (A.1/2.n)k  approaches  0.  The  process  converges  to  provide  an  estimate  of 
w  which  is  proportional  to  the  maximum  eigenvector  xn.  A  scaling  factor  must  be  introduced 
in  the  iterations  to  avoid  uk+1  from  growing  very  large  or  very  small;  uk  should  be  divided  by 
its  first  component  ukl  before  taking  the  next  step.  This  results  in  a  power  method  iteration 
which  gives  uk+1  =  Auk/ukl,  which  has  uk+1  converging  to  and  ukl  converging  to  *mBX. 

The  power  method  generates  A*;  when  A  -  R0_1Rjas  defined  above,  Ak  approaches  a  diagonal 
matrix  with  the  first  diagonal  element  approaching  the  SINR  +  1,  and  all  other  diagonal 
elements  approaching  1  in  nominal  communication  environments.  For  large  values  of  maximum 
attainable  SINR  this  procedure  should  converge  very  quickly.  Conditions  may  arise  in  which 
the  eigenvalue  spreads  are  not  conducive  to  fast  convergence  in  some  cases  (see  section  4.3.2.6 
for  transient  interference  "collision"  analysis).  To  speed  the  convergence  in  a  transient 
interference  environment,  or  for  other  such  conditions,  a  factor  u  is  introduced  as  follows:  A 
=  Ro^Rj  -  ul.  In  this  case,  Ak  approaches  a  diagonal  matrix  with  the  first  diagonal  element 
approaching  the  SINR  -  u,  and  all  other  diagonal  elements  approaching  1  -  u.  Selection  of  u  can 
be  customized  for  a  specific  jammer  threat  or  conditions.  Under  normal  conditions  the 
algorithm  may  converge  slightly  more  quickly  than  with  u  *  0;  for  certain  special  cases  the 
algorithm  will  converge  much  more  quickly.  Further  discussion  of  convergence  time  is  provided 
in  section  4. 3.2.5. 

Stochastic  Block  Updates  - 

In  the  actual  implementation  of  this  algorithm,  A  =  R0-1Rj;  data  f or  R0_1  is  collected  during  an 
entire  timeslot  which  is  different  than  the  timeslot  used  to  collect  data  for  Rr  Furthermore, 
the  data  for  Rk  is  collected  over  several  blocks  of  data  samples  within  one  timeslot.  The  data 
collection  process  is  described  in  detail  in  the  Architecture  section  of  this  report,  section  4.3.3. 
The  impact  to  the  power  method  is  that  the  matrix  A  changes  slightly  with  each  block  update, 
even  though  the  input  signals  have  not  changed  statistically  between  blocks.  If  the  integration 
time  for  each  block  is  long  enough,  the  deviation  between  sequential  A  matrices  will  be 
negligible.  Even  with  modest  integration  times  (16  samples),  the  affect  to  the  performance  of 
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4. 3. 2. 3  BASELINE:  DOMINANT  MODE  BLOCK  UPDATE  (continued) 

Stochastic  Block  Updates  (continued!  - 

the  algorithm  is  negligible  (verified  by  simulation,  see  section  4.3.4.).  This  is  because  multiple 
iterations  of  uk+1  =  Auk/u  smooths  the  variances  of  the  A  matrix  from  block  to  block.  The 
power  method  does  not  require  a  deterministic  A  matrix  without  change,  but  will  converge  with 
a  stochastic  A  matrix  taken  from  a  statistically  unchanging  RF  environment. 

4. 3. 2. 4  DOMINANT  MODE  DERIVATION 

The  objective  of  the  Dominant  Mode  Block  Update  algorithm  is  to  maximize  the  Signal  to 
Interference  plus  Noise  Ratio  (SINR)  at  the  output  of  the  AAP.  The  AAP  output  signal  is  y(t), 
and  is  equal  to  the  sum  of  the  weighted  input  signals.  Therefore, 

y  =  Ewk*xk(t)  =  wHx(t) 

where  ()H  denotes  Hermitian  (conjugate  transpose)  operation. 

The  AAP  output  power  may  be  written  as  Pout  =  <|y(t)|2>T  =  wH  <x(t)xH(t)>Tw  =  w^^’w 
where  <  >x  denotes  time  averaging  over  an  interval  of  duration  T  and  R^’  is  the  sample 
autocorrelation  of  x(t).  If  the  time  interval  is  long  and  x(t)  is  stationary  over  the  time  interval, 
then  R^’  can  be  approximated  by  its  expected  value  Rxx=  EfR^’J. 

Since  the  AAP  output  is  merely  a  weighted  summation  of  all  the  input  signals,  the  components 
of  the  AAP  output  depends  on  the  signals  present  at  the  input  when  the  measurements  are 
taken.  Let  S  represent  the  desired  signal,  or  Signal  Of  Interest  (SOI),  I  represent  the  total 
interference,  or  Signals  Not  Of  Interest  (SNOI),  and  N  represent  the  noise  generated  internally 
by  the  AAP  circuits.  In  an  RF  environment  with  active  jammers,  it  is  assumed  there  are  always 
interferes  present.  When  the  AAP  is  operational,  its  internal  noise  will  always  be  present.  The 
SOI,  however,  may  or  may  not  be  present.  Two  general  output  power  equations  may  be  written 
which  describe  these  conditions: 

pouto  *=  I  +  N.  and  Pouu  =  S  +  I  +  N. 

The  subscript  1  or  0  designates  whether  the  SOI  is  present  or  not.  Forming  the  ratio  Pouti/Pouto 
gives  (S  +  I  +  N)  /  (I  +  N),  which  is  equal  to  the  SINR  +  1.  Obviously,  the  weights  which 
maximize  SINR  +  1  will  result  in  a  maximum  SINR. 

The  presence  or  absence  of  signals  is  included  in  the  variable  x,  so  the  covariance  matrix  Rn 
will  depend  on  whether  the  SOI  is  present  or  not.  Representing  R^  as  R0  when  the  SOI  is  not 
present,  and  as  R,  when  the  SOI  is  present,  the  power  ratio  can  be  written  as 

F(w)  =  SINR  +  1  =  whRjw/whR0w. 

The  objective  of  the  adaptation  algorithm  is  to  maximize  this  ratio.  The  weight  vector  for 
which  this  function  is  a  maximum  can  be  designated  w„..„  and  has  the  following  solution: 

Vaxwm»x=  R0'lRlwn»r 

This  solution  can  be  found  by  setting  the  gradient  of  F(w)  equal  to  zero  and  solving  for  w.[20] 


Advanced 

CounterMeasure 

Systems 


16 


4. 3. 2. 4  DOMINANT  MODE  DERIVATION  (continued) 

The  same  solution  is  obtained  by  manipulating  F(w)  into  a  ratio  of  quadratics,  specifically  the 
Rayleigh’s  Quotient,  then  applying  Rayleigh’s  Principle  to  find  the  maximum  and  minimum 
eigenvalues  and  corresponding  eigenvectors.[33] 

Allowing  the  following  change  of  variables: 

z  =  R0^aw  and  A  =  Ro'^RjRq"1/2  (Note:  w  =  R0_1/Jz), 

F(w)  takes  the  form  of  the  Rayleigh’s  quotient  R(z): 

F(w)  =  whRjW/ wHR  0w 

=  (R0-1/,z)hR1(R0-1/sz)/(R0-1/2z)hR0(R0-1/*z) 

=  zhR0-1*,*R1R0-1/»z/zhR0-1/*R0R0-1/*z 

=  z^z/z^z  =  R(z). 

A  singular  value  decomposition  of  A  simplifies  the  quadratic  expressions  in  R(z).  The 
decomposition  sets  A  =  UAUH,  where  A  is  a  diagonal  matrix  whose  diagonal  elements  are  the 
eigenvalues  of  A  and  U  is  a  matrix  whose  columns  are  the  orthonormal  eigenvectors 
corresponding  to  those  eigenvalues.  The  corresponding  eigenequation  is  Xkuk  =  Aut  Applying 
this  singular  value  decomposition,  the  ratio  of  quadratics  R(z)  takes  the  form: 

R(p)  -  pHUHAUp/pHUHUp 

=  pHUHUAUHUp/pHUHUp 

«=  p^p/pHp 

-  (X^2  +  X2p22  +  ...  +  Xmpro2)/(pj2  +  p22  +  ...  +  pm2) 
where  z  ■  Up. 

The  Rayleigh  principle  states  that  the  Rayleigh  quotient  is  never  below  Xj  and  never  above  Xm; 
its  minimum  is  at  the  eigenvector  zx  and  its  maximum  is  at  zw  Intermediate  eigenvectors  are 
saddle  points.  The  eigenvalues  are  ordered:  zm>  ...  z  z2  z  z,. 

Working  backwards  to  return  to  the  original  variables,  the  equivalent  eigenequation  which 
results  is  XkR0wk  =  RjW^  These  transformations  require  that  R0  be  positive  definite  such  that 
R(j  2  exists  and  is  Hermitian.  Then,  the  maximum  value  of  F(w)  is  equal  to  Xm  and  is  achieved 
by  setting  w  proportional  to  vrm 

This  result  can  also  be  achieved  by  noting  that  the  non-Hermitian  matrix  B  *  Ro'^i  also  has 
a  singular  value  decomposition  given  by  WAA.  W  is  a  matrix  comprised  of  the  vectors  wk,  A  is 
a  matrix  comprised  of  the  vectors  a^  and  they  are  a  pair  of  reciprocal  basis  vectors  which 
satisf y  the  relationships  that  ak  =  Rqw^  wkHaj  =  1  when  k  =  1  and  wkHa,  =  0  when  k  *  I.  These 
results  imply  that  WHR0W  =  Im  and  =  A  The  diagonal  of  A  are  the  real  eigenvalues  Xk. 
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4. 3. 2. 4  DOMINANT  MODE  DERIVATION  (continued) 

Substituting  w  =  Wp  into  F(w)  gives 

F(w)  =  (Wp)HR1(Wp)/(Wp)HR0(Wp) 

=  p^W^^p/p^W^oW)? 

=  p^p/p^ 

which  has  an  unambiguous  maximum  at  p  =  ex,  yielding  =  Wj  and  Fmax  F(w„..J  =  A.r 

4.3. 2.5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS 

Note:  This  analysis  uses  slightly  different  notations  than  is  used  elsewhere  in  this  report. 

A  convergence  analysis  for  the  actual  (stochastic)  dominant-mode  algorithm  is  beyond 
the  scope  of  this  report.  However,  a  convergence  proof  is  obtainable  for  the  nonstochastic 
algorithm  that  maximizes  the  power  ratio 


TO 


&BR  w. 

&BQ  w 


(i) 


where  R  and  £  do  not  change  between  iterations  (JR  corresponds  to  the  on  time  autocorrelation 
matrix,  £  corresponds  to  the  off  time).  The  nonstochastic  update  procedure  is  as  follows: 


Initialization: 


d( o)  =  mo) 


(2) 


(3) 


Update 


k  -  k+1 : 


m  - 

d(0)  - 

u  =  L 


m 


(4) 

(5) 


(6) 


a  -  Ru 


(7) 


l£(*+l)  =  M  -  ni£(*) 


(8) 
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4.3.2.5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


i(fc+l)  =  y  -  v-dQc) 

(9) 

g(k+ 1)  = 

(10) 

v.  n  ^l) 

y&k+l)  =  - 

y/g(k+l) 

(ID 

_  /-s 

i  4 

II 

/—S 

1-H 

4 

(12) 

X()t+l)  =  )Ju.Hx 

(13) 

The  objective  of  the  algorithm  is  to  set  w(k)  equal  to  the  dominant  mode  of  the  eigen-equation 

A-Ctt  =  Ryl 

(14) 

as  k  -  without  instability  (|w(k)|  <  »  Vk).  In  particular,  the  algorithm  is  designed  to  yield 

*1 

m)  -  ~ 

&(k)  •*«!  X  y/IJ’ 

Uk)  -  A,j 

as  k  -  where  {Am,  w^,  arc  the  eigenvalues  and  corresponding  right-hand  and  left-hand 
eigenvectors  of  the  new  matrix  T  «= 


K*k  =  1*4 
*k*  =  f* 


***«/  =  «*-/ 


This  choice  of  normalization  is  made  to  ensure  that  the  ouptut  signal  power  is  normalized  to 
unity  over  the  signal  on  time  (wHRw). 


-  t  =  EHA  =  4 EH=  Lm 

E  =  bu,  •  •  •  jeJ 
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4.3.2.S  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


4  =  [«,  •• 


A  = 


*1 

0 


and  where  A.m  are  real  and  ordered,  with  distinct  maximum,  such  that  A.x  >  A2  2  •  •  •  2  A.M. 
As  shown  in  section  4.3. 2.4,  the  power  ratio,  F(w),  is  maximized  at  F(wmax)  =  Xj  for  wmax  =  Wj. 


The  convergence  proof  is  facilitated  by  proving  the  following  lemma: 
Lemma:  Given  (A,(k),  g(k),  w(k)}  in  the  equations  above. 


1  d(k)  =  R&k)  for  every  k  2  0 
2.  yvH(k)d(k)  =  1  for  every  k  2  0 


Proof  bv  induction: 

1.  Both  properties  are  true  fot  k*0 


#0) 


MQ) 

^"(0)4(0) 


4(0) 


4(0) 

vte*(0)4(0) 


=  g&O) 


-  *"(0)4(0)  =  1 

2.  Assume  both  properties  are  true  for  k,  then  show  they  must  be  true  for  k+1 
a)  Equations  (6)-(7)  and  property  1  of  the  lemma  yield 

u  =  g-'m 
=  g-lRm 
=  tm 

(15) 

y  =  ECld(k) 

=  Rim) 

(16) 


Advanced 

CounterMeasure 

Systems 


4.3.2.5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


20 


b)  Equations  (8)-(  12)  and  ( 1 5)-(  1 6)  yield 


=  «(*+l)[lf  -  infill,  a(£+l)  =  — 
=  a(k+  1)[Z  -  ii/]i&(*) 

d(k+ 1)  =  a(fc+l)[Y  -  \L&(k)] 

=  a(k+l)[gtyy(k)  ~  Ji  £*(*)] 

=  fia(*+i)iz  -  vnm) 

=  P-%k+ 1) 


(17) 


(18) 


c)  Equations  (1 1)-(12)  yield  condition  2  for  any  choice  of  w(k+l)  and  g(k+l),  as 
long  as  wH(k+l)  &(k+l)  *  0,  which  holds  for  any  positive  definite  R,  since  w 
and  q_  satisfy  condition  1.  “ 


Given  this  Lemma,  then  ( 1 7)-( 1 8)  can  be  used  to  express  w(k)  and  &(k)  directly  in  terms  of  w(0), 

(£  -  nz>*  m) 

u£k)  =  - —  — 

0)(r"  -  -  nD*  i£(0) 

(19) 

d(k)  =  RMk)  (20) 


-  £(*)  =  J&B(k-l) C'RC-'&ik-l) 


N 


i£g(0)(g"  -  y.tf-lfiRl(Z  -  tifl*-1  m) 

-  \xDk-lg(x  -  ix/)*-1  m) 


(21) 


Equations  (19)-(21)  can  be  expressed  in  terms  of  the  eigenvalues  and  eigenvectors  of  T,  by 
setting  “ 


using  the  identities 


m)  =  Km 

KHA  =  AWh=  I 


(22) 


(23) 


4  HtE  =  A  ;  4  =  QE 


(24) 


Advanced 

CounterMeasure 

Systems 


21 


4.3.2.S  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


-  AB(X-  nZ)ff  =  A  -  \il 


then  (19)-(21)  reduce  to 


EH&E  =  4 


E( 4  -  vD'ttt Q) 

^*(0)4(4  -  iiZ)2*  *(0) 

&*)  =  44cW 


*(*)  = 


N 


pa(0) 43(4  -  (iZ)"‘2  ^(0) 


£>)4  (4  -  n*)2*'2  *«>) 


,2Jk-2 


(25) 

(26) 

(27) 

(28) 

(29) 


Correct  convergence  is  demonstrated  by  fi(k)  -  a as  k  - 

Algorithm  convergence  can  be  measured8  using  the  generalized  correlation  coefficient 


S(k)  = 


1^)4^  I2 

[*"(*)£ 


(30) 


8  Note  that  stability  is  guaranteed  at  every  value  of  k,  since 

XBmRu(k)  -  1  V  * 
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4.3.2.5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


*1 

|j2*(0)(A  - 

1  £*(0)A(A  -  iiZ)S(O) 
_  X1q1  -  ^Ip.CQ)!2 

£*„(*»  ~  ^“W0)!2 


M 

1  +  E 

m-2 


a.  -  n'” 


,  M°)l  j  Ui  J(  Ai  -  •*  J 


-l 


(31) 


Equation  (31)  is  monotonically  convergent  to  unity  as  long  as  Ip^O))  >  0  (i.e.,  as  long  as  w(0) 
and  Wj  are  aox  orthogonal)  and 


< . 


A.j  +  A, 


Vm=2,...M 


>  p.  and  A.j  >  A,,  Vm>l 


(32) 


(33) 


The  speed  of  this  convergence  is  also  given  by  the  ratio 


r«mv  =  ““ 


>2-  M2 

lUi  -  . 

which  constitutes  the  largest  of  the  M-l  decay  constants  in  (31). 

Similarly,  the  eigenvalue  estimate  i(k)  also  converges  to  A.j  as  k  -  •»,  since 


m  = 


M 


tWoil’A.d.  -  n): 

_1 

M 


Ji-2 


Eb.«»l2U*.  -  M): 


2k- 2 


(34) 
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4.3.2. 5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


=  X 


l 


lp.(0)l2 

(X  V 

m 

(K  -  n] 

^2A-2 

b,(Q)l2 

IM 

A,  -  t*  J 

M°)l2 

(K) 

(K  -  v) 

2*-2 

Lpx(0)P 

IM 

A,  -  1*  J 

-A,  as  k  -  » 


(35) 


if  {'*rn}!Jl  satisfies  (33).  However,  this  estimate  is  n^t  monotonically  convergent  in  general,  since 
the  numerator  and  denominator  of  (35)  converge  at  slightly  different  rates. 

It  can  also  be  shown  that  the  power  ratio  F(w)  converges  to  Aj  as  k  -  <*,  since 


Fbaik)]  = 


eH(k)Rz(k) 

ll*(*)l2 


lg(0)A(A  -  MiP2*  1(0) 

&H(0)  (a  - 


=  Xj  x 


(i 

m 

‘ 

M\2* 

2  bi(0)!2 

l» 

1  > 

u  - 

2  lPx(0)  |2 

7k 

=  Xj  x 


s  X, 


it  lp.co>l2r 


1  - 


*  Lp,(0)  (' 


1  - 


'A.  -  !»?» 
.Ai  - 


,  .  glf.WlVA.  -  nf* 

2  1p,(0)IjI  a,  - 


with  equality  achieved  as  k  -  «  if  {Am)Y  satisfies  (33). 
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4.3. 2.5  ADAPTATION  CONVERGENCE  AND  STABILITY  ANALYSIS  (continued) 


As  shown  i-  ne  next  section,  the  convergence  factor  given  in  (34)  can  be  strongly  affected  by 
the  choice  of  p.  In  particular,  if  the  {Am}^  are  given  by 


i  .  1 1  + 

m  \  1 


m  =  ll 

m  >  1J 


where  Ymax  is  the  maximum  attainable  SINR  of  the  array,  then  setting  p  =  1  can  greatly 
reduce  the  ratio  given  in  (34). 


4. 3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS 

It  is  beyond  the  scope  of  this  report  to  determine  the  performance  of  the  Dominant  Mode 
algorithm  under  finite  time  conditions.  However,  useful  results  can  be  obtained  by  replacing 
sample  autocorrelation  matrices  in  the  true  algorithm  with  idealized  (infinite  time  average) 
statistics.  The  discrepancy  between  the  actual  and  idealized  quantities  should  be  small  if  the 
averaging  times  are  long  in  the  actual  algorithm. 

The  analysis  conducted  here  is  limited  to  determination  of  the  dominant  mode  of  the 
generalized  eigenequation  Aw  =  R0'1R1w  where  R0“lis  the  autocorrelation  of  the  received  data 
measured  over  the  off  interval  and  Ri  is  the  data  autocorrelation  measured  over  the  on 
interval.  This  analysis  is  conducted  for  the  following  scenarios: 

o  The  nominal  scenario,  where  both  Rq'1  and  Rj  are  properly  computed. 

o  Timing  error  scenarios,  where  the  on  and  off  time  data  collection  intervals  are  out  of 
synch  with  actual  SOI  timing. 

o  Collision  scenarios,  where  the  undesired  signals,  signals  not  of  interest  (SNOIs),  turn 
on  or  off  during  the  data  collection  intervals. 

The  analysis  which  follows  will  show  the  dominant  mode  eigenvalues  to  converge  to  the 
maximum  attainable  SINR  beamformer  in  all  the  nominal  and  timing  error  scenarios,  and  in 
many  of  the  collision  scenarios.  Consequently,  the  algorithm  is  expected  to  be  very  robust  in 
actual  practice. 

Nominal  Scenario  - 

The  instantaneous  SOI  power  is  illustrated  below  for  the  nominal  communication  scenario. 
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4. 3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 

The  off  time  and  on  time  measurements  are  made  over  the  correct  intervals,  yielding  the 
relations  x0g(t)  =  iofl(t)  and  xon(t)  =  ion(t)  +  as(t),  where  iofl(t)  and  ion(t)  have  identical 
autocorrelation  data  Ru.  The  off  time  and  on  time  autocorrelation  matrices  are  therefore  given 
by  R0=  Ruand  Rj  =  Ra  +  aaHRM.  Substituting  these  relations  into  the  generalized  eigenequation 
yields  m-1  solutions  with  aHwk  =  0  and  Xk  =  1,  and  one  solution  with  Wj  «*  R^aR,,  =  w,,^  and 
Xj  =  1  +  aHRy'1aR><  =  1  +  Ymax-  The  maximum  SINR  beamformer  is  wmiLX  and  the  maximum 
attainable  SINR  of  the  array  is  Ymax- 

From  the  preceding  section  (4. 3.2. 5),  it  can  be  seen  that  the  Dominant  Mode  algorithm  should 
converge  quickly  for  u  =  0  and  Ymax  large,  or  for  u  =  1  and  virtually  any  value  of  ymax. 
Furthermore,  this  mode  shall  correspond  closely  to  the  maximum  SINR  beamformer.  The 
dominant  mode  can  also  be  used  as  a  quality  statistic,  to  determine  whether  or  not  the 
algorithm  has  successfully  captured  the  SOI. 


;nanos  - 


The  instantaneous  SOI  power  is  illustrated  below  for  the  timing  error  scenarios. 

SCENARIO  #  1  SCENARIO  #  2 


These  figures  show  the  off  and  on  time  data  collection  intervals  to  be  out  of  synch  with  the 
true  SOI  on  time.  If  the  on  time  data  collection  interval  is  enabled  late,  some  of  the  on  time 
data  is  collected  during  the  off  interval;  if  its  enabled  early,  the  SOI  is  absent  over  part  of  the 
on  interval. 


In  the  first  timing  error  scenario,  the  off  and  on  autocorrelation  matrices  arc  approximated  by 
R0  =  Rji  +  eaaHR„  and  Rt  =  Ru  +  aaHRM.  The  eigenequation  has  m  -  1  solutions  with  aHwk  =  0 
and  Xk  =  1,  and  one  solution  with  Wj  «  and  Xj  -  (1  +  Ymax)/(1  +  eYmax)-  F°r  very  large 
SINRs,  Xj  approaches  1/e. 

The  second  scenario  gives  autocorrelation  matrices  R0  =  Ru  and  Rj  =  Ru  +  (l-e)aaHRfr  The 
eigenequation  also  has  m-1  solutions  with  a®wk  =  0  and  Xk  =  1,  and  one  solution  with  Wj  « 
and  Xj  =  1  +  (1  -  e)ymiuc- 

As  demonstrated  above,  in  both  cases  the  dominant  mode  of  Xw  =  Rq^RjW  corresponds  to  the 
maximum  SINR  beamformer.  These  calculations  can  be  verified  by  substitution  of  the 
applicable  autocorrelation  matrices  into  the  general  eigenequation. 


These  results  have  important  ramifications  for  the  Dominant  Mode  algorithm.  The  perf ormance 
of  the  algorithm  is  nearly  independent  of  the  synchronization  error  in  the  receiver.  The 
algorithm  should  be  able  to  capture  the  SOI  even  with  very  large  errors  in  timing.  A  simple 
post-capture  synchronization  circuit  can  then  be  designed  to  remove  this  timing  error  and 
optimize  performance  after  the  removal  of  the  other  stationary  interference. 
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4.3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 

Timing  Error  Scenarios  tcontlnuedl  - 

Timing  error  can  affect  the  convergence  speed  of  the  algorithm.  The  spread  between  the 
dominant  eigenvalue  and  the  lesser  eigenvalues  is  reduced  when  timing  errors  exist.  The  effect 
is  especially  pronounced  in  the  first  case:  is  limited  by  1/e  as  the  SINR  approaches  infinity. 
The  Dominant  Mode  algorithm  should  therefore  converge  much  more  slowly  in  this  case. 
However,  it  is  expected  that  the  convergence  of  this  algorithm  can  be  accelerated  by  setting  u 
=  1  in  the  algorithm  since  this  can  yield  an  arbitrarily  low  value  of  rconv  as  the  averaging  time 
grows  to  infinity. 

In  a  ground  based  Army  tactical  radio  application  the  timing  errors  analyzed  above  would 
occur  only  before  the  host  radio  has  been  acquired  by  the  network.  Once  in  the  network  the 
host  radio  would  provide  the  necessary  timing  information  to  the  AAP  to  allow  it  to  collect  off 
time  and  on  time  data  as  required  by  its  adaptation  algorithm.  When  the  host  radio  is  not  in  the 
network  it  monitors  the  activity  on  a  guard  band  channel  specifically  designated  for  acquiring 
synchronization  with  the  network.  The  AAP  could  easily  alter  its  adaptation  algorithm  to 
improve  the  host  radio’s  likelihood  of  acquiring  synchronization  if  the  AAP  is  inf ormed  by  the 
host  that  it  is  not  synchronized.  The  AAP  could  even  determine  this  condition  on  its  own  if  it 
monitored  the  consecutive  revisits  to  the  same  guard  band  frequency.  Once  the  AAP  knows  the 
host  radio  is  not  in  the  network,  it  could  increase  its  adaptation  time:  this  would  provide 
rejection  of  stationary  jammers  and  have  less  affect  on  the  SOIs.  The  jammers  would 
continuously  be  received  from  the  same  directions,  so  a  long  adaptation  time  would  suffice. 
The  SOIs,  however,  would  be  coming  from  different  directions,  or  would  occur  infrequently, 
so  they  would  not  be  rejected  as  quickly.  Jammer  rejection  would  provide  less  interference 
output  power;  SOIs  from  various  directions  would  decrease  the  likelihood  that  they  all  fall  in 
the  directions  of  inadvertant  nulls.  Therefore,  the  AAP  would  be  providing  improvement  of 
the  SINR  available  to  the  host  radio  relative  to  operation  without  the  AAP  most  of  the  time. 

An  extension  of  the  timing  error  scenario  is  the  condition  that  an  expected  SOI  is  not  present. 
The  AAP  expects  a  SOI  during  an  on  time  data  collect,  but  in  network  communications  this  is 
not  guaranteed  to  occur.  The  reasons  are  many  and  varied:  All  of  the  transmit  time  periods 
available  may  not  be  used  by  the  network;  the  radio  transmitting  during  a  specific  timeslot  may 
be  so  far  away  from  the  radio  which  is  host  to  the  AAP  that  the  received  signal  level  is  too  low 
to  be  detected;  or  the  AAP  may  be  shielded  from  some  SOIs.  In  these  cases,  the  AAP  should 
maintain  its  previously  developed  weights  rather  than  adapt  to  an  erroneous  environment.  The 
Dominant  Mode  Block  Update  algorithm  handles  this  circumstance  very  well:  the  weight  update 
is  based  on  calculations  of  R0_1Ri.  Since  there  is  no  SOI  present  during  Rx  data  collection,  and 
the  jammers  are  assumed  stationary,  the  Rt  data  should  exhibit  the  same  statistics  as  the  R0 
data  in  this  case.  Calculations  of  R^Rj  should  repeatedly  approach  the  identity  matrix  I,  and 
the  new  weight  estimate  would  equal  the  old.  Simulations  verify  this  expected  result.  In 
addition,  the  algorithm  could  easily  determine  that  no  SOI  is  present  and  discontinue 
adaptation  before  degrading  previously  formed  weights. 
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4.3.2. 6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 
Collision  Scenarios  - 
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A  collision  is  said  to  occur  when  a  signal  not  of  interest  (SNOI)  changes  state  during  a  data 
collection  interval.  The  four  most  likely  collision  scenarios  are:  1)  A  SNOI  turns  off  during  an 
off  time  data  collect,  2)  a  SNOI  turns  on  during  an  off  time  data  collect,  3)  a  SNOI  turns  off 
during  an  on  time  data  collect,  and  4)  a  SNOI  turns  on  during  an  on  time  data  collect. 

SCENARIO  #  1  SCENARIO  #2 
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SCENARIO  #3 


SCENARIO  #4 
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The  general  equations  used  in  the  analysis  of  these  collision  scenarios  are: 


t 


i+(t)  •=  u(t)  +  «il(t) 

A  «  [  a,  ] 

R0  *=  Ruu  +  ASjjA®  or  R0  »  R++  +  AS„j+Ah 

R|  “  Run  +  AS^A®  Ri  *  R++  + 

The  S  matrices  depend  on  which  of  the  collision  scenarios  are  being  analyzed.  The  other  terms 
are  defined  by: 

Ruu  «  autocorrelation  matrix  of  steady  state  interferers  only 

R  ,  |  *>  autocorrelation  matrix  of  steady  state  interferers  and  transient  interferer 

a j  =  aperature  of  transient  interferer 

a,  =  aperature  of  the  desired  signal 

Rii  ■=  covariance  of  transient  interferer 

R„  *=  covariance  of  the  desired  signal 


Substituting  the  first  definitions  for  R0and  Rj  into  the  general  eigenequation  gives 

MRUU  +  AS^V  -  [R„u  +  ASWAV 
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4.3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 

Collision  Scenarios  (continued)  - 

The  solutions  to  this  eigenequation  includes  m  -  2  solutions  with  A®wk  =  0  and  Xk  ■  1,  and  two 
solutions  with  wm  *  Ruu“1Apm  where  pm  satisfies  the  eigenequation 

A[I  +  S^A^^AHp  =  [I  +  S^A^^AjJp 

The  other  resulting  general  eigenequation  is 

X[R++  +  ASoff+AH]w  =  [R++  +  ASon+AH]w 

In  this  case,  the  solutions  of  interest  are  =  R++_1Apm+  where  pm+  satisfies  the  eigenequation 

MI  +  Soff+(AHR++'iA)]p+  =  [I  +  Son+(AHR++*1A)]p+. 

First  case  -  In  this  case  a  SNOI  turns  off  during  an  off  time  data  collect.  The  resulting 
autocorrelation  matrices  are  R0  =  Ruu  +  and  Ri  =  Rnu  +  Substituting  these 

into  the  generalized  eigenequation  gives 

X(RUU  +  eajaiHRH)w  =  (Ruu  +  a,a.HRM)w. 

This  eigenequation  has  m-2  solutions  with  Xk  =  i,  a^*,  and  a,11^  =  0;  the  other  two  solutions 
are  the  solutions  of  interest: 

1)  The  maximum  (dominant)  eigenvalue  solution,  which  gives  w,  <*  R,,,,"^,  and  Xj  *  1  +  Yimwt- 

2)  The  minimum  eigenvalue  solution,  which  gives  wm  *  Ruu^aj  and  Xm  •  1/(1  +  eYim8X). 

These  results  will  hold  true  for  any  value  of  e  between  0  and  1,  even  in  the  extreme  case  when 
the  interferer  turns  off  at  exactly  the  same  time  that  the  SOI  turns  on.  The  algorithm  should 
behave  robustly  in  the  presence  of  any  collision  which  actually  occurs  in  the  off  interval  (i.e., 
where  the  SNOI  energy  does  not  actually  interfer  with'the  SOI  waveform).  This  suggests  that 
the  algorithm  should  be  compatible  with  dense  FHMA  and  TDMA  communication  formats,  as 
long  as  the  signals  are  not  time  coincident.  This  underscores  the  applicability  of  the  algorithm 
to  multi-node  communication  networks. 

However,  the  convergence  speed  of  the  Dominant  Mode  algorithm  can  be  affected  by  these 
collisions,  and  the  choice  of  u  becomes  more  significant  in  this  environment.  In  particular,  the 
choice  of  u  *  1  can  cause  the  antenna  array  to  beam  up  on  the  interferer  if  the  maximum 
attainable  SINR  of  the  SOI  is  very  low.  If  the  magnitude  of  Xm  -  u  is  greater  than  the 
magnitude  of  Xx  -  u,  or  equivalently,  if  (Xx  +  Xm)/2  is  greater  than  u,  then  the  algorithm  will 
capture  the  minimum  mode  and  reject  the  SOI.  Note  the  algorithm  will  not  capture  the 
interferer  since  it  is  off  before  the  on  time  data  collection  interval  when  the  beamformer 
weights  are  actually  computed  and  used. 

The  algorithm  will  capture  the  minimum  mode  only  if  y(  satisfies  the  inequality 


Y,  <  eYj/(l  +  «Y;)  -  2(1  -  u). 
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4. 3. 2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 
Collision  Scenarios  ^continued)  - 
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This  implies  y,  <  2u  -  1  and  Yj  >  f  1  /(2u  -  1  -  ys)  -  1  ]/e.  This  condition  can  be  precluded  by 
setting  u  s  i/2,  and  is  not  a  threat  if  Y*max  is  greater  than  unity.  However,  if  u  =  1  is  chosen, 
then  the  convergence  time  may  be  dominated  by  1  -  Xm  rather  than  X2  -  1. 

Second  case  -  In  this  case  a  SNOI  turns  on  during  an  off  time  data  collect.  The  resulting 
autocorrelation  matrices  are  R0  =  R++  -  (1  -  and  R2  =  R++  +  a,a,HR>g.  Substituting 

these  into  the  generalized  eigenequation  gives 

X(R++  *  (1  -  e^^R^w  -  (R++  +  atatHR„)w. 

This  eigenequation  has  two  solutions  of  interest: 

1)  The  SOI  capture  solution,  which  gives  w,  «  R++_1a,  and  X,  =  1  +  y8+- 

2)  The  transient  interferer  solution,  which  gives  Wj  «  R+^^iand  X;  =  1/(1  -  (l-e)Yi+)- 

In  these  equations:  1)  Yi+  is  the  maximum  attainable  SINR  over  the  on  time  and  considers  the 
SOI  to  be  the  signal  rather  than  an  interferer.  2)  Yj+  is  a  function  of  Yj,  the  maximum 
attainable  SOI-free  Interference  to  Noise  Ratio  (INR),  which  ignores  the  SOI  and  treats  the 
transient  interferer  as  the  desired  signal.  The  term  Yi+  is  derived  as  follows: 

Yi+  =  a^R^ajRii 

“=  +  aiaiHRiil'laiRii 

-  aiHRuu'laiRii/(i  +  aiHR»u'laiRii) 

=  Yj/O  +  Yi). 

Note  that  while  Yj  has  a  range  from  0  to  infinity,  the  range  of  Yj+  is  limited  between  0  and  1. 

For  this  scenario,  it  is  possible  for  the  dominant  mode  of  the  eigenequation  to  capture  the 
interferer:  both  eigenvalues  can  have  a  magnitude  greater  than  unity,  and  the  eigenvalue 
associated  with  the  capture  of  the  SOI  can  be  less  than  the  eigenvalue  associated  with  the 
capture  of  the  interferer.  By  comparing  the  associated  eigenvalues,  the  magnitude  of  Xt  is 
greater  the  magnitude  of  X§  if  yi+  satisfies 

(1  *  I/O  +  Yi+))/(l  -  e)  <  Yi+  <  (1  +  l/d  +  Yi+))/(l  -  *) 

This  yields  two  conditions: 

Yi+  <  (1  -  e)Yi+/(l  *  (1  -  «)Yi+)  where  Yi+  <  1/(1  *  «)  and 

Y.+  <  (2  -  (1  -  e)Yi+)/(l  -  e)Yi+  -  1)  where  1/(1  -  e)  <  Yi+  <  2/(1  -  e). 
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4.3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 
Collision  Scenarios  (continued)  - 
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The  critical  value  of  y»+  is  plotted  as  a  function  of  y-l+  in  figure  l.a.  The  critical  value  spans 
the  entire  range  of  yi+,  approaching  infinity  for  Yj+  =  1/(1  -  e)  However,  only  a  portion  of  this 
region  is  pertinent,  since  Yi+  is  bounded  by  0  and  1.  In  particular,  Yj+  is  less  than  1/(1  -  e)  for 
every  value  of  e  (0  <  e  <  1).  Using  these  relations,  a  bound  on  yi+  that  is  directly  dependent 
on  Yi  (the  true,  maximum  attainable  INR)  can  be  found: 

Y»+  «=  (e'1  -  l)cYi/(l  +  «Yi). 

This  threshold  is  plotted  as  a  function  of  Y;  in  figure  l.b. 

As  is  evident  in  figure  1,  the  Dominant  Mode  algorithm  is  guaranteed  to  capture  the  SOI  if  the 
maximum  attainable  SOI  SINR  is  greater  than  e"1  -  1.  Assuming  that  t  is  uniformly  distributed 
between  0  and  1,  then  the  median  value  of  critical  y1+  is  given  by  Yi/(2  +  Yi)  which  is  less  than 
1  for  all  Y;- 

The  practical  significance  of  this  result  is  that  this  collision  scenario  will  result  in  capture  of 
the  SOI  in  over  50  %  of  all  collisions  if  the  maximum  attainable  SINR  of  the  SOI  is  greater 
than  0  dB.  The  SOI  will  be  captured  in  over  75  %  of  these  collisions  if  its  maximum  attainable 
SINR  is  greater  than  5  dB. 

This  collision  scenario  will  directly  affect  the  convergence  time  of  the  Dominant  Mode 
algorithm,  by  decreasing  the  spread  between  the  two  largest  eigenvalues.  This  convergence  time 
will  not  be  significantly  affected  by  choice  of  u  in  the  algorithm. 

Third  and  fourth  cases  -  The  remaining  two  collision  scenarios  yield  a  mixture  of  the  SOI  and 
interference  waveforms.  In  scenario  #  3,  S<m+  and  Sofr+  are  given  by 


K  o 

°  tRu 

which  yields  the  eigenequation 


ZOFF+ 


0  0 
0  0 


and  S. 


hon* 


ePu.Ra  tPu*Ra 
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INTERFERENCE  _ /  ~ ;  — 

CAPTURE  REGION  fc  * 

FIGURE  1JL  THRESHOLD  VALUE  OF  HWVS.*U . 


CAPTURE  REGION 

FIGURE  I.B.  THRESHOLD  VALUE  OF  VS. 


FIGURE  1.  COLLISION  THRESHOLD  VALUES 
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4.3.2.6  RECEIVE  ADAPTATION  PERFORMANCE  ANALYSIS  (continued) 
In  scenario  #  4,  the  corresponding  S  matrices  are  given  by 


*OFF 


0  0 
0  0 


HoN 


K  o 

0  (1  -  t)Ru 


ZOFF* 


0  0 
0  -R, 


yielding  the  eigenequations 


;  and  £ 


ON* 


*.  0 

0  zRu 


P„R„  P«R* 

(1  -  t)pisRu  (1  -  e)paKaj 


1  0 

P«*K 

Psi*R . 

o  1  -  e/y 

P*  - 

*Pl»*RU 

tPu*Ru 

Both  of  these  eigenequations  apply,  and  are  appropriate  at  different  times  in  the  collection 
interval. 

It  is  beyond  the  scope  of  this  report  to  analyze  these  scenarios  in  any  detail.  However,  it  is 
expected  that  the  dominant  mode  weight  vector  will  be  weighted  towards  the  SOI  if  e  is  small 
in  scenario  #  3  and  if  e  is  large  in  scenario  #  4.  It  is  also  expected  that  this  weighting  shall  be 
dependent  on  the  maximum  attainable  interference  SINR.  However,  this  affect  should  be 
limited,  especially  for  scenario  #  3,  since  Yj+  is  limited  between  0  and  1  but  y»  is  unbounded. 

In  any  case,  some  measures  can  be  taken  to  mitigate  against  these  collisions.  The  last  two 
collision  scenarios  can  all  be  detected  by  monitoring  the  power  level  of  the  beamformer  output 
signal.  In  LPI  communications,  if  the  SOI  has  been  captured  on  previous  hops,  the  previous 
aperature  estimates  can  also  be  used  to  detect  capture  conditions:  if  the  current  aperature 
estimate  is  vastly  different  than  the  history  of  aperature  estimates,  it  can  be  assumed  that  the 
algorithm  is  acquiring  the  wrong  signal,  and  that  timeslot’s  processing  can  be  discarded. 

4.3.2.7  DIRECTIVE  /  RETRODIRECTIVE  COMMUNICATION 

The  receive  mode  adaptation  simultaneously  steers  nulls  in  the  directions  of  jammers  and  a 
beam  in  the  direction  of  the  desired  signal.  This  receive  antenna  pattern,  which  is  customized 
to  the  RF  signal  environment,  provides  tremendous  improvement  to  the  SINR  seen  by  the 
receiver  compared  to  omni-directional  antenna  performance.  This  same  information  can  be 
used  to  improve  the  transmit  antenna  pattern  as  well.  Two  methods  of  using  the  receive 
statistics  and  weights  to  direct  the  transmit  energy  are  considered:  retrodirective  schemes  and 
directive  schemes. 
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4. 3.2. 7  DIRECTIVE  /  RETRODIRECTIVE  COMMUNICATION  (continued) 

Retrodirective  transmission  schemes  apply  the  complex  conjugates  of  the  receive  weights  to  the 
transmit  weighting  elements.  This  is  described  by  wtx  =  w^*  where  wtx  is  the  transmit  weight 
vector  and  wra  is  the  receive  weight  vector  formed  by  adaptation.  If  the  receive  and  transmit 
paths  are  precisely  calibrated,  this  will  result  in  a  transmit  antenna  pattern  which  is  identical 
to  the  receive  pattern.  The  advantage  of  this  approach  is  that  very  little  power  is  transmitted 
in  the  direction  of  the  jammers  since  nulls  have  been  formed  in  those  directions. 
Simultaneously,  there  is  high  gain  (a  beam)  in  the  direction  of  the  last  signal  received  at  the 
current  channel  frequency. 

Directive  transmission  schemes  ignore  the  jammer  locations,  and  attempt  to  maximize  the 
power  transmitted  in  the  direction  of  the  desired  signal.  The  nulls  in  the  directions  of  jammers 
are  not  maintained,  so  the  jammers  will  see  more  transmit  power  in  this  case  than  for 
retrodirective  cases,  but  more  power  can  be  provided  in  the  desired  direction  and  less  in  all 
other  directions.  To  achieve  directive  communications,  the  transmit  weights  are  set  to  the 
complex  conjugate  of  the  receive  aperature  estimates  determined  during  receive  adaptation. 
The  pertinent  equation  is  wu  = 

The  retrodirective  transmission  weights  are  very  precise  since  they  are  based  on  the  high  level 
jammers  present  during  reception.  Very  deep,  accurately  directed  nulls  are  provided,  and  the 
desired  signal  direction  has  relatively  high  gain.  The  directive  transmission  weights  are  not  as 
precise,  at  least  not  with  aperature  estimates  based  on  a  limited  number  of  receive  data  samples. 
In  this  case,  the  jammer  power  corrupts  data  being  used  to  estimate  the  desired  signal’s 
direction  of  arrival.  Inaccuracies  in  this  estimate  directly  impacts  the  quality  of  the  directive 
transmit  antenna  pattern.  This  performance  can  be  greatly  improved,  however,  by  averaging 
a  large  number  of  aperature  estimates  taken  over  a  long  period  of  time.  There  are  a  variety  of 
ways  this  averaging  can  be  accomplished;  some  schemes  are  outlined  below: 

o  The  receive  aperature  estimates  for  each  frequency  hop  can  be  stored.  This  historical 
data  can  be  averaged  to  form  the  transmit  weights  for  a  given  frequency. 

o  The  history  of  receive  aperature  estimates  for  each  frequency  can  be  averaged,  then 
used  to  interpolate  the  aperature  estimates  for  frequencies  not  previously  visited  or  rarely 
visited. 


o  The  receive  aperature  estimates  for  each  frequency  hop  can  be  included  in  the 
averaging  of  the  aperature  estimates  for  some  of  the  adjacent  channels  as  well  as  its  own 
frequency  channel.  This  will  reduce  the  number  of  visits  to  each  f requency  required  to  obtain 
an  accurate  transmit  beam. 

o  The  output  SINR  provided  by  receive  adaptation  can  be  calculated  by  the  AAP.  This 
information  can  be  used  to  weight  the  aperature  estimate  calculated  during  that  particular 
timeslot.  For  example,  if  the  output  SINR  is  very  high,  the  aperature  estimate  is  more  accurate 
than  a  low  SINR  output  would  provide,  so  the  weight  of  this  estimate  should  be  relatively  large. 
These  weighted  aperature  estimates  would  then  be  averaged  to  determine  the  transmit  weights. 

Certainly  many  other  techniques  could  be  suggested,  and  the  techniques  listed  above  could  be 
combined  to  provide  more  options. 
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4.3.2.7  DIRECTIVE  /  RETRODIRECTIVE  COMMUNICATION  (continued) 

Directive  transmission  techniques  can  be  combined  with  receive  mode  adaptive  processing  and 
careful  transmit  power  level  control  to  provide  a  directive  communication  link  which  is  very 
difficult  to  jam  or  intercept.  This  anti-jam  (A/J)  and  low  probability  of  intercept  (LPI)  link 
would  operate  as  described  below: 

1)  Radio  1  transmits  omni-directionally,  intending  to  communicate  with  Radio  2. 

2)  Radio  2  enables  the  receive  adaptation  of  its  adaptive  array  processor.  This  rejects 
any  interference  present  and  assures  high  gain  in  the  direction  of  Radio  1.  Even  if  the  radios 
are  not  in  synch,  the  algorithm  can  be  adjusted  such  that  the  signal  is  acquired. 

3)  After  Radio  2  has  acquired  the  signal  from  Radio  1,  and  has  formed  nulls  in  the 
directions  of  jammers,  uses  this  information  to  determine  its  transmit  weights.  These  weights 
are  designed  to  provide  high  gain  in  the  direction  of  Radio  1  and  little  gain  elsewhere. 

4)  Radio  1  will  employ  its  AAP  to  reject  interferers  and  beam  up  on  the  signal  from 
Radio  2.  It  will  then  be  capable  of  directive  transmissions  as  well. 

5)  At  this  point  the  link  between  Radio  1  and  Radio  2  is  established.  The  transmit 
antenna  patterns  are  optimized  to  transmit  most  of  their  transmit  energies  in  the  direction  of 
each  other;  their  receive  antenna  patterns  are  optimized  to  provide  the  maximum  attainable 
SINR  to  their  demodulators.  Each  radio  could  calculate  their  SINR,  and  knowing  the  SINR 
required  for  reliable  communications  they  could  determine  how  much  margin  was  available. 
This  information  could  be  transmitted  back  to  the  other  radio,  which  would  then  know  how 
much  less  power  it  could  transmit  and  still  maintain  a  reliable  link. 

With  transmissions  at  minimum  power  and  antenna  gains  beamed  up  for  the  link,  and  jammers 
nulled  during  reception  while  simultanously  providing  gain  for  the  link,  a  very  robust 
communication  link  which  is  difficult  for  the  enemy  to  detect  is  provided. 

4.3.2.8  DIRECTIVE  COMMUNICATION  PERFORMANCE  ANALYSIS 

Evaluation  of  the  directive  transmission  performance  depends  on  the  application;  two 
applications  are  considered  in  this  study: 

LEI- 

The  power  radiated  from  the  array  is  calculated;  it  is  highly  dependent  on  direction.  A  locus 
of  points  around  the  transmitter  are  calculated  which  indicate  the  geographical  area  within 
which  the  radiated  power  is  sufficient  for  detection  by  a  reconnaissance  receiver.  The 
boundary  of  this  region  is  identical  to  the  transmit  antenna  pattern,  only  depicted  on  a  linear 
scale  rather  than  logarithmic,  and  includes  a  linear  scaling  factor  which  reflects  the  sensitivity 
of  the  recon  receiver.  This  locus  of  points  establishes  a  threshold  of  transmit  power;  receivers 
within  this  boundary  will  detect  the  transmissions,  while  receivers  outside  this  boundary  will 
not.  The  area  inside  this  boundary  is  referred  to  here  as  the  intercept  footprint. 

The  directive  transmissions  will  result  in  some  intercept  footprint  which  greatly  reduces  the 
intercept  footprint  relative  to  isotropic  transmissions,  but  which  may  not  be  optimal.  The 
optimal  pattern  would  minimize  the  intercept  footprint  area  while  maintaining  consistent  gain 
in  the  direction  of  the  desired  link.  Thus  the  figure  of  merit  for  LPI  communications  is  the 
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4.3. 2.8  DIRECTIVE  COMMUNICATION  PERFORMANCE  ANALYSIS 
LPI  (continued)  - 

areas  of  the  intercept  footprint.  The  optimum  radiated  pattern  is  derived  below,  and  is  used 
for  comparison  to  the  directive  patterns  generated  via  receive  adaptation  in  the  Simulation 
section  of  this  report  (4.3.4).  Isotropic  data  are  also  used  in  the  comparisons. 

Derivations  -  The  terms  used  in  the  derivation  of  the  optimal  transmit  antenna  weights  and 
pattern,  and  in  the  computation  of  the  intercept  footprint  areas,  are  defined  below: 

dICV  =  distance  from  xmitter  to  revr 
exmt  =  direction  from  xmitter  to  revr 

Brecon  =  distance  from  xmitter  to  recon  system 
® recon  ■  direction  from  xmitter  to  recon  system 

Ylpi  =  required  revr  SINR 
Ydet  =  required  recon  SINR 

NFKy  m  revr  noise  figure 
NFn con  ■  recon  noise  figure 

&xmt(0)  *  emitter  array  manifold 
£xmt  “  xmitter  array  weightings 

«  number  of  xmitter  sensors 
Mkv  *  number  of  revr  sensors 

When  both  the  transmitter  and  receiver  are  isotropic,  the  ratio  of  distances  that  the 
reconnaissance  and  desired  receivers  are  located  from  the  transmitter,  both  capable  of 
detecting  the  transmission,  is  given  by 


^ recon  _ 

fYLP/) 

A.  \ 

When  the  transmitter  is  isotropic  but  the  desired  receiver  employs  an  AAP  and  has  established 
the  optimum  antenna  pattern,  this  ratio  becomes 


^ recon 

l 

f  V  LPI 

dec 

=  N 

\(  NF. 


rev 


\  NF  recon  j 


v  e 


If  both  the  transmitter  and  receiver  employ  antenna  arrays,  and  directive  transmissions  are 
provided,  the  distance  ratio  is 
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LPI  (continued^  - 
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A, recon 
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The  area  of  the  intercept  footprint  corresponding  to  the  isotropic  transmission  and  adaptive 
reception  is  given  by 


where  M,.ec 


A 


iso 


for  non-adaptive  reception. 


1 

Y  LPI 

NF  \ 

rev 

NF 

V  recon  / 

The  area  in  the  directive  transmission  case  is 


=  A 


8xmt 


8xmft  ncot) 


iso 


Is 


where 


;(0)  0^(6)  dQ 


-7t 


The  optimum  transmit  weights  for  LPI  directive  communications  would  be  those  weights  which 
minimize  the  intercept  footprint  area.  The  solution  to  the  minimization  of  Atg^)  is 


8 opt  ** 

The  resulting  minimum  intercept  footprint  area  is  then 


A/J  - 


The  concept  for  A/J  applications  is  very  similar  to  LPI  except  the  reconnaissance  receiver  is 
replaced  by  a  jammer,  and  the  location  of  the  jammer  relative  to  the  receiver  is  of  importance. 
The  transmit  power  is  assumed  to  be  fixed  -  the  transmitter  will  transmit  at  its  maximum  level 
during  all  transmissions,  disregarding  the  likelihood  of  being  detected  by  a  reconnaissance 
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4. 3. 2. 8  DIRECTIVE  COMMUNICATION  PERFORMANCE  ANALYSIS 
A/J  (continued)  - 

receiver.  The  directive  communication  system  must  use  the  adaptive  antenna  arrays  at  the 
receiver  and  transmitter  to  provide  the  maximum  power  into  the  receiver  given  the  fixed 
transmit  powe.  limitation.  Thus,  in  the  A/J  case,  the  figure  of  merit  is  the  SOI  output  power 
of  the  receiver’s  AAP  (given  the  fixed  transmit  power  constraint). 

The  calculations  in  the  A/J  case  are  very  similar  to  the  LPI  case.  As  was  done  for  the  LPI 
analysis,  the  optimum  radiated  pattern  is  derived  below,  and  is  used  for  comparison  to  the 
directive  patterns  generated  via  receive  adaptation  in  the  Simulation  section  of  this  report 
(4.3.4).  Isotropic  data  are  also  used  in  the  comparisons. 

Derivations  -  The  terms  used  in  the  derivation  of  the  optimal  transmit  antenna  weights  and 
pattern,  and  in  the  computation  of  the  receiver  input  power,  are  defined  below: 

rfrev  =  distance  between  xmitter  and  rcvr 
6rcv  =  DOA  of  xmitter  at  rcvr 

d jam  =  distance  between  jammer  and  rcr 
0jam  *=  DOA  of  jammer  at  rcvr 

Ya/j  *  required  rcvr  SINR 

ailo  =  isotropic  jammer-free  link  margin 

Pxmt  =  total  xmitter  ERP 
Pjam  *  total  jammer  ERP 

flrcv(0)  ■  rcvr  array  manifold 
£xmt  =  emitter  array  weightings 

Afxmt  *  number  of  xmitter  sensors 
A/rev  *=  number  of  rcvr  sensors 

For  A/J  applications,  the  distances  of  interest  are  the  jammer  and  the  desired  receiver  relative 
the  the  transmitter;  the  analysis  uses  the  ratio  of  these  two  distances.  When  both  the  transmitter 
and  receiver  are  isotropic,  this  distance  ratio  is  given  by 


fal)  (Pjm' 

^  “  aiL  l 


When  the  transmitter  is  isotropic  but  the  desired  receiver  employs  an  AAP  and  has  established 
the  optimum  antenna  pattern,  this  ratio  becomes 
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4.3. 2.8  DIRECTIVE  COMMUNICATION  PERFORMANCE  ANALYSIS 


d  \2  /P.  \  I 

~~  =  MrCyaMyaJj  -jr 

°rcv  rxmt 


_ _jSL _  _  1 

1  ' 


G,„  - 

rev 


In  the  fully  adaptive  case,  both  the  transmitter  and  the  receiver  employ  antenna  arrays,  and 
directive  transmissions  are  provided.  The  resulting  distance  ratio  is 


^!=)2  =  J-^]max|o, - 2s= - - 

“J  *v-j  1  - 


The  optimum  transmit  weights  for  A/J  directive  communications  would  be  those  weights  which 
maximize  the  receiver  input  power.  These  weights  are  given  by  the  solution  to 


“a*  I sL  axm$ ntt)  I  mch  that  te«J2  = 


results  in  the  optimal  A/J  transmit  weights 


xmt^  pl/2 
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4.3.3  ARCHITECTURE  DEVELOPMENT 
4.3.3. 1  INTRODUCTION 

To  implement  the  algorithm  developed  in  Phase  I,  the  architecture  must  include  specific 
features;  different  algorithms  will  require  other  features.  The  intent  of  the  ACM  Systems 
architecture  development  was  to  provide  all  the  functions  necessary  to  implement  a  large 
variety  of  algorithms,  and  simultaneously  allow  compatibility  with  an  existing  Army  tactical 
radio.  One  of  the  long  term  goals  of  this  study  is  to  lead  to  the  hardware  development  of  an 
A  AP  unit  which  could  be  produced  in  high  volume  and  deployed  with  Army  radios  already  in 
the  field.  The  corresponding  constraints  on  the  AAP  unit  requires  the  architecture  to  be  small, 
have  low  complexity,  low  cost,  and  accomodate  a  range  of  applications. 

To  allow  adaptation  of  the  weights,  a  minimum  number  of  RF  statistical  measurements  are 
required.  Depending  on  the  algorithm  being  used,  the  statistics  may  be  accumulated  in  vector 
form  (such  as  the  input  to  output  cross-correlation  vector)  or  in  matrix  form  (such  as  the  input 
covariance  matrix).  If  the  algorithm  is  minimizing  output  power  when  it  is  known  that  only 
interference  plus  noise  is  present  (no  signal),  then  "off  time"  data  only  is  required.  However, 
if  the  algorithm  is  maximizing  the  Signal  to  Interference  plus  Noise  Ratio  (SINR),  then  data 
must  be  acquired  when  the  signal  is  present  "on  time”  and  when  not  present  "off  time".  As 
discussed  in  the  Algorithm  sections,  "on  time”  and  "off  time"  refers  to  the  expected  presence  of 
the  signal  of  interest.  This  section’s  notation  for  the  correlation  vectors  and  matrices 
corresponding  to  off  time  and  on  time  data  collects  are  defined  below: 


Off  time  cross-correlation  vector  *=.  R^  = 

1  E(xjy) 

1 

1  1 

1 

1 

1 

(measured  when  signal  is 

1 

1 

known  not  to  be  present) 

1  E(xmy) 

1 

On  time  cross-correlation  vector  =  R^j  = 

1  E(xiy) 

i 

1  1 

1 

1 

1 

(measured  when  signal  is 

1 

1 

expected  to  be  present) 

i  E(*my> 

1 

1 

l 

E(x1x1) 

E(x,xm) 

1 

Off  time  covariance  matrix  =  R0  *  | 

1 

1 

(measured  when  signal  is  ! 

known  not  to  be  present)  | 

E(*mxl) 

E(xmxm) 

I 

1 

EU^j) 

E(xxxm) 

1 

On  time  covariance  matrix  =  Rj  *  j 

I 

1 

(measured  when  signal  is  | 

expected  to  be  present)  | 

E(*mxl) 

E(xmxm) 

In  the  definitions  above,  E( )  denotes  the  expected  value  of  the  product  inside  the  brackets.  The 
baseline  architecture  design  is  capable  of  measuring  all  of  the  statistics  identified  above;  the 
baseline  algorithm  assumed  for  subsequent  discussion  is  the  Dominant  Mode  Block  Update. 
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4.3.3  ARCHITECTURE,  INTRODUCTION  (continued) 

The  minimum  statistical  data  required  f or  the  off  time  is  the  collection  of  covariance  data  and 
the  calculation  of  Rq'1.  During  the  on  time,  the  requirements  are  more  demanding:  A  block  of 
on  time  cross-correlation  data  is  measured.  This  is  combined  with  R0_1  and  the  last 
calculated  weight  estimate  w(m)  to  calculate  an  update  of  the  weight  estimate  w(m+l).  This 
procedure  is  iterated  with  sequencial  blocks  of  on  time  data,  resulting  weights  which  converge 
to  values  which  provide  the  maximum  attainable  output  SINR.  Several  iterations  are  made 
during  each  timeslot;  as  the  weights  are  updated,  the  receive  antenna  gain  pattern  changes, 
ultimately  steering  nulls  in  the  directions  of  jammers  and  a  beam  in  the  direction  of  the  desired 
signal. 

In  a  frequency  hopped  system  the  on  time  is  entirely  defined  by  the  tactical  radio  the  AAP  is 
operating  with;  there  is  no  flexibility  in  when  on  time  data  may  be  collected.  The  off  time  data, 
howe'er,  may  be  separated  from  the  on  time  either  spectrally  or  temporally  (by  frequency  or 
time).  This  is  illustrated  in  figure  2.  Spectrally  separated  off  time  data  is  collected  by  sampling 
the  RF  statistics  in  a  frequency  bin  adjacent  to  the  on  time  frequency,  and  is  done 
simultaneously  with  on  time  processing.  Temporally  separated  off  time  data  is  collected  at  the 
same  frequency  as  the  on  time  signal,  but  at  a  time  when  the  signal  is  not  on.  The  baseline 
selected  by  ACM  for  temporally  separated  data  collection  uses  the  immediately  preceeding 
timeslot,  and  is  therefore  called  anticipative. 


RANDOM  FREQUENCY  HOPPING 


777, 
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7//7, 

yy// 

'///, 
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time  slots 


y//y\  =  ON  TIME  —  a  1  s  +  I  +  N 

|  1  =  OFF  TIME  — -  &  Q  — *  |  +  N 

OPTIONS  FOR  COLLECTING  OFF  TIME  DATA: 

1)  SPECTRAL  SEPARATION 

OR 

2)  ANTICIPATIVE  TIME  SEPARATION 


Figure  2.  Data  Collection  Timing. 

The  spectral  and  temporal  means  of  off  time  data  collection  provide  two  candidate  top  level 
architectures;  these  are  both  presented  and  a  baseline  is  selected  and  justified.  This  leads  to 
examination  and  selection  of  each  feature  of  the  architecture  implementation.  Topics  discussed 
below  includes  the  key  electrical  functions  (complex  weighting  elements,  statistical 
measurements,  and  weight  calculation),  other  f unctions  specific  to  the  chosen  architecture  (RF 
assembly,  IF  chains,  frequency  synthesizer,  interface  and  control,  etc.),  the  necessity  of  some 
calibration  and  proposed  methods  to  achieve  it,  a  complete  baseline  receive  AAP  block  diagram, 
and  some  comments  regarding  the  transmit  architecture  required  for  directive  or  retrodirective 
communications. 
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4.3.3.2  DATA  COLLECTION  AND  TOP  LEVEL  ARCHITECTURES 
Introduction  - 

As  discussed  above,  two  methods  of  separating  off  from  on  time  data  are  considered  -  spectral 
and  temporal.  A  section  below  is  dedicated  to  each  method,  and  includes  block  diagram  of  the 
top  level  architecture  capable  of  implementing  them.  In  both  cases,  the  following  key  function 
implementations  are  assumed:  The  complex  weights  are  realized  at  IF,  the  statistical 
measurements  are  made  at  both  IF  and  baseband,  and  the  weight  calculations  are  done  digitally. 

Spectral  Separation  - 

A  block  diagram  showing  how  on  time  and  off  time  data  can  be  collected  simultaneously  using 
spectral  separation  is  given  in  figure  3.  The  basis  for  this  technique  are  two  statistical 
measurement  circuits  in  each  channel  operating  in  parallel.  ("Channel"  refers  to  the  •'ifferent 
antenna  elements  and  their  associated  circuits).  The  data  provided  to  the  two  statistical 
measurement  circuits  are  from  different  frequency  bands  by  virtue  of  the  filtering  provided 
by  FLTR  0  (off  time)  and  FLTR  1  (on  time). 

The  incoming  RF  signals  are  sensed  by  the  antenna  elements  and  processed  by  the  receiver 
front  ends  which  filter  out  signals  beyond  the  operating  range  of  the  host  radio  and  establish 
a  low  noise  figure  for  the  AAP  system.  A  Local  Oscillator  signal  (LO)  is  provided  to  the 
downconverters  which  hops  synchronously  with  the  host  radio  to  cause  the  signal  of  interest 
(SOI)  to  always  fall  within  the  fixed  IF  passband  provided  by  FLTR  1.  The  same 
downconverter  is  used  to  simultaneously  cause  an  adjacent  RF  band  to  f all  at  an  offset  IF  band 
which  is  passed  by  FLTR  0.  The  circuits  following  these  two  highly  selective  IF  filters  all  use 
wideband  devices  capable  of  processing  both  bandpass  signals  equally  well,  and  are  therefore 
identical.  The  IF  feedfoward  chains  are  required  to  process  the  IF  signals  so  they  are  within 
the  most  linear  region  of  the  statistical  measurement  circuits’  operating  range.  The  statistical 
measurement  circuits  provide  multiplications,  integration,  and  A/D  conversions  as  required 
(discussed  in  later  sections),  and  foward  digitized  statistical  data  to  the  digital  circuits  used 
to  calculate  the  complex  weights.  The  calculated  weight  estimates  are  applied  to  the  weighting 
elements  of  each  channel,  and  the  weighted  results  from  each  channel  arc  summed  to  form  a 
composite  signal.  The  on  time  band  is  filtered  and  output  to  the  host  radio  (after  coherent 
upconversion  to  RF),  and  is  fedback  to  the  on  time  statistical  measurement  circuits  for 
measuring  input  to  output  cross-corrlation  data.  In  addition,  a  parallel  feedback  path  is 
provided  for  off  time  data  collection.  Power  monitors  are  provided  in  both  feedback  channels 
to  allow  measurement  of  the  output  powers  in  both  on  and  off  time  bands. 

Comments  regarding  this  architecture  - 

o  Inverting  a  complex  matrix  such  as  R0  is  time  consuming.  If  the  algorithm  being  used 
requires  R0~\  only  a  few  (perhaps  just  one)  updates  of  the  weight  estimates  may  be  available 
in  a  timeslot.  One  way  around  this  limitation  is  to  delay  applying  the  off  time  data  until  the 
following  timeslot,  which  is  essentially  anticipative  processing. 

o  Since  the  on  and  off  time  signals  are  not  processed  by  the  same  circuits,  errors  will 
be  introduced:  measurements  are  made  over  two  frequency  bands  selected  by  different  filters, 
the  signals  are  processed  by  different  feedfoward  and  feedback  circuits,  and  the  actual 
measurements  are  performed  by  different  circuits. 
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FIGURE  3.  BLOCK  DIAGRAM,  SPECTRAL  SEPARATION  FOR  OFF  TIME  DATA 
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4.3.3  ARCHITECTURE,  DATA  COLLECTION  (continued) 

Comments  regarding  spectral  separation  (continued)  - 

o  The  off  time  data  may  not  include  the  same  jamming  signals  as  the  on  time  data  since 
they  are  from  different  frequency  bands.  For  the  same  jammers  to  be  present  over  adjacent 
frequency  bands  simultaneously  requires  an  assumption  that  the  jammers  are  wideband. 

o  Using  data  collected  from  different  frequency  bands  to  perform  weight  estimate 
calculations  may  introduce  frequency  dependent  errors.  If  operating  with  a  radio  network 
having  a  large  message  bandwidth  (such  as  EPLRS  which  uses  DSSS),  some  form  of  frequency 
compensation  may  be  required. 


Anticinative  - 

A  block  diagram  showing  temporally  separated  collection  of  on  and  off  time  data  is  shown  in 
figure  4.  Two  downconversions  are  required  in  this  case;  one  is  tuned  to  the  present  frequency 
being  received  by  the  host  radio,  the  other  is  tuned  to  the  next  frequency  to  be  received  by  the 
host  radio.  The  latter  sample  provides  off  time  data  in  advance  of  the  on  time  data  collect 
(anticipative).  Since  these  frequency  conversions  are  done  simultaneously  in  parallel,  two 
frequency  synthesizers  are  required. 

Figure  4  shows  that  the  majority  of  circuits  are  duplicated;  although  they  do  not  include  all 
the  circuits  required  to  be  stand  alone  AAPs,  they  are  labeled  "ODD  AAP"  and  "EVEN  AAP". 
An  example  of  which  circuits  are  used  in  what  capacity  for  a  given  timeslot  is  described: 
Assume  the  synthesizer  associated  with  the  odd  AAP  is  tuned  to  correspond  to  the  host  radio’s 
present  receive  frequency.  The  even  AAP  is  therefore  tuned  to  the  next  frequency  to  be 
received  by  the  host  radio.  In  this  case,  the  odd  AAP  is  processing  on  time  data,  the  even  AAP 
is  collecting  off  time  data  in  the  frequency  band  to  be  received  in  the  next  timeslot.  The  even 
AAP  collects  a  long  sample  of  off  time  data,  forms  the  matrix  Rq,  performs  the  calculations 
necessary  to  form  R0_1,  and  stores  this  data  for  use  during  the  on  time  to  follow  in  the 
subsequent  timeslot.  The  odd  AAP,  having  its  R0_1  available  from  the  last  timeslot,  collects  a 
block  of  on  time  data  (perhaps  16  or  32  bits  long)  by  measuring  R^j,  then  calculates  the  next 
weight  estimates  using  the  equation  (or  one  very  similar):  w(m+l)  ■  R0JRnv  Since  R^  ■ 
Rtw(m),  the  new  weight  estimates  are  updates  of  the  old  estimates,  but  are  calculated  without 
requiring  forming  the  entire  covariance  matrix  Rt  A  vector  measurement  (R^)  is  used  in 
place  of  a  matrix  measurement  (Rt)  and  matrix  operation,  thereby  reducing  the  processing  time 
for  generating  each  weight  update. 

All  the  critical  circuits  for  each  AAP  (even  and  odd)  are  provided  separately.  Each  AAP 
collects  its  own  off  and  on  time  data,  using  the  same  circuits  operating  at  the  same  frequency 
for  each  measurement  This  eliminates  many  of  the  errors  introduced  in  the  spectral  separation 
architecture.  The  weights  are  applied  and  updated  in  the  AAP  which  is  processing  on  time  data 
only;  the  resulting  composite  waveform  is  selected  by  the  output  switch,  coherently 
upconverted  and  passed  on  to  the  host  radio  at  its  RF.  The  AAP  processing  off  time  data  may 
not  use  the  weighting  elements  during  its  off  time,  but  will  in  the  subsequent  timeslot  (which 
will  be  its  on  time).  The  weighting  elements  of  the  AAP  processing  off  data  may  be  updated 
in  concert  with  the  on  time  weight  updates  to  provide  a  real  time  estimate  of  the  output 
interference  plus  noise  power;  this  may  be  useful  to  the  adaptation  algorithm  or  may  be  used 
for  a  coarse  verification  of  adaptation  performance.  Other  blocks  in  figure  4  are  similar  to 
blocks  in  figure  3  and  are  described  in  the  "Spectral  Separation"  paragraphs  above. 
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FIGURE  H.  BLOCK  DIAGRAM,  ANTICIPATIVE  OFF  TIME  DATA 
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4.3.3  ARCHITECTURE,  DATA  COLLECTION  (continued) 

Comments  regarding  the  anticipative  architecture  - 

o  Simultaneous  processing  of  on  and  off  time  data  from  different  RF  bands  requires 
duplicate  hardware  (compared  with  the  spectral  separation  approach),  including  a  synthesizer, 
frequency  conversions,  and  complex  weights. 

o  Providing  dedicated  circuits  for  each  AAP,  particularily  complex  weighting, 
downcon  version  and  IF  processing,  and  statistical  measurements,  allows  independent  adaptation 
and  eliminates  errors  which  result  in  attempting  to  share  these  functions  between  two  AAPs. 

o  The  on  and  off  time  data  is  derived  from  the  identical  frequency  band,  therefore 
frequency  compensation  is  not  required. 

o  The  on  time  is  split  into  several  blocks  over  which  data  samples  are  accumulated,  then 
weight  calculations  and  updates  are  provided.  The  off  time  has  the  same  (or  more)  time 
available,  but  only  one  set  of  statistical  measurements  are  required  (Ro).  This  allows  for  longer 
averaging  and  time  to  invert  this  matrix  (R0_1  is  usually  required). 

o  Since  the  off  time  statistics,  in  the  form  of  R0"\  are  stored  and  ready  for  use  in  the 
adaptation  equation,  and  since  the  on  time  measurements  require  the  vector  R^  only  rather 
than  a  matrix,  the  on  time  calculations  are  quick  and  simple,  allowing  multiple  iterations 
within  the  same  timeslot. 


Combination  of  Approaches  - 

It  is  possible  to  use  both  spectral  and  temporal  separation  simultaneously  for  collecting  off  time 
data.  This  may  be  beneficial  to  the  extent  that  more  data  bits  are  collected,  thereby  providing 
a  better  average  of  the  actual  environment.  However,  since  the  off  time  allows  ample  time  for 
averaging  a  large  number  of  bits  on  its  own,  and  in  light  of  the  added  complexity  for  allowing 
both  methods,  it  is  recognized  that  a  combination  of  the  two  approaches  is  not  cost  effective. 


Selection  - 

It  is  evident  in  the  proceeding  discussions  that  the  optimal  configuration  is  an  anticipative 
architecture.  This  is  the  baseline  selected  for  f urther  architecture  development  presented  in  the 
rest  of  this  report.  The  reasons  for  this  choice  are  briefly  summarized  below: 

o  Long  off  time  allows  R0_1  calculation. 

o  Availability  of  Ro'1,  ease  of  measuring  R_j,  and  simple  calculations  allow  several 
iterations  of  the  weight  estimates  within  each  timeslot. 

o  Duplicate  circuits  add  little  risk  or  cost. 

o  Same  circuits  used  for  all  measurements,  etc.  thereby  eliminating  errors, 
o  On  and  off  time  data  from  same  frequency  band,  thereby  eliminating  errors. 
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4.3.3.3  COMPLEX  WEIGHTING 
Introduction  - 

Implementing  the  complex  weights  can  be  accomplished  a  number  of  ways.  The  primary 
variables  are  whether  they  are  realized  at  RF,  IF  or  baseband  (analog  or  digital),  and  whether 
the  magnitude  and  phase  (polar)  components  or  the  real  and  imaginary  (quadrature) 
components  are  controlled.  These  options  are  illustrated  in  figure  5;  control  of  the  magnitude 
and  phase  of  the  baseband  signals  is  not  presented  since  it  is  not  a  reasonable  option.  Also,  the 
baseband  approach  shows  analog  multiplications,  although  an  actual  system  would  probably 
digitize  these  signals  and  perform  the  calculations  digitally. 

These  options  are  presented  and  discussed  in  the  following  sections.  The  baseline  approach  for 
this  report  is  then  selected  and  justified. 

At  RF  with  Quadrature  Weighting  - 

This  option  is  shown  in  figure  5a;  the  device  which  performs  this  function  is  the  Complex 
Phasor  Modulator  (CPM).  It  consists  of  a  quadrature  power  divider,  two  linear  voltage  variable 
attenuators  (VVA),  and  a  power  combiner.  The  quadrature  divider  splits  the  signal  into  in- 
phase  (I)  and  quadrature  (Q)  components.  Each  of  these  components  are  multiplied  in 
proportion  to  the  weights  applied  to  their  control  ports.  These  weights  are  based  on  correlation 
data  and  are  calculated  in  the  digital  circuits  of  the  AAP.  The  weighted  in-phase  and 
quadrature  RF  signals  are  then  added  together  in  the  power  combiner  to  form  a  resultant  signal 
which  is  coherent  to  the  original  signal  but  can  have  any  phase  or  amplitude  relative  to  the 
original  signal.  When  combined  with  the  phase  and  amplitude  adjusted  resultant  RF  signals 
from  the  other  channels  the  composite  signals  may  be  maximized  or  minimized  for  a  given 
input  signal  (steering  of  a  beam  for  signals  of  interest  or  nulls  for  jammers). 

Comments  regarding  RF,  quadrature  implementation  -  The  critical  components  of  the  CPM  are 
the  VVAs;  parameters  of  interest  for  these  devices  include: 

o  Since  they  are  applied  at  RF,  their  operating  frequency  band  must  exceed  the  host 
radio’s  RF  bandwidth.  The  amplitude  and  phase  adjustments  of  the  resultant  signals  must  be 
consistent  regardless  of  which  RF  channel  is  being  adapted,  otherwise  the  weights  calculated 
digitally  must  be  corrected  depending  on  the  center  frequency  being  adapted. 

o  Since  the  CPMs  follow  the  front  end  circuits  of  the  AAP,  they  must  exhibit  a  large 
dynamic  range  to  avoid  degrading  the  system’s  composite  dynamic  range.  If  the  VVAs  cause 
nonlinearities  for  moderately  high  level  signals,  the  amount  of  gain  in  the  front  end  must  be 
reduced,  potentially  degrading  the  overall  noise  figure  of  the  system.  If  the  VVAs  are  very 
lossy,  their  output  signals  may  drop  to  levels  below  the  internal  noise  of  following  active 
devices,  thereby  directly  degrading  the  dynamic  range  of  the  system.  The  VVAs  must  accept 
high  level  signals  and  be  low  loss  devices. 

o  The  weights  are  calculated  digitally  then  applied  to  the  CPMs  in  an  iterative  fashion. 
The  VVAs  must  respond  to  the  weight  updates  quickly  so  they  do  not  delay  adaptation.  This 
control  voltage  to  RF  output  level  settling  time  is  often  related  to  the  allowable  input  RF 
power  levels,  so  a  design  trade-off  between  speed  and  dynamic  range  may  be  necessary. 
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FIGURE  5.  COMPLEX  WEIGHTING  OPTIONS 
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4.3.3.3  ARCHITECTURE,  COMPLEX  WEIGHTING  (continued) 

Comments  regarding  RF,  quadrature  implementation  (continued)  - 

o  The  transfer  function  of  RF  gain  versus  control  voltage  must  be  nearly  perfectly 
linear  to  allow  weights  calculated  digitally  to  be  applied  at  RF  and  provide  the  desired  results. 
This  is  very  important  and  very  difficult  to  achieve  using  commercially  available  devices.  One 
solution  to  this  problem  is  to  provide  a  custom  correction  term  to  each  VVA  control  word.  This 
is  discussed  further  in  the  calibration  section  of  this  report. 

o  To  achieve  360  degrees  of  phase  shift  thru  the  CPMs,  they  must  be  capable  of 
weighting  both  the  I  and  Q  components  with  either  positive  or  negative  weights.  This  is 
accomplished  by  the  VVAs  if  they  are  bipolar  devices.  The  magnitude  of  the  gain  thru  the 
VVA  must  be  the  same  for  a  given  negative  voltage  as  it  is  for  the  same  positive  voltage,  but 
must  introduce  a  phase  shift  of  180  degrees. 

o  The  VVAs  must  provide  a  large  range  of  attenuation;  to  support  eight  bit  digital 
weight  calculations  an  RF  power  attenuation  range  of~48  dB  is  required. 


At  RF  with  Polar  Weighting  - 

This  option  is  illustrated  in  figure  5b;  the  gain  and  phase  are  controlled  directly  and 
independently  (whereas  the  quadrature  approach  requires  controlling  both  I  and  Q  components 
to  accomplish  either  a  gain  or  phase  adjustment).  This  circuit  could  be  realized  a  number  of 
ways: 


o  The  attenuator  could  be  a  bipolar  VVA  identical  to  the  devices  described  for 
quadrature  weighting,  thereby  requiring  a  phase  shifter  with  just  180  degrees  of  adjustment 
range. 


o  The  attenuator  could  be  a  unipolar  VVA  which  controls  the  magnitude  only;  the 
companion  phase  shifter  would  provide  a  full  360  degrees  of  adjustment  range. 

o  It  is  unlikely  that  the  attenuator  could  be  a  digitally  controlled  device  due  to  the 
resolution  which  would  be  required.  For  example,  eight  bit  resolution  would  require  the 
minimum  step  size  to  be  less  than  0.034  dB  of  full  scale. 

o  The  phase  shifter  could  be  digitally  controlled,  but  the  resultant  errors  due  to  phase 
resolution  limitations  are  even  more  sensitive  than  those  due  to  amplitude  resolution 
limitations. 


Comments  regarding  RF,  polar  implementations  -  Regardless  of  the  exact  implementation  of 
this  circuit,  the  composite  attenuator  and  phase  shift  function  must  meet  the  following 
requirements: 

o  As  in  the  RF/quadrature  case,  these  devices  must  operate  over  a  frequency  band 
which  exceeds  the  host  radio’s  RF  bandwidth.  See  discussion  above. 
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4.3.3  ARCHITECTURE,  COMPLEX  WEIGHTING  (continued) 

Comments  regarding  RF,  polar  implementations  (continued)  - 

o  As  in  the  RF/quadrature  case,  since  these  devices  follow  the  front  end  circuits  of  the 
AAP,  they  must  exhibit  a  large  dynamic  range  to  avoid  degrading  the  system’s  composite 
dynamic  range.  See  discussion  above.  In  addition,  since  these  devices  operate  in  series,  they 
both  must  have  high  dynamic  range,  with  the  second  device  in  the  chain  requiring  the  higher 
intercept  point  of  the  two. 

o  As  in  the  RF/quadrature  case,  these  devices  must  respond  to  digital  control  signals 
quickly  and  not  limit  the  adaptation  time  of  the  AAP.  See  discussion  above. 

o  Both  the  magnitude  and  phase  transfer  functions  versus  control  voltage  must  be  very 
precise  to  assure  the  weights  calculated  digitally  are  achieved  at  RF.  The  phase  resolution 
required  will  be  even  more  difficult  to  achieve  than  I  and  Q  magnitude  resolution.  To  achieve 
the  required  accuracy  and  maintain  it  over  temperature,  a  correction  factor  customized  for 
each  device  in  every  channel  would  be  required. 

o  As  in  the  RF/quadrature  case,  the  total  range  of  the  attenuation  must  be  at  least  48 
dB  to  accomodate  8  bit  weight  resolution.  The  combination  of  phase  shift  range  between  the 
VVA  and  the  phase  shifter  must  be  at  least  360  degrees. 

At  IF  with  Quadrature  Weighting  - 

This  option  is  illustrated  in  figure  5c;  it  is  identical  to  quadrature  weighting  at  RF  only  it  is 
realized  at  a  lower  frequency  (IF)  and  occuts  after  dehopping  so  the  required  bandwidth  is  less. 
Most  of  the  circuit  description  and  component  performance  requirements  are  identical  to  those 
discussed  for  RF,  quadrature  weighting  above,  with  the  following  differences: 

o  Since  the  signal  at  IF  is  dehopped,  the  bandwidth  requirement  for  this  circuit  needs 
to  just  exceed  the  message  bandwidth  of  the  host  signal.  If  a  70  MHz  IF  is  selected,  and  the 
message  bandwidth  is  3  MHz  (as  for  EPLRS),  the  percent  bandwidth  is  only  4.3%.  There  may 
be  no  significant  frequency  dependent  variations  in  performance  over  this  bandwidth. 

o  Since  this  implementation  occurs  after  significant  RF  and  IF  signal  processing,  the 
AAP  system  dynamic  range  will  be  less  sensitive  to  dynamic  range  limitations  of  these  circuits. 

o  Since  implemented  at  a  lower  frequency,  the  devices  used  will  be  less  sensitive  to 
secondary  factors  such  as  impedance  matching,  dielectric  variations  in  the  circuit  cards,  cross 
coupling  of  signals,  etc.  Also,  at  lower  frequencies  there  will  be  more  commercially  available 
devices  to  choose  from,  typically  with  faster  delivery  and  less  expense. 

o  The  number  of  components  in  the  RF  input  to  RF  output  signal  path  is  greatly 
increased  in  this  case.  This  is  undesirable  since  these  more  narrowband  components  are  likely 
to  introduce  amplitude,  phase  and  delay  distortions  to  the  signal.  This  may  adversely  effect  the 
host  radio  system’s  performance. 
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4.3.3  ARCHITECTURE,  COMPLEX  WEIGHTING  (continued) 
At  IF  with  Polar  Weighting  - 
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This  option  is  illustrated  in  figure  5d;  it  is  identical  to  polar  weighting  at  RF  except  at  the 
lower  frequency,  smaller  bandwidth  IF.  The  advantages  and  disadvantages  of  IF  versus  RF 
implementations  discussed  for  the  quadrature  weighting  approach  also  holds  true  for  polar 
weighting;  similarily,  the  advantages  and  disadvantages  of  polar  versus  quadrature 
implementations  discussed  for  the  RF  approaches  also  holds  true  at  IF.  See  the  discussions 
presented  in  the  preceeding  paragraphs. 

At  Baseband  with  Quadrature  Weighting  - 

This  option  is  illustrated  in  figure  5e;  it  is  significantly  more  complex  than  the  other  options 
described  above  because  it  requires  quadrature  downconversion  to  get  to  baseband,  then  must 
be  quadrature  upconverted  to  return  to  IF.  The  only  advantage  to  this  approach  is  that  the 
weights  can  be  implemented  digitally  (if  message  bandwidth  is  small  enough),  thereby 
eliminating  the  errors  introduced  applying  precisely  calculated  digital  weights  to  RF 
components  which  are  not  ideal. 

Comments  regarding  baseband  implementation  •  Problems  associated  with  this  approach  are 
listed  below; 

o  The  amplitude,  phase  and  delay  distortions  of  the  desired  signal  by  frequency 
conversion  to  and  from  IF  would  be  much  more  severe  if  converted  all  the  down  to  baseband 
and  back  up  to  RF.  In  addition  to  the  distortion  from  more  components  with  more  narrow 
bandwidth,  the  quadrature  conversions  would  introduce  phase  errors  between  in-phase  and 
quadrature  components  of  the  signal. 

o  The  benefit  of  applying  the  weights  digitally  introduces  a  limitation  to  the  bandwidth 
of  the  message  to  be  processed.  If  spread  spectrum  signals  must  be  processed,  depending  on  the 
extent  of  frequency  spreading,  it  may  become  difficult  to  perform  digitization  of  the  signal 
at  a  sufficient  rate. 

Sfiksliiia- 

The  quadrature  approach  is  preferable  to  the  polar  approach  because  VVAs  are  more  widely 
available,  exhibit  less  variations  and  require  less  control  resolution  than  phase  shifters. 
There  are  several  manufacturers  of  integrated  circuits  that  can  achieve  the  basic  VVA 
requirements,  particularly  if  implemented  at  IF.  Phase  shifters  manufacturers  are  more 
specialized;  more  options  are  available  at  microwave  frequencies  than  VHF  or  UHF.  The 
electronically  controlled  phase  shifters  are  based  on  tuning  varactor  diodes  -  components  which 
introduce  the  following  problems:  high  control  voltages  (exceeding  24  volts)  are  often  required 
for  a  wide  range  of  phase  shift;  the  transfer  function  is  nonlinear  near  the  ends  of  the  phase 
shift  range;  there  is  significant  drift  in  phase  shift  over  temperature;  the  bandwidth  of 
operation  is  usually  limited  to  ten  percent;  and  variations  are  expected  as  a  function  of 
frequency,  even  over  this  limited  range.  Similar  limitations  can  be  identified  for  VVAs,  but 
they  are  not  as  severe  and  can  be  corrected  for  more  easily.  Finally,  the  control  resolution  of 
the  phase  shifter  would  have  to  be  four  times  more  sensitive  than  required  for  the  VVAs  used 
in  the  quadrature  approach;  that  is,  the  phase  shifters  would  require  two  more  bits  of  control 
resolution  to  achieve  the  same  accuracy  (null  depth)  as  the  VVAs. 
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4.3.3  ARCHITECTURE,  COMPLEX  WEIGHTING  (continued) 

The  choice  between  RF  and  IF  implementations  is  more  difficult  to  make.  Assuming  the 
statistical  measurements  originate  at  IF  (discussed  later),  implementing  the  complex  weights 
at  IF  would  provide  more  flexibility  than  the  RF  approach.  All  significant  AAP  functions 
would  be  included  in  modules  operating  after  downconversion  from  RF  -  the  AAP  could  be 
applied  to  many  different  radio  networks  by  merely  providing  the  necessary  frequency 
down/up  conversions  and  timing  data  for  specific  applications.  More  options  exist  at  IF  than 
RF  for  the  components  to  be  used  in  the  weighting  circuits,  but  components  are  readily 
available  for  implementations  below  500  MHz.  The  deciding  factor,  however,  will  depend  on 
the  host  radio  system  the  AAP  is  to  be  deployed  with.  The  additional  delays  caused  by  IF 
processing  may  degrade  the  radio  system’s  time  of  arrival  calculations,  a  critical  parameter  in 
position  location  systems  such  as  EPLRS.  The  amplitude  and  phase  distortions  introduced  by 
IF  processing  which  are  not  required  in  RF  processing  may  degrade  the  radio  system’s 
performance  to  such  a  degree  that  IF  weighting  is  prohibited.  Until  the  specific  host  radio  is 
selected,  and  these  concerns  are  proven  to  be  negligible,  RF  weighting  is  selected  as  the  baseline 
since  it  appears  to  present  the  lower  risk  to  the  composite  (radio  /  AAP)  system  performance 


4.3.3.4  STATISTICAL  MEASUREMENTS 
Introduction  - 


The  information  required  to  calculate  the  weights  which  result  in  the  forming  of  nulls  in  the 
directions  of  jammers  while  simultaneously  forming  a  beam  in  the  direction  of  a  signal  of 
interest  is  the  correlation  data  between  each  of  the  antenna  elements  of  the  array.  This  data 
is  used  to  form  the  covariance  matrix  R^  discussed  previously  in  other  sections.  Also  of 
interest  is  the  correlation  between  each  of  the  antenna  element  inputs  and  the  AAP  output, 
which  may  be  measured  directly  as  the  cross-correlation  vector  R  or  calculated  knowing  the 
covariance  matrix  and  the  weight  vector  from  R_  =  R„w.  The  measurement  of  these 
correlation  statistics  is  the  topic  of  this  section. 


T  e  basis  of  these  measurements,  the  correlation  function,  requires  common  mathematical 
operations  to  be  performed  on  the  signals  present  at  the  AAP  inputs  and  output  - 
multiplication,  integration,  addition  and  subtraction.  These  functions  are  not  optional,  they  are 
required.  At  what  frequency  they  are  implemented  is  the  option  to  be  selected  in  this  section. 
The  options  are  RF,  IF  or  baseband,  or  combinations  of  these. 

RF  Approach  - 

The  RF  option  is  eliminated  immediatelv  because  the  statistics  arc  least  corrupted  when  made 
only  over  the  bandwidth  of  the  messat  To  accomplish  this  at  RF  would  require  a  very 
precisely  tuned  preselector  filter  which  ollows  the  hops  of  the  host  radio  system  (and  is 
therefore  incredibly  fast),  and  has  extremely  steep  stopband  roll-off  performance.  Such  devices 
do  not  exist.  The  superhet  receiver  solves  this  problem  by  downconverting  the  message 
bandwidth  to  a  fixed  IF  f requency;  this  allows  highly  selective  filtering  and  requires  only  that 
the  frequency  synthesizer  be  fast  enough  to  follow  the  hopping  radio.  This  downconversion 
scheme  is  employed  in  the  AAP  whether  the  statistics  are  measured  at  IF  or  baseband. 
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4.3.3  ARCHITECTURE,  STATISTICAL  MEASUREMENTS  (continued) 

If  RF  statistical  measurements  were  attempted,  all  the  jammer  power  present  over  the  entire 
RF  operating  band  would  be  included,  as  well  as  the  internal  noise  of  the  AAP  within  the  RF 
band.  The  message,  however,  would  be  contained  entirely  within  its  bandwidth,  which  is  a 
fraction  of  the  RF  bandwidth.  Obviously,  this  results  in  a  very  dramatic  reduction  in  the 
Signal  to  Interference  plus  Noise  Ratio  (SINR).  This  approach  would  also  require  all 
components  to  be  capable  of  operation  over  the  RF  bandwidth  -  a  much  more  difficult 
requirement  than  operation  over  the  relatively  narrow  message  bandwidth. 

IF  Approach  - 

This  option  is  selected  if  the  multiplication  occurs  at  IF;  the  desired  output  of  the  multipliers 
is  at  baseband,  so  low  pass  filters  are  employed  to  separate  the  baseband  signals  from  the 
signals  p*  twice  the  input  frequency.  These  low  pass  signals  are  then  integrated  to  acquire  the 
desired  statistical  information.  This  IF  approach  is  illustrated  in  the  block  diagram,  figure  6a. 

Comments  regarding  IF  statistical  measurements  * 

o  Compared  to  implementation  at  baseband,  the  IF  approach  requires  only  half  as  many 
multiplications.  This  is  because  each  IF  signal  includes  in-phase  and  quadrature  components, 
whereas  each  baseband  signal  has  its  in-phase  and  quadrature  components  separated,  essentially 
forming  two  separate  signals.  This  is  evident  in  the  following  equations: 

^  At  IF  - 

(xk)(y)  «=  (xkicos(wet)  -  xkqsin(wct))(yiCOS(wct)  -  yqsin(wct)) 

After  low  pass  filtering,  (xk)(y)  «=  ((xki)(yj)  +  (xkq)(yq))/2 

Similarly,  (xk)(y’)  •=  ((xki)(yq)  -  (xkq)(yi»/2 

where  y’  is  the  Hilbert  transform  of  y 

i>  At  baseband  - 

(xk+)(y+)  *  (xki  +  jxkq)(yi  +  jyq) 

*  (x^y;  +  xkqyq)  +  j(xkqyj  -  xkjyq) 

>  Comparison  of  the  above  calculations  reveals  the  following  - 

After  filtering,  the  IF  multiplication  (xk)(y)  =  Re[(xk+)(y+)at  baseband]/2 
and  (xk)(y’)  *=  -  Im[(xk+)(y+)  at  baseband]/2 

o  Both  multiply  and  integrate  outputs  are  used  to  determine  the  desired  correlation  of 
one  antenna  element  input  to  the  output.  This  is  demonstrated  in  the  following  equation: 
E[xky*]  -  2(E[xky]  -  jE[xky’] 

o  The  number  of  IF  multiplications  required  to  calculate  the  input  to  output  cross- 
correlation  vector  R,^  is  equal  to  2m,  where  m  is  the  number  of  antenna  array  elements.  A  four 
clement  array  would  require  8  multipliers,  which  is  a  reasonable  number. 

o  The  number  of  IF  multiplications  required  to  calculate  the  input  covariance  matrix 
R^  is  equal  to  2m2  +  m.  The  same  four  element  array  would  take  36  multipliers,  a  quantity 
which  would  be  difficult  to  achieve  with  reasonable  size,  cost  and  power  dissipation. 
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6.a.  IF  APPROACH 


6.b.  BASEBAND  APPROACH  (W/DIGITAL  CALCULATIONS) 


6.C.  COMBINATION  IF  &  BASEBAND 
FIGURE  6.  STATISTICAL  MEASUREMENT  OPTIONS 
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4.3.3  ARCHITECTURE,  STATISTICAL  MEASUREMENTS  (continued) 

Comments  regarding  IF  statistical  measurements  (continued)  - 

o  Based  on  the  number  of  multiplications  required,  if  only  R^  were  needed  an  IF 
scheme  would  be  very  attractive.  Only  R^  is  needed  for  some  algorithms,  but  others  (including 
Dominant  Mode  Block  Update)  require  R^  to  maintain  the  flexibility  to  implement  a  variety 
of  algorithms  using  the  same  hardware,  limiting  statistical  measurements  to  IF  is  prohibitive. 

Baseband  Approach  - 

This  approach  requires  quadrature  downconversion  to  baseband,  and  is  illustrated  in  figure  6b. 
Once  at  baseband,  the  multiplication  and  integration  can  be  done  with  analog  or  digital 
processing.  To  be  achieved  digitally,  the  baseband  signal  would  need  to  be  sampled  at  a  rate 
exceeding  twice  the  message  bandwidth,  then  a  large  set  of  samples  would  have  to  be  summed 
and  divided  by  the  number  of  samples  to  determine  the  integrated  result.  This  implies  a  large 
number  of  digital  calculations  at  potentially  high  clock  rates,  depending  on  the  bandwidth  of 
the  message  spectrum  after  spreading.  The  analog  approach  would  perform  all  functions  thru 
integration,  then  sample  the  resulting  integrated  output.  This  requires  a  much  lower  sample  rate 
and  relieves  the  speed  requirements  of  the  digital  circuitry.  The  analog  calculation  circuitry 
is  presented  in  block  diagram  form  in  figure  7. 

Comments  regarding  baseband  statistical  measurements  - 

o  Quadrature  downconversion  of  each  antenna  element’s  signal  is  required.  The  circuits 
required  to  perform  this  downconversion  are  identical  to  the  circuits  used  to  perform  the  IF 
multiplications  discussed  above. 

o  All  the  information  required  to  calculate  the  elements  of  the  covariance  matrix  is 
available  from  the  baseband  data.  The  baseband  data  is  of  the  form  xk  *  xki  +  xkq.  If  formed 
digitally,  the  covariance  between  channels  k  and  m  requires  integration  of  the  following 
calculation:  xkxm*  «  (xkiXmi  +  xkqxmq)  +  jU^x^  -  xkixmq). 

o  Having  all  the  covariance  matrix  data  available  allows  for  the  most  flexibility  in 
algorithm  selection.  Few  algorithms  can  be  used  if  only  vector  measurements  are  available. 

o  The  speed  required  of  the  digital  calculations  is  greatly  reduced  by  performing  the 
multiply  and  integrate  functions  in  analog  circuits.  The  integration  will  be  performed  over  a 
period  corresponding  to  at  least  16  data  samples;  for  EPLRS,  having  a  chip  rate  of  about  3 
MHz,  this  requires  a  digital  sampling  rate  of  once  every  5.33  microseconds  (187.5  KHz).  The 
same  EPLRS  data,  if  sampled  at  baseband  before  integration,  would  need  to  be  sampled  at  least 
twice  during  each  chip,  or  at  least  every  167  nanoseconds.  To  avoid  aliasing,  this  rate  would 
need  to  be  higher  than  this;  assume  every  125  nanaseconds  (8  MHz).  The  sampling  rate  of 
EPLRS  statistics  is  therefore  43  times  faster  if  integration  is  done  digitally.  For  systems  which 
do  not  use  DSSS,  the  factor  of  43  still  holds  true,  but  the  sample  rates  would  be  much  slower 
(and  within  reason)  since  the  message  bandwidth  is  so  much  less. 

o  The  precision  of  the  digital  calculations  is  superior  to  either  IF  or  baseband  analog 
methods.  In  any  case,  a  calibration  sequence  will  be  needed;  this  is  described  later  in  a  section 
dedicated  to  calibration  issues,  section  4.3. 3. 7. 
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4.3.3  ARCHITECTURE,  STATISTICAL  MEASUREMENTS  (continued) 
Combination  IF  &  Baseband  - 
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As  discussed  above,  both  the  IF  and  baseband  approaches  offer  advantages  for  accumulating 
different  sets  of  data.  All  the  matrix  data  is  available  using  baseband  measurements,  but  this 
is  slower  and  more  cumbersome  than  IF  vector  measurements.  Realizing  that  many  algorithms 
require  the  inverted  covariance  matrix  for  off  time  data  only,  and  that  several  iterations  of 
weight  estimates  are  desired  during  on  time  processing,  it  is  apparent  that  both  techniques 
could  be  employed  in  an  optimal  system,  providing  the  added  circuit  complexity  is  not 
prohibitive. 

A  means  of  providing  fast  statistical  vector  measurements  at  IF  and  complete  statistical  matrix 
measurements  more  slowly  at  baseband,  with  the  addition  of  only  an  IF  switch,  is  illustrated 
in  figure  6c.  When  off  time  data  is  to  be  collected,  the  IF  switch  selects  the  fixed  oscillator  as 
the  input  to  the  multipliers.  The  quadrature  power  splitter  provides  the  cosine  and  sine  terms 
required  for  quadrature  downconversion  to  baseband.  As  described  above,  this  baseband  data 
can  be  processed  with  analog  or  digital  circuits  for  forming  the  desired  statistics.  The  entire 
timeslot  may  be  used  for  calculating  the  off  time  statistics,  since  the  weights  are  not  iterated 
during  this  period.  The  last  portion  of  the  off  timeslot  could  be  used  to  invert  the  matrix 
formed  from  off  time  measurements,  then  the  R0_1  matrix  will  be  ready  for  use  in  calculating 
updates  of  the  weights  during  on  time  processing. 

The  on  time  data  requires  only  the  cross-correlation  vector;  by  using  the  IF  switch  to  select  the 
output  (feedback)  IF  signal,  the  multipliers  have  the  input  signals  required  to  perform  these 
measurements.  In  this  case,  the  quadrature  power  splitters  after  the  IF  switch  provide  the 
output  y  and  the  Hilbert  transform  of  y  which  are  required  for  cross-correlation  measurements. 
The  on  time  data  can  be  collected  and  integrated  in  blocks  of  16  samples  or  more,  and  an 
update  of  the  weight  estimates  can  be  calculated  based  on  each  of  these  blocks.  The  calculations 
can  be  performed  quickly  in  the  digital  calculator,  since  R0''  is  stored  in  memory  and  R^  is 
updated  every  16  chip  samples  or  so.  In  the  Dominant  Mode  Block  Update  algorithm,  the 
equation  being  implemented  is  w(m+l)  = 

Selection  - 

Based  on  the  discussion  above,  the  obvious  baseline  choice  is  to  implement  the  statistical 
measurements  both  at  IF  and  baseband  since  the  only  added  complexity  is  an  IF  switch. 


4.3.3.S  WEIGHT  CALCULATIONS 
Introduction  - 

To  be  compatible  with  multiple  algorithms,  the  architecture  must  provide  for  digital 
calculation  of  the  complex  weights.  If  a  specific  algorithm  is  to  be  implemented,  and  the 
hardware  can  be  customized  for  that  particular  type  of  algorithm,  then  a  non-digital  approach 
may  be  justified.  For  general  applications,  such  as  the  type  addressed  in  this  study,  digital 
weight  calculation  is  mandatory. 

The  digital  implementation  is  described  below,  with  indication  of  its  power  and  flexibility.  An 
example  of  a  non-digital  approach,  which  has  been  employed  and  proven  successful,  is  then 
introduced.  The  limitations  of  this  technique  are  discussed  and  revealed  to  be  too  restrictive. 
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4.3.3  ARCHITECTURE,  WEIGHT  CALCULATIONS  (continued) 
Digital  Approach  - 
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Calculating  the  weights  digitally,  especially  using  a  microprocessor  based  design,  provides  a 
great  deal  of  flexibility  to  the  AAP  system.  The  inputs  to  the  weight  calculations  are  the 
digitized  correlation  measurements  and  some  power  level  measurements  made  at  IF.  Equations 
are  used  to  manipulate  this  data  and  generate  estimates  of  the  optimal  weights;  the  equations 
used  and  the  details  of  the  manipulations  are  dictated  by  the  algorithm  selected.  Different 
calculations,  all  using  the  same  data,  can  be  performed  in  a  microprocessor  based  system  by 
installing  different  firmware  programs.  Since  most  algorithms  use  the  same  data,  but  perform 
different  calculations,  reprogramming  the  firmware  allows  implementing  different  algorithms. 
The  outputs  are  the  digital  representations  of  the  complex  weight  estimates  for  each  antenna 
element;  these  will  be  converted  to  analog  voltages  within  the  CPM  assemblies  for  application 
to  the  WAs’  control  ports.  Since  the  correlation  data  available  digitally  includes  both  on  and 
off  time  data,  the  AAP  will  steer  nulls  in  the  direction  of  jammers  while  simultaneously 
steering  a  beam  in  the  direction  of  the  desired  signal. 

Additional  benefits  of  using  digital  weight  calculations:  I)  Including  digital  circuits  between 
the  IF  measurements  and  the  VVA  controls  allows  the  weights  calculated  at  each  timeslot  to 
be  stored;  this  allows  adaptation  for  each  timeslot  to  begin  with  weights  which  were  optimized 
at  the  last  visit  to  that  frequency  bin.  2)  Digital  circuits  can  be  implemented  on  a  monolithic 
chip  which  requires  very  little  space  and  dissipates  minimal  power.  These  features  could  be 
exploited  if  development  proceeds  to  full  scale  production  activities. 

Use  of  a  microprocessor  based  system  for  calculating  the  weights  allows  flexibility  in  other 
aspects  of  the  AAP  as  well.  The  same  microprocessor  could  be  used  for  interfacing  to  external 
devices  and  controlling  other  internal  functions  besides  the  complex  weights.  A  block  diagram 
of  a  digital  controller  capable  of  such  functions  is  presented  in  figure  8;  the  basic  functions 
of  this  controller  are  discussed  under  the  paragraph  "Digital  Controller",  part  of  the  "Baseline 
AAP  Architecture"  section  (4.3.3.6). 

Non-digital  Approaches  - 

One  example  of  a  non-digital  approach  is  the  hybrid  analog/digital  approach  described  below: 
An  LMS  steepest  descent  algorithm  may  be  selected  which  attempts  to  minimize  the  output 
power  of  the  AAP  when  it  was  known  that  no  signals  of  interest  were  present.  The  only 
statistical  data  required,  using  the  terminology  associated  with  this  study,  is  the  off  time  cross- 
correlation  vector  data.  This  measurement  could  be  performed  at  IF,  resulting  in  an  analog 
voltage  proportional  to  the  degree  of  correlation  between  the  output  error  and  each  input 
channel.  This  analog  output  (rather  than  being  digitized  for  weight  calculations)  is  used  to 
generate  a  digital  bit  stream  whose  duty  cycle  is  proportional  to  the  voltage  level.  This  bit 
stream  provides  the  up/down  control  to  a  digital  counter;  the  output  of  the  digital  counter  is 
the  weighting  of  one  component  (I  or  Q)  of  one  of  the  antenna  elements’  signals.  If  the 
component  of  a  particular  complex  weight  is  optimal,  the  correlation  between  the  output  error 
and  that  channel’s  input  would  be  zero,  resulting  in  a  fifty  percent  duty  cycle,  causing  the 
control  word  to  the  weighting  element  to  toggle  by  its  LSB  about  its  optimal  weight.  Similarily, 
if  there  is  a  portion  of  the  output  error  which  correlates  with  the  input  of  that  channel,  the 
duty  cycle  would  not  be  fifty  percent,  so  the  counter  output  would  either  increase  or  decrease, 
thereby  changing  the  weight  in  a  way  which  will  reduce  the  output  power. 
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4.3.3  ARCHITECTURE,  WEIGHT  CALCULATIONS  (continued) 
Non-digital  Approaches  (continued)  - 
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The  hybrid  analog/digital  approach  has  been  used  and  worked  precisely  as  it  was  designed. 
This  approach  has  limitations,  however.  First,  since  an  algorithm  based  on  the  minimization 
of  output  power  in  the  presence  of  interference  plus  noise  only  is  used,  the  architecture  is 
optimized  for  off  time  cross-correlation  data  collection  only.  Any  adaptive  processing 
performed  using  off  time  data  only  can  steer  only  nulls;  the  signal  of  interest  may  be  arriving 
at  an  angle  at  which  the  AAP  has  inadvertently  steered  a  null,  resulting  in  unintentional 
degradation  in  output  SINR.  By  allowing  for  on  time  data  collection  and  processing  as  designed 
and  simulated  in  this  study,  a  beam  is  always  provided  in  the  direction  of  the  desired  signal. 
Secondly,  the  weight  updates  are  provided  by  a  closed  loop  integration,  without  dedicated 
calculations.  Alternative  equations  for  calculating  and  updating  the  weights  are  not  provided 
for.  The  digital  calculations  provided  by  this  study’s  architecture  allows  much  more  flexibility 
in  selecting  and  optimizing  adaptation  algorithms. 


4.3.3.6  BASELINE  AAP  ARCHITECTURE 
Introduction  - 

This  section  begins  by  compiling  all  the  architectural  choices  made  in  the  preceeding 
discussions  (including  RF  weighting)  to  form  our  first  choice  baseline  architecture.  This 
baseline  will  be  further  developed  in  this  section  by  detailed  discussion  of  each  major 
subassembly,  and  completed  by  a  discussion  of  the  calibration  circuits  and  sequences  required 
for  optimal  performance  (section  4.3.3.7).  Since  an  IF  implementation  of  the  weights  is  such  an 
attractive  option,  it  is  also  presented  as  a  baseline,  but  is  identified  at  this  time  as  our  second 
choice.  The  discussions  regarding  further  development  of  the  architecture  is  generally 
applicable  to  both  baseline  choices,  and  includes  design  details  for  the  RF  Assembly,  the  IF 
Chains,  the  Frequency  Synthesizers,  the  Digital  Controller,  the  DC  Power  Supplies,  and  the 
Antenna  Array. 

Baseline  Choice  #1  - 

Based  on  the  selections  made  in  sections  4.3. 3.2  thru  4.3. 3. 5  above  -  off  time  data  collected 
during  anticipative  timeslots,  the  complex  weights  implemented  at  RF,  a  combination  of  IF  and 
baseband  statistical  measurements,  and  digital  weight  calculations  -  a  block  diagram  of  the 
baseline  AAP  system  can  be  drawn.  This  is  shown  in  figure  9. 

The  input  signals  come  from  each  element  of  the  antenna  array.  A  receiver  front  end 
conditions  this  signal  for  the  following  circuits,  establishing  the  AAP  noise  figure.  The  signal 
is  power  split  to  form  two  identical  signals.  Both  signals  will  undergo  on  and  off  time 
processing,  but  during  opposite  timeslots  from  one  another.  The  timeslots  alternate  between  odd 
and  even  timeslots;  allow  the  path  which  processes  on  time  data  during  odd  timeslots  to  be 
called  the  odd  AAP  and  the  path  which  processes  on  time  data  during  even  timeslots  to  be 
called  the  even  AAP.  This  implies  that  the  odd  AAP  processes  off  time  data  during  even 
timeslots,  and  the  even  AAP  processes  off  time  data  during  odd  timeslots.  Components  are 
grouped  and  labeled  "odd  AAP"  and  "even  AAP"  in  figure  9  to  reflect  this  convention. 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 

Both  the  odd  and  even  AAPs  split  the  RF  signal  again,  providing  identical  signals:  one  is 
applied  to  the  CPMs  for  complex  weighting,  the  other  is  downconverted  to  IF  for  conditioning 
prior  to  being  applied  to  the  statistical  measurement  circuits.  A  chain  of  IF  components  provide 
this  signal  conditioning. 

During  odd  timeslots,  the  output  of  each  odd  AAP  CPM  is  a  phase  and  amplitude  adjusted 
version  of  the  signal  present  at  its  antenna  element;  a  composite  signal  is  formed  by  combining 
all  those  weighted  RF  signals.  When  properly  adapted,  this  signal  will  have  the  maximum 
attainable  SINR,  meaning  the  jammers  have  been  minimized  by  steering  nulls  in  their 
directions  while  the  desired  signal  has  been  maximized  by  steering  a  beam  in  its  direction.  This 
composite  signal  is  power  split  to  provide  an  output  to  the  host  radio  and  a  feedback  signal  for 
further  AAP  processing.  The  feedback  signal  is  downconverted  and  processed  for  application 
to  the  statistical  measurement  circuits.  A  chain  of  IF  components  very  similar  to  those  used  for 
processing  the  feedforward  IF  signal  is  used  in  the  feedback  path  as  well.  During  odd  timeslots, 
the  statistical  measurement  circuits  in  the  odd  AAPs  perform  the  correlation  measurements 
between  the  inputs  and  the  output  directly  at  IF,  resulting  in  Rxyl- 

During  odd  timeslots,  the  even  AAPs  are  used  only  to  collect  off  time  data.  The  complex 
weights  in  these  paths  could  be  set  to  maximum  attenuation  and  ignored  for  this  timeslot,  or 
could  be  used  for  secondary  data  collection  such  as  estimates  of  Pout.  The  feedforward  signal 
is  downconverted  and  processed  by  its  IF  chain,  ultimately  being  applied  to  the  statistical 
measurement  circuits.  In  this  case,  the  output  is  not  applied  to  the  other  inputs  of  the  statistical 
measurement  circuits,  but  rather  a  sinusoidal  waveform  is  applied  via  the  feedback  IF  chain. 
This  allows  the  IF  multipliers  to  perform  a  quadrature  downconversion,  thereby  generating  the 
baseband  data  for  each  channel.  This  data  contains  all  the  information  required  to  calculate 
and  construct  the  off  time  covariance  matrix  Rp  This  data  is  accumulated  and  stored  during 
off  times  for  use  in  the  adaptation  algorithm  during  the  subsequenct  on  time. 

During  even  timeslots,  the  operation  of  the  odd  and  even  AAPs  are  reversed  from  those 
described  in  the  preceding  paragraphs.  The  odd  AAP  does  off  time  processing  while  the  even 
AAP  performs  on  time  operations.  In  order  to  provide  a  continous  RF  output  to  the  host  radio, 
the  outputs  of  the  odd  and  even  AAPs  are  alternately  selected  in  sequential  timeslots.  This 
selection  is  achieved  by  the  RF  switch  at  the  AAP  unit  output. 

In  this  top  level  block  diagram  the  digital  f unctions  are  shown  simply  as  boxes.  These  functions 
are  described  completely  in  sections  4.3. 3.5.  and  4. 3.3.6. 
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This  architecture  is  similar  to  choice  #1  except  the  complex  weighting  is  realized  at  IF  rather 
than  RF.  If  the  host  system  does  not  require  the  minimum  distortions  and  delays  provided  by 
RF  weighting,  this  baseline  would  be  more  strongly  recommended.  IF  weighting  would  enable 
all  functions  which  are  purely  AAP  related  to  be  based  on  IF  signals.  This  would  allow 
application  of  the  AAP  to  communication  systems  having  different  operating  bands  by  simply 
replacing  the  AAP  front  ends  and  synthesizers.  These  assemblies  could  be  plug  in  modules  and 
easily  interchanged.  The  block  diagram  for  this  architecture  is  shown  in  figure  10.  Only  the 
architectural  differences  caused  by  IF  rather  than  RF  weighting  will  be  discussed  below. 

For  IF  weighting,  the  RF  signal  must  be  downconverted  before  application  to  the  CPMs.  To 
accomplish  this,  the  downconverter  used  in  the  RF  weighting  approach  is  moved  in  front  of 
the  power  splitter  which  provides  a  feedforward  path  and  an  input  to  the  CPMs. 

Since  the  outputs  of  the  odd  and  even  AAPs  are  at  IF,  which  is  compatible  with  the  feedback 
signal  requirements,  the  downconverters  provided  in  the  feedback  path  of  the  RF  weighting 
approach  are  no  longer  necessary.  Instead,  to  translate  the  output  to  the  host  radio’s  RF  band, 
an  upconverter  is  required  after  the  output  switch  which  interleaves  the  even  and  odd  time 
AAP  outputs. 

RF  Assembly  - 

By  implementing  the  weighting  at  RF,  the  number  of  components  which  process  the  desired 
signal  are  minimized,  and  are  selected  to  operate  over  the  widest  reasonable  bandwidth  to 
minimize  delay  and  distortion.  A  block  diagram  of  the  RF  Assembly  is  given  in  figure  11. 

The  signals  from  each  antenna  element  is  input  to  the  equivalent  of  a  receiver  front  end:  RF 
power  limiters  for  protection,  bandpass  filters  for  establishing  the  operating  frequency  range 
(for  example,  matched  to  EPLRS,  approximately  423  to  450  MHz),  and  low  noise  amplifiers 
which  impact  the  dynamic  range  of  the  AAP  system. 

It  is  important  to  select  components  with  low  insertion  losses  for  each  component  preceding  the 
RF  amplifier,  and  a  low  noise  amplifier  must  be  used,  to  minimize  the  AAP  noise  figure.  As 
an  example,  a  20  dB  coupler  may  be  used,  giving  a  theoretical  "coupling"  loss  of  .044  dB,  and 
an  expected  total  insertion  loss  of  just  0.25  dB.  All  components  except  the  filter  will  have  an 
operating  frequency  range  greatly  exceeding  the  EPLRS  RF  bandwidth,  thereby  minimizing 
the  accumulated  phase  and  amplitude  distortion  and  insertion  delay.  The  limiter  -bandpass 
filter  will  limit  up  to  10  watts  of  incoming  RF  signal  levels  to  leakage  of  less  than  +13  dBm, 
while  providing  steep  roll-off  via  the  filter  function.  The  tradeoff  between  desirable  passband 
performance  (low  insertion  loss,  low  group  delay  peaking  near  3  dB  points,  low  phase  ripple, 
etc.)  and  required  stopband  rejection  (particularity  image  rejection)  will  be  evaluated  based  on 
their  anticipated  impact  on  system  performance.  The  RF  amplifier  will  be  selected  based  on 
its  overall  dynamic  range  capability;  the  low  end  is  determined  by  its  noise  figure,  the  high  end 
by  its  third  order  intercept  point. 

The  only  "custom"  component  in  the  RF  Front  End  is  the  limiter  -  bandpass  filter;  all  other 
components  will  be  available  as  catalog  items.  Although  the  limiter  -  bandpass  filter  will 
require  a  custom  design,  this  is  standard  procedure  for  RF  filters  and  can  be  accomplished  in 
8  weeds  without  premium  charges. 
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The  primary  function  included  in  the  RF  Assembly  is  the  implementation  of  complex 
weighting  of  each  antenna  element’s  output.  The  baseline  approach  employs  a  complex  phasor 
modulator  (CPM)  which  allows  the  quadrature  components  of  the  RF  signal  to  be  amplitude 
weighted  independently.  The  result  is  a  replica  of  the  input  signal  which  may  be  manipulated 
to  have  any  desired  amplitude  and  phase  relative  to  the  other  antenna  elements.  The  CPMs  were 
discussed  in  detail  in  section  4.3.3. 3;  they  must  include  precise  error  correction  capability  to 
assure  accurate  application  of  the  weights  calculated  in  the  digital  circuits  of  the  AAP  system. 
Discussion  of  this  calibration  technique  is  included  in  section  4. 3. 3. 7. 

The  outputs  of  the  CPMs  are  combined  to  form  a  resultant  signal  which  has  the  maximum 
attainable  SINR  for  the  given  RF  environment.  This  output  is  power  split  to  provide  a 
feedback  signal  used  in  calculating  the  weights,  and  an  output  which  is  switched  between  the 
odd  and  even  AAP  circuits  to  form  the  final  AAP  system  RF  output  to  the  host  EPLRS  radio. 

IF  Assembly  - 

The  IF  Assembly  will  include  the  IF  chains  (with  downconverters)  and  the  statistical 
measurement  circuitry.  A  block  diagram  of  these  functions  is  provided  as  figure  1 1. 

IF  Chains  -  The  inputs  to  this  assembly  are  at  RF  frequencies;  one  of  the  first  functions 
performed  in  the  feedback  and  feedfoward  chains  is  to  downconvert  the  message  band  from 
RF  to  a  convenient  IF  (such  as  70  MHz).  This  downconversion  is  accomplished  at  the  mixers 
of  the  IF  Chains;  although  the  choice  of  a  mixer  stems  insignificant,  its  selection  is  critical  to 
overall  AAP  system  performance.  The  mixer  may  generate  many  undesired,  spurious  signals 
which  could  fall  within  the  3  MHz  IF  passband,  causing  erroneous  correlation  measurements. 
Sources  of  these  spurious  signals  include  mixer  cross-products  and  two  tone  intermodulation 
products.  To  minimize  these  spurs,  a  termination  insensitive  mixer  could  be  used.  These  include 
additional  circuitry  which  reduces  the  mixer’s  dependence  upon  things  such  as  drive  level,  IF 
port  termination  or  RF  source  impedances,  and  provides  better  impedance  matches  for 
interfacing  components.  The  mixers  must  be  selected  to  assure  they  do  not  limit  the  high  end 
of  the  AAP  system  dynamic  range  because  of  an  inadequate  third  order  intercept  point.  Higher 
LO  drive  mixers  may  be  used  to  relieve  this  constraint. 

The  bandpass  filter  preceding  the  mixer  in  the  IF  chain  is  to  reject  the  LO  leakage  from 
flowing  in  the  opposite  direction  of  the  intended  signal  path  and  being  radiated  via  the 
antenna  elements.  Although  this  is  not  a  critical  component,  it  may  prove  necessary  to  meet 
FCC  regulations  for  spurious  transmissions. 

Selection  of  the  LO  and  IF  frequencies  must  be  supported  by  a  spur  analysis  which  will 
indicate  which  mixer  spurs  will  fall  inside  the  IF  passband.  The  spurs  must  be  kept  to  levels 
which  are  insignificant  compared  to  the  expected  SOI  levels;  this  is  difficult  when  jammers 
exceeding  the  desired  signals  by  30  dB  or  more  are  present,  and  are  therefore  determining  the 
high  end  of  the  dynamic  range.  Careful  selection  of  the  operating  frequencies  and  components 
will  optimize  this  difficult  performance  requirement.  The  IF  passband  and  stopband 
performance  will  be  determined  by  the  SAW  filters. 
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FIGURE  12.  IF  ASSEMBLY 


67 


4.3.3  ARCHITECTURE,  BASELINE  (continued) 

IF  Assembly.  IF  Chains  (continued)  - 

Surface  Acoustic  Wave  (SAW)  filters  are  used  to  reject  all  signals  except  those  within  the 
message  bandwidth  (as  an  example,  the  EPLRS  message  has  a  3  MHz  bandwidth).  SAW  devices 
are  used  here  because  they  provide  very  steep  stopband  rolloff  while  providing  very  linear 
insertion  phase.  The  only  disadvantage  to  using  SAW  filters  is  their  insertion  loss;  losses  of 
about  25  dB  are  expected.  The  impact  of  this  high  loss  can  be  minimized,  however,  by  careful 
signal  level  management  to  optimize  the  chain’s  dynamic  range.  The  amplifiers  on  either  side 
of  the  filter  help  accomplish  this,  and  provide  a  buffer  between  the  SAW  device  and  other 
circuit  components  (SAW  devices  generally  provide  po  >r  input  and  output  impedances). 

The  level  of  the  IF  signals  are  adjusted  using  digitallv  controlled  step  attenuators  to  fall  within 
the  most  linear  operating  range  of  the  key  components  in  the  correlator  assemblies:  the  four 
quadrant  multipliers.  The  step  attenuators  must  have  at  least  32  dB  of  range  and  respond  to  the 
digital  command  in  less  than  one  microsecond.  Control  transient  signals  leaking  through  this 
device  to  its  output  must  not  fall  in  the  IF  passband;  this  can  be  avoided  by  proper 
specification  and  design  of  the  step  attenuator. 

The  output  of  the  IF  chain  is  coupled  off  to  a  power  monitor  circuit.  This  circuit  will  provide 
a  digital  output  word  proportional  to  the  IF  chain’s  output  power  level.  This  word  may  be  used 
by  the  adaptation  algorithm  in  calculation  of  the  next  weight  estimate,  or  just  used  to  adjust 
the  step  attenuator.  The  power  monitor  circuit  could  be  designed  and  constructed  using  catalog 
ICs  to  form  a  successive  detection  logarithmic  IF  strip. 

The  last  component  in  an  IF  chain  is  a  delay  line  which  must  be  adjusted  at  subsystem 
integration  to  align  the  delays  of  each  feedfoward  path  relative  to  one  another  to  within  a 
fraction  of  one  period  at  the  IF  frequency.  For  example,  at  70  MHz  one  cycle  takes  about  14.3 
nanoseconds.  The  delays  can  be  adjusted  to  within  less  than  one  nanosecond,  or  less  than  30 
degrees.  Using  .085"  diameter  semirigid  cable,  a  miniature  spool  of  coax  cable  could  provide 
ample  adjustment  range  and  resolution  in  a  reasonable  volume.  The  coax  exhibits  about  1 
nanosecond  of  delay  per  foot. 

Statistical  Measurement  Circuits  - 

The  key  components  of  these  circuits  are  the  four  quadrant  multipliers.  These  multipliers 
operate  at  IF,  and  perform  one  of  two  functions:  1)  For  input  to  output  cross-correlation 
measurements  the  multipliers  are  simply  that  -  multipliers.  Their  output  in  this  case  is  "x  times 
y";  after  low  pass  filtering  and  integration  this  becomes  the  expected  value  of  x  times  y,  or  the 
cross-correlation  between  the  input  and  the  output.  2)  For  correlation  measurements  between 
the  four  antenna  elements,  these  multipliers  perform  quadrature  downconversion  from  IF  to 
baseband.  This  is  realized  by  switching  the  feedback  signal  from  the  output  signal  to  an 
internally  generated  CW  signal  at  the  IF  frequency.  These  baseband  signals  can  then  be  used 
to  calculate  all  the  expected  values  required  to  form  the  covariance  matrix,  using  either  analog 
or  digital  processing. 

Cross-correlation  measurements  require  that  the  feedback  IF  (which  contains  the  same  narrow 
band  signals  as  the  RF  output)  be  correlated  with  each  feedfoward  IF  (same  narrow  band 
signals  as  each  input).  To  enable  this  correlation,  the  feedback  IF  switch  is  used  to  select  the 
signal  coming  f rom  the  IF  feedback  chain.  This  signal  is  power  split  to  provide  four  identical 
signals,  one  for  correlation  with  each  of  the  feedfoward  signals.  Each  cha-'^l’s  operation  is 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 

IF  Assembly.  Statistical  Measurements  (continued)  - 
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identical.  The  feedback  IF  is  quadrature  split  to  provide  I  and  Q  representations.  The 
feedfoward  IF  is  power  di  /ided  to  create  two  identical  signals  for  multiplication  with  the  I  and 
Q  feedback  representations.  The  IF  multiplication  is  performed  by  a  linear  four  quadrant 
multiplier  IC  which  is  available  from  multiple  sources.  The  dynamic  range  of  these  multipliers 
are  limited,  however,  which  is  why  step  attenuators  are  used  in  the  feedfoward  and  feedback 
paths  to  adjust  their  levels.  The  step  attenuator  settings  are  taken  into  account  in  the  digital 
calculators.  The  outputs  of  the  multipliers  are  low  pass  filtered  to  reject  the  spectrum  centered 
at  twice  the  IF  frequency.  These  signals  are  then  integrated  (which  calculates  the  expected 
value  of  input  times  output,  I  and  Q)  and  sampled  as  dictated  by  the  controls  provided  by  the 
AAP  microprocessor.  The  sampled  data  is  held  in  registers  until  retrieved  by  the 
microprocessor. 

Covariance  measurements  use  the  same  circuits,  but  with  different  functions.  The  feedback  IF 
switch  is  used  to  select  a  precision  oscillator  rather  than  the  output  IF  signal.  When  quadrature 
split  and  multiplied  with  the  feedfoward  signal  the  signals  out  of  the  lowpass  filters  are  the 
I  and  Q  baseband  components  of  the  message  bandwidth.  The  availability  of  these  baseband 
signals  allow  calculation  of  any  covariance  terms  desired;  this  allows  calculation  of  the 
covariance  matrix  required  by  the  baseline  algorithm  for  off  time  statistics. 

The  baseband  I  and  Q  signals  can  be  used  to  calculate  covariance  terms  either  digitally  or  with 
analog  circuits.  For  digital  calculations,  both  I  and  Q  signals  must  be  sampled  at  rate  greater 
than  twice  the  message  bandwidth.  For  message  bandwidths  of  3  MHz  (such  as  EPLRS)  this  is 
easily  realizable,  although  it  places  a  burden  on  the  digital  calculator  circuits.  The  integration 
is  done  digitally  by  an  accumulator  which  adds  the  current  sample  to  the  previous  sum,  then 
after  many  samples  divides  the  total  by  the  number  of'samples.  As  in  the  IF  correlation  case, 
these  integrated  results  are  the  I  and  Q  components  of  the  expected  values. 

If  analog  circuits  are  used  to  operate  on  the  baseband  data  to  calculate  the  covariance  statistics, 
the  high  speed  requirements  of  the  digitization  and  digital  calculations  are  relieved.  The  analog 
approach  allows  the  integration  to  be  done  before  sampling,  resulting  in  a  much  slower  sample 
rate.  The  analog  calculator  is  shown  in  figure  7.  The  circuits  need  only  support  the  message 
bandwidth  and  are  comprised  of  linear  multipliers,  analog  summers,  and  integrators.  The 
integrator  timing  and  digital  sampling  are  controlled  by  the  microprocessor  as  was  the  case  for 
IF  correlations.  The  digital  samples  are  held  in  output  registers  until  retrieved  by  the 
microprocessor. 

Regardless  of  which  correlation  method  is  used,  the  accuracy  of  the  measurement  is  crucial  to 
the  AAP  system  performance.  Calibration  circuits  must  be  used  to  provide  the  required 
accuracy;  the  calibration  circuits  and  sequences  are  described  in  section  4.3. 3. 7. 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 

Synthesizer  - 

The  Synthesizer  Assembly  responds  to  the  commands  of  the  digital  controller  (based  on 
frequency  tuning  data  from  the  host  radio)  to  generate  a  fast  hopping  Local  Oscillator  (LO) 
signal  used  in  the  superheterodyne  function  of  each  IF  chain.  The  synthesizer  must  be  designed 
to  tune  over  a  range  of  frequencies  compatible  with  the  host  radio’s  operating  frequency  range 
in  much  less  time  than  is  alotted  for  one  timeslot.  Using  EPLRS  as  an  example,  the  synthesizer 
must  tune  over  21  MHz  in  significantly  less  than  1.5  milliseconds,  preferably  less  than  60 
microseconds.  It  must  settle  to  within  less  than  5  KHz  of  its  final  frequency  in  this  time.  To 
allow  the  AAP  system  to  be  applied  to  a  variety  of  Army  ;adio  systems,  the  synthesizer  design 
should  support  much  smaller  steps  in  hop  frequencies  than  the  3  MHz  step  size  of  the  EPLRS 
network.  A  step  size  of  25  KHz  is  easily  achievable  and  would  be  compatible  with  other  Army 
radios  such  as  Sincgars.  In  a  single  phase  locked  loop  (PLL)  design,  smaller  step  size  requires 
longer  settling  time,  a  trade-off  which  can’t  be  made  in  this  system.  The  solution  is  to  employ 
a  dual  PLL  synthesizer  with  the  fine  step  size  determined  by  a  very  fast  direct  digital 
synthesizer  capable  of  DC  to  15  MHz  operation.  The  actual  range  used  might  only  be  13.75  to 
15.00  MHz.  The  coarse  tuned  PLL  has  a  fixed  1.25  MHz  reference,  which  allows  wide  loop 
bandwidths  and  correspondingly  fast  settling  time.  To  provide  very  fast  tuning,  even  over  the 
entire  operating  band,  a  coarse  tune  digital  control  will  pre-steer  both  the  VCOs  to  nearly  the 
correct  frequency,  then  the  PLLs  will  be  enabled  to  allow  accurate  and  fast  settling  to  the 
desired  frequency.  A  block  diagram  of  the  synthesizer  is  shown  in  figure  13. 

Digital  Controller  - 

As  explained  in  section  4.3.3.S,  digital  calculation  of  the  weights  using  a  microprocessor  based 
design  provides  a  great  deal  of  flexibility  to  the  AAP  system.  The  algorithm  selected  for 
calculating  these  weights  can  be  changed  by  reprogramming  the  microprocessor,  assuming  the 
rest  of  the  AAP  hardware  supports  the  algorithm  as  well.  The  digital  controller  provides  all 
interface  and  control  functions  between  the  AAP  system  and  external  systems  such  as  the  host 
radio  and  a  personal  computer.  It  also  interfaces  the  digital  assembly  with  all  other  AAP 
subsystems.  Using  instructions  and  information  from  external  sources,  it  collects  the  necessary 
data  from  the  other  internal  subsystems  to  perform  the  algorithm’s  calculations,  then  controls 
the  other  subsystems  for  null  and  beam  steering  which  optimizes  the  output  SINR.  The  weight 
calculation  requirements  of  the  digital  controller  were  discussed  previously  (4. 3. 3. 5);  this 
section  will  address  the  other  digital  controller  functions. 

The  same  microprocessor  used  for  calculating  the  weights  could  be  used  for  interfacing  to 
external  devices  and  controlling  other  internal  functions  as  well.  A  block  diagram  of  a  digital 
controller  capable  of  such  functions  was  presented  in  figure  8;  functions  other  than  weight 
calculations  are  discussed  below. 

Interface  and  Control  -  To  perform  all  the  functions  required  to  support  adaptation,  several 
sources  of  information  must  be  input  to  the  digital  controller.  After  manipulating  this 
information,  the  results  must  be  output  to  several  destinations.  To  accomplish  these  tasks, 
several  interface  ports  should  be  supported,  including: 

o  A  personal  computer  (PC)  interface  could  be  provided  which  might  be  used  for 
loading  algorithms,  storing  and  displaying  adaptation  data,  single  stepping  for  debug  and 
trouble  shooting,  or  enabling  full  speed  operation  with  or  without  software  breakpoints. 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 
Digital  Controller  (continued)  - 
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o  An  interface  compatible  with  a  host  radio  must  be  included.  This  allows  transfer  of 
frequency  tuning  and  system  timing  information  from  the  radio  to  the  AAP. 

o  A  bus  internal  to  the  AAP  would  allow  interfacing  to  all  the  other  AAP  subsystems. 
Tuning  of  the  synthesizers,  retreival  of  the  correlation  data,  storing  and  recalling  the  weight 
estimates,  sampling  the  IF  power  monitors,  and  controlling  step  attenuators  are  examples  of  the 
functions  provided  by  this  interface. 

o  Miscellaneous,  special  function  controls,  such  as  RF  switches  and  calibration 
components,  may  require  a  dedicated  port  with  custom  capabilities. 

o  If  properly  designed,  the  digital  controller  would  allow  for  future  expansion.  This 
could  include  more  I/O  ports,  and  digital  to  analog  and  analog  to  digital  converters. 

AAP  Functions  -  The  specific  functions  required  for  AAP  operation  are  summarized  below. 

o  Tune  the  frequency  synthesizers  per  the  information  provided  by  the  host  radio. 

o  Collect  off  time  statistical  data  sufficient  to  form  the  covariance  matrix  Ra 

o  Form  the  off  time  covariance  matrix  and  perform  a  matrix  inversion,  yielding  R0_1. 

o  Sample  IF  power  monitor  outputs.  Based  oh  these  power  levels,  control  step  attenuators 
to  ensure  the  signal  levels  into  critical  circuits  are  within  their  most  linear  range. 

o  Establish  the  timing  for  adaptation  sequences,  such  as  integration  time,  propagation 
delay  through  a  chain  of  RF  and  IF  components,  RF  or  IF  path  selection  of  switches,  time 
delay  beteen  control  signals  and  output  settling,  etc. 

o  Collect  blocks  of  on  time  statistical  data  sufficient  to  form  the  cross-correlation 
vector  R^i  several  times  within  one  timeslot. 

o  Calculate  new  estimates  of  the  weight  vector  based  on  the  equation  w(m+ 1 )  «=  R0_lRxyt 

o  Update  the  weight  estimates  applied  to  the  CPMs  as  they  are  calculated. 

o  Store  and  recall  the  weight  estimates  and  other  pertinent  data  at  the  beginning  or  end 
of  each  timeslot. 

o  Control  the  calibration  components  and  sequence  as  required. 

o  Perform  Built-In-Test  f unctions  which  provide  a  level  of  confidence  that  the  circuits 
are  operating  per  their  design  expectations. 

An  in  depth  explanation  of  how  these  digital  controller  functions  are  used  to  enable  the  AAP 
to  adapt  to  its  RF  environment  is  provided  next. 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 
Digital  Controller  ^continued!  - 
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Generically  speaking,  the  digital  controller  must  perform  different  operations  depending  on 
whether  off  time  or  on  time  processing  is  being  executed.  These  concepts  were  introduced  in 
this  section  under  paragraphs  describing  "Baseline  Choice  #1".  The  functions  performed  by  the 
digital  controller  are  described  in  terms  of  the  tasks  required  to  generate  nulls  and  a  beam  at 
the  AAP  system  output.  Although  this  explanation  is  given  with  a  particular  algorithm  in  mind, 
the  same  sequence  or  slightly  modified  sequence  could  be  used  for  a  number  of  algorithms. 

Off  time  tasks  -  At  the  beginning  of  the  timeslot,  the  frequency  control  data  for  one  advanced 
timeslot  is  input  from  the  host  radio.  The  controller  immediately  tunes  the  off  time  synthesizer 
to  this  frequency.  (The  on  time  AAP  for  this  timeslot  was  tuned  to  the  correct  frequency  at  the 
beginning  of  the  last  timeslot;  each  AAP  stays  at  the  same  frequency  for  two  timeslots).  While 
the  synthesizer  is  settling,  the  final  weights  calculated  by  this  AAP  in  the  last  timeslot  are 
stored  for  retrieval  at  the  next  visit  to  that  timeslot.  The  weights  are  then  set  to  zero  for  off 
time  operation.  Once  the  synthesizer  is  settled  and  the  input  signals  are  allowed  time  to  reach 
the  power  monitor  circuits  and  accurate  power  measurements  are  made,  the  power  monitor 
outputs  arc  sampled.  Based  on  these  measurements,  the  feedfoward  step  attenuators  are  set  to 
provide  signal  levels  into  the  correlator  circuits  which  are  within  their  linear  range.  The 
feedback  IF  oscillator  is  selected  by  controlling  the  feedback  switch.  After  these  settings  are 
settled,  the  power  monitor  levels  are  verified.  If  suitable,  the  correlators’  integrations  are 
enabled  for  a  long  integration  period  (long  relative  to  the  on  time  cross-correlation  integration 
times).  While  the  integration  is  occuring,  the  microprocessor  is  free  to  accomplish  any  overhead 
functions  such  as  monitoring  the  power  monitor  outputs  to  verify  a  consistent  RF  environment. 
After  integration,  the  correlation  results  are  sampled  by  the  digital  controller,  and  used  to  form 
the  covariance  matrix.  For  a  system  with  four  antenna  elements,  an  8  by  8  covariance  matrix 
is  constructed  (complex  calculations).  This  matrix  is  then  inverted;  the  inverted  matrix  will  be 
used  in  calculating  the  new  weight  estimates  during  the  on  time.  At  the  conclusion  of  the  off 
time  timeslot,  preparations  are  made  for  on  time  processing. 

On  time  tasks  -  The  CPM  weights  are  set  to  the  optimal  weights  as  estimated  the  last  time  this 
AAP  was  tuned  to  this  frequency  by  recalling  those  weights  from  memory.  The  input  and 
output  powers  are  then  sampled  at  the  power  monitors,  allowing  step  attenuator  adjustments 
to  be  made  to  place  the  signal  levels  into  the  correlator  in  its  linear  range.  The  feedback  switch 
is  controlled  to  select  the  IF  version  of  the  output,  thus  allowing  cross-correlation 
measurements  rather  than  baseband  conversions.  Once  the  power  levels  into  the  correlator 
circuits  are  verified,  the  integrators  are  enabled  for  a  period  of  about  128  data  bits  (assuming 
EPLRS  is  the  host  radio).  During  this  integration  period,  the  microprocessor  is  free  to  perform 
limited  overhead  functions.  At  the  end  of  the  integration  period,  the  results  arc  sampled  and 
retrieved.  These  measurements  are  used  as  the  cross-correlation  vector  for  that  block  of  data. 
This  vector,  along  with  the  inverted  covariance  matrix  from  off  time  processing,  is  used  in 
calculating  an  update  of  the  weight  estimates.  The  new  weight  estimates  are  applied  to  the 
CPMs,  causing  changes  in  the  antenna  pattern  and  output  SINR.  The  entire  on  time  sequence 
described  above  is  repeated  for  several  blocks  of  on  time  data.  After  the  last  block  of  on  time 
data,  the  weights  are  stored  for  initialization  at  the  next  visit  to  that  frequency. 
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4.3.3  ARCHITECTURE,  BASELINE  (continued) 
DC  Power  Supply  - 
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The  DC  voltages  required  by  all  the  AAP  circuits  could  be  generated  from  a  single  110  VAC 
input  within  this  DC  Power  Supply  Assembly.  The  designs  of  the  other  circuits  would  be  well 
estalished  so  an  accurate  estimate  of  the  current  require  at  each  supply  voltage  could  be 
formulated.  The  total  current  requirements,  with  considerations  for  isolating  the  supplies  of 
some  circuits  from  others,  will  help  determine  which  AC  to  DC  power  converters  should  be 
used.  Protection  features  should  include  input  overvoltage,  input  and  output  short  circuit, 
reverse  voltage  and  transient  suppression.  Converters  which  include  most  of  the  regulation  and 
protection  circuits  required  for  most  of  the  AAP  circuits  are  available. 

Some  AAP  circuits,  such  as  the  synthesizers,  the  CPMs,  and  the  correlators,  are  sensitive  to 
noise  or  transients  on  their  DC  supply  voltages  and  will  require  additional  filtering  and 
regulation.  These  additional  circuits  may  be  included  in  the  DC  Power  Supply  Assembly  or  may 
be  located  with  the  circuits  they  are  supplying.  Depending  on  the  power  dissipation  of  the 
system,  the  power  supplies  may  need  heat  sinking  and  forced  air  flow  to  maintain  reliable 
operation.  A  thermal  analysis  would  need  to  be  conducted  to  determine  those  requirements. 


Antenna  Array  - 

The  ability  of  the  AAP  to  provide  spatially  distributed  nulls  in  the  direction  of  jammers  and 
a  beam  in  the  direction  of  the  signal  of  interest  is  based  upon  the  delay  experienced  by  these 
incoming  signals  depending  on  their  direction  of  arrival  (DOA).  This  DOA  dependent  delay 
is  provided  by  the  antenna  array  via  the  physical  Separation  between  elements.  As  a  result,  the 
location  of  each  element  relative  to  one  another  (the  array  geometry)  is  critical  to  the  AAP 
performance.  To  provide  flexibility  in  the  location  of  individual  elements,  the  array  could 
allow  for  "floating"  antenna  mounts.  The  precise  location  of  each  element  could  be  adjusted  for 
a  specific  test,  allowing  empirical  determination  of  the  array  geometry.  Subsequent 
developments  could  apply  this  experimental  data  in  selecting  permanent,  fixed  locations. 

The  design  of  an  antenna  array  is  very  dependent  on  the  host  radio  the  AAP  is  operating  with. 
The  elements  of  the  array  could  be  a  version  of  the  omni-directional  whip  antenna  used  with 
the  host  radio.  Using  the  expected  antenna  gain  pattern  of  these  individual  elements,  the  array 
geometry  could  be  customized  based  upon  computer  simulated  performance  and  empirical  field 
test  data.  Using  these  constraints,  a  mechanical  engineer  would  design  a  support  structure  for 
the  antenna  elements. 
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4.3.3.7  CALIBRATION 
Introduction  - 

Circuitry  must  be  included  for  calibrating  or  correcting  three  critical  AAP  functions: 

o  Complex  Weights  -  The  weights  are  calculated  digitally,  but  implemented  at  radio 
frequencies.  Any  errors  between  the  calculated  and  the  applied  weights  will  degrade  system 
performance.  An  em  .  correction  term  may  be  generated  at  the  CPMs  which  will  optimize  the 
match  of  the  applied  and  calculated  weights. 

o  Statistical  Measurement  Circuits  -  All  weight  estimates  are  based  on  correlator 
measurements;  errors  in  these  measurements  will  result  in  inaccurate  weight  estimates  and 
degraded  system  performance.  The  correlation  circuits  could  include  self  calibration  capability, 
or  the  correction  factors  must  be  measured  and  taken  into  account  by  the  digital  controller. 

o  Amplitude,  Phase  and  Delay  Variations  -  Any  amplitude,  phase,  or  delay  variations 
in  the  RF  and  IF  chains  leading  to  the  correlator  functions  will  result  in  erroneous  correlation 
measurements  and  degraded  system  performance.  The  addition  of  an  RF  calibration  signal 
source  at  the  AAP  inputs  will  allow  these  variations  to  be  measured  in  the  correlator  circuits 
and  taken  into  account  by  the  digital  controller  for  weight  calculations. 

Each  of  these  calibration  or  correction  techniques  are  described  in  the  following  paragraphs. 

Complex  Weights  -  A  block  diagram  for  the  complete  CPM  function,  including  error  correction, 
is  shown  in  figure  14.  The  calculated  I  and  Q  weights  are  updated  by  the  digital  controller  then 
transfered  to  the  VVAs.  The  algorithm  assumes  that  the  exact  weight  calculated  is  realized  in 
the  CPM;  any  variation  f rom  this  ideal  degrades  the  jammer  rejection.  For  example,  if  the  error 
were  just  one  LSB  for  both  I  and  Q  weights,  an  eight  bit  system  would  yield  a  39  dB  null 
against  a  CW  jammer.  If  this  accuracy  were  degraded  to  2  LSBs,  the  null  would  degrade  to  just 
31  dB.  Since  the  transfer  functions  of  the  VVAs  (RF  gain  versus  control  voltage)  are  not 
perfectly  linear,  and  the  bipolar  change  in  phase  is  not  exactly  an  abrupt  180  degrees,  a 
correction  factor  is  required. 

The  sources  of  errors  in  the  VVAs  are  illustrated  as  part  of  figure  14.  They  include: 

o  A  highly  non-linear  saturation  region  in  the  V^Vj  vs.  VctJ  transfer  function  for  high 
values  of  control  voltages. 

o  Any  variations  from  the  ideal  linear  slope  in  the  V0/Vj  vs.  Veil  transfer  function 
within  the  linear  range  of  control  voltages. 

o  Any  mismatch  in  the  slopes  of  the  V0/Vj  vs.  Vct]  transfer  function  for  positive  and 
negative  control  voltages. 

o  Any  drift  in  the  null  offset  voltage;  the  null  should  occur  at  precisely  zero  volts  Vctl. 

o  Any  phase  error  resulting  from  the  insertion  phase  approaching  the  desired  values  of 
0  and  180  degrees  asymptotically.  The  ideal  values  of  precisely  0  or  180  degrees  are  never  quite 
achieved. 

o  Phase  errors  resulting  f rom  a  non-instantaneous  transition  between  0  and  180  degrees. 
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FIGURE  14.  CALIBRATION  OF  COMPLEX  WEIGHTS  (AT  RF  OR  IF) 
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4. 3.3. 7  CALIBRATION  (continued) 
Complex  Weights  (continued) 
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These  errors  can  be  corrected  for  by  adding  a  PROM  to  each  control  port  of  CPM  which  maps 
the  desired  (calculated)  weight  to  the  required  control  word  to  achieve  that  weight.  To  allow 
for  temperature  variations  at  the  VVAs,  a  sensor  can  be  integrated  into  the  CPM  assembly 
which  generates  a  digital  word  representing  its  operating  temperature.  This  word  can  be 
combined  with  the  calculated  control  word  to  address  the  correction  PROM.  The  entire  CPM 
circuit  could  be  characterized  over  temperature  by  an  automated  test  program  which  would  test 
selected  voltages  and  temperatures,  then  interpolate  to  get  the  data  points  which  fall  between 
measured  values. 

SlMlitkai  Measurements  -  The  calibration  circuits  are  included  in  the  block  diagram  for  the 
statistical  measurement  circuit  shown  in  figure  15.  The  basic  circuit  is  used  to  either  measure 
the  correlation  between  two  IF  signals  or  to  provide  quadrature  downconversion  of  an  IF  signal 
to  baseband.  The  calibration  technique  must  therefore  be  capable  of  measuring  and  correcting 

errors  in  either  of  these  modes.  The  circuit  shown  accomplishes  this;  its  operation  is  explained 
below. 

To  calibrate  the  correlation  of  two  IF  signals,  the  feedback  switch  selects  the  termination 
thereby  setting  one  of  the  two  IF  signals  to  zero.  Regardless  of  the  signal  present  at  the 
feedforward  port  of  the  multiplier,  its  output  should  be  zero.  To  assure  operation  in  the 
multiplier’s  linear  range,  however,  a  test  signal  should  be  provided  and  applied  at  the 
feedforward  port.  The  analog  switches’  positions  will  be  selected  to  provide  integration  of  the 
lowpass  filter  output,  as  is  the  case  when  IF  correlation  is  desired.  The  integrated  output  is  then 
converted  to  a  digital  word.  During  normal  operation  this  is  the  correlator  output-  for 
calibration  it  is  a  test  signal  which  indicates  the  amount  of  error  accumulated  in  all  the 
associated  circuitry.  This  error  is  input  to  a  PROM,  which  determines  the  polarity  and  amount 

°LCOIreCtu°n  Wh'Ch  mus.t  bc  apphed  t0  cancel  the  detected  error.  This  correction  factor  is 
added  to  the  total  corrections  accumulated  by  previous  samples.  The  composite  correction  term 
is  converted  to  an  analog  voltage  which  reduces  the  total  correlation  error  by  summing  with 
the  output  of  the  lowpass  filter  which  follows  the  multiplier. 


To  calibrate  the  baseband  measurement  accuracy,  the  feedback  switch  selects  the  IF  oscillator. 
A  quadrature  power  sphttei ,  not  shown  in  f  igure  1 5,  provides  both  cos(w„.t)  and  sin(w11Pt)  terms 
for  the  quadrature  downconversion.  In  this  case,  the  output  of  the  lowpass  filter  is  either  the 
in-phase  or  quadrature  baseband  representation  of  the  signal  present  at  the  feedforward  IF 
port  of  the  multiplier.  When  the  feedfoward  signal  is  disabled  (using  an  IF  switch  or  step 
attenuator  not  shown  in  figure  15),  the  baseband  signal  should  be  zero.  If  not  the  same 
calibration  process  discussed  above  is  used  to  adjust  the  calibration  signal  until  the  composite 
output  is  forced  to  zero.  In  this  case,  the  analog  switches  bypass  the  integration  function  just 
as  they  do  for  baseband  data  collection. 

The  techniques  described  above  are  simple,  easy  to  implement,  and  inexpensive.  Even  so,  an 
alternative  solution  is  available:  the  calibration  measurements  could  be  sampled  by  the  digital 
controller  and  the  errors  could  be  accounted  for  in  its  subsequent  calculations 
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FIGURE  15.  CALIBRATION  OF  STATISTICAL  MEASUREMENTS 
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4.3.3.7  CALIBRATION  (continued) 

Amplitude.  Phase  and  Delay  Variations  -  The  RF  and  IF  paths  leading  to  the  statistical 
measurement  circuits  must  be  precisely  calibrated  otherwise  they  may  introduce  errors  in  the 
statistical  measurements.  The  correlation  of  input  signals  relative  to  one  another  is  measured 
for  construction  of  the  off  time  covariance  matrix.  The  correlation  of  each  input  signal  to  the 
output  signal  are  measured  for  construction  of  the  on  time  cross-correlation  vector.  These 
correlation  measurements  assume  they  are  made  using  data  directly  from  the  antenna  elements 
and  directly  from  the  output.  Any  amplitude,  phase  or  delay  errors  between  the  inputs  and  the 
statistical  measurement  circuits,  or  the  output  and  the  statistical  measurement  circuits,  will 
result  in  erroneous  correlation  data.  As  an  example,  the  measurements  and  resulting  weight 
estimates  are  based  on  the  relative  separation  between  each  of  the  antenna  elements.  If  a  phase 
shift  or  delay  is  introduced  in  one  feedforward  IF  chain  leading  to  the  correlators  which  is  not 
introduced  in  all  other  feedforward  IF  chains,  the  data  for  that  one  channel  will  be  in  error. 
Each  channel  is  expected  to  have  different  and  unpredictable  amplitude,  phase  and  delay  terms 
in  their  overall  transf er  functions.  These  must  be  calibrated  to  provide  the  required  correlation 
data  accuracy  necessary  to  support  null  and  beam  steering. 

Since  the  data  required  to  construct  the  correlation  matrix  is  collected  as  laseband  data  while 
the  cross-correlation  vector  data  is  collected  as  IF  correlation  measurements,  two  separate 
calibrations  are  required  f or  these  two  procedures.  The  measurements  required  to  perform  these 
calibrations  are  available;  there  are  a  number  of  sequences  which  could  provide  the  necessary 
calibration  accuracy.  Example  calibration  procedures  for  baseband  data  collection  and  for  IF 
correlation  measurements  are  presented.  These  procedures  may  not  be  optimal  (due  to  Phase  I 
SBIR  time  and  budget  constraints),  but  are  representative  of  the  types  of  measurements  and 
procedures  which  are  available  with  the  proposed  architecture. 

Calibrating  baseband  data  measurements  -  Figure  16  shows  a  simplified  block  diagram  of  the 
AAP  in  which  the  RF  and  IF  paths  which  are  being  calibrated  arc  highlighted.  An  RF 
generator  has  added  which  is  coupled  into  each  of  the  antenna  elements’  signal  paths.  All  the 
CAL  signals  into  each  channel’s  front  end  are  assumed  to  be  amplitude  and  phase  matched  to 
one  another.  The  components  in  the  paths  from  the  antenna  input  to  the  feedfoward  port  of 
the  statistical  measurement  circuits  introduce  amplitude  and  phase/delay  variations  between 
channels.  One  procedure  for  measuring  these  variations,  determining  correction  factors,  and 
applying  these  corrections  to  measured  data  is  outlined  below. 

1.  Enable  the  IF  Oscillator  in  the  feedback  path,  tune  the  input  RF  Cal  Generator  to  the 
desired  frequency  and  enable  its  output,  and  tune  the  LO  for  proper  downconversion.  All  of 
these  signals  are  derived  from  the  same  source  and  are  therefore  coherent. 

2.  Use  the  statistical  measurement  circuits  to  measure  the  phases  and  amplitudes  of  each 
feedfoward  signal  relative  to  the  reference  signal  from  the  feedback  path.  The  outputs  of  the 
quadrature  downconverter  provides  this  information  in  rectangular  coordinates.  The  digital 
controller  can  sample  this  data  and  convert  it  to  polar  coordinates.  For  example,  the  data  for 
channel  k  would  have  the  form  RkA0k. 

3.  Choose  one  channel  to  be  the  reference  channel;  all  other  channels  will  be  calibrated 
to  match  the  amplitude  and  phase  of  that  reference  channel.  Assuming  channel  one  is  the 
reference  channel,  the  correction  factors  are  ClkA0lk  *  R^Rfc  _ek). 
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FIGURE  16.  CALIBRATION  OF  BASEBAND  DATA  COLLECTION 
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4.  Measure  the  actual  data  (this  occurs  after  the  cal  sequence  determines  the  correction 
factors).  Process  this  data  to  form  the  elements  of  the  covariance  matrix  Rq,  which  are  the 
expected  values  of  (xkX|).  Note  that  the  correction  factors  associated  with  each  component  of 
the  expected  values  are  deterministic,  and  therefore  can  be  moved  outside  the  expected  value 
operation.  This  allows  the  correction  factors  to  be  applied  to  the  elements  of  the  covariance 
matrix  after  integration  rather  than  directly  to  each  sample  of  baseband  data.  The  number  of 
calculations  (and  therefore  time  delay)  is  greatly  reduced  for  the  formation  of  Rq.  The 
covariance  ma'rix  elements  are  in  rectangular  coordinate  format,  so  they  must  be  converted 
to  polar  coordinates.  Assume  this  data  is  MklA.7kl. 

*5.  For  each  element  of  the  covariance  matrix,  the  correction  factor  is  provided  by  the 
calibration  sequence:  CklA.0kl  =  (ClkA.6lk)(C1IA6u)  =  ClkCuA.(lk  +  1().  This  multiplication  factor 
is  applied  to  the  measured  data  MklA.7kl  to  provide  the  calibrated  covariance  data  to  be  used 
by  the  algorithm  for  adaptation  calculations.  This  data  is  used  to  construct  a  calibrated  R,,. 

Calibrating  IF  correlation  measurements  -  This  procedure  involves  calibrating  two  sets  of  paths 
and  applying  the  calibration  factors  in  different  ways  to  achieve  the  desi-cd  results.  The  goal 
of  this  IF  correlation  calibration  procedure  is  to  calibrate  out  amplitude  and  phase/delay  errors 
in  the  cross-correlation  data  associated  with  on  time  processing.  This  correlation  data  is  used 
in  vector  form  (designated  R^)  by  the  adaptation  algorithm;  the  elements  of  this  vector  are 
E[xky*]  where  E[]  denotes  an  expected  value.  Calibrating  these  measurements  can  be 
accomplished  by  a  two  step  process:  Step  1  calibrates  all  the  direct  signal  paths  from  the  input 
to  the  output.  As  explaned  below,  this  effectively  provides  a  calibrated  output  term  y  *  Ewkxk 
for  the  adaptation  algorithm.  Step  2  calibrates  the  inputs  (xk)  to  one  another;  when  steps  1  and 
2  are  combined,  each  product  xky*  is  calibrated. 

The  first  step  in  this  process  calibrates  the  paths  shown  in  figure  17;  it  recognizes  that  the 
theoretical  output  of  the  A  AP  is  a  summation  of  all  the  weighted  inputs  (y  «  Ewkxk)  and  that 
one  calibration  of  th"  ^  iput  only  is  not  sufficient  unless  the  paths  from  input  (x)  to  output 
(y)  are  identical.  Since  this  entire  procedure  is  necessary  because  those  paths  are  not  identical, 
a  means  of  calibrating  each  wkxk  term  is  required.  The  first  step  of  this  two  step  calibration 
is  designed  to  apply  a  different  correction  factor  to  each  complex  weight  (the  wk  terms). 
Calibration  data  is  accumulated  by  measuring  the  amplitude  and  phase/dclay  variations 
between  the  direct  signal  paths  from  the  inputs  to  the  output.  After  the  power  combiner  which 
sums  the  weighted  RF  signals  from  each  channel,  the  path  to  the  feedback  port  of  the 
statistical  measurement  circuits  is  common  to  all  channels.  Therefore,  measurements  made  at 
th:  IF  correlators  indicate  the  variations  between  the  RF  paths  only.  A  reference  signal  for 
measuring  each  channel’s  amplitude  and  phase  is  provided  by  the  feedforward  signal  from  one 
of  the  channels.  The  channel  associated  with  the  correlator  circuit  selected  is  the  reference 
channel  for  these  measurements;  as  done  above,  a  correction  factor  will  be  determined  for  all 
other  channels  to  result  in  a  set  of  data  which  is  f ully  calibrated  with  one  another.  An  example 
procedure  for  accomplishing  this  calibration  is  listed  below. 

1.  Tur.e  and  enable  the  RF  calibration  generator  and  the  LO.  Select  the  FB  IF  signal 
rather  than  the  IF  oscillator  signal  using  the  feedback  switch.  The  RF  calibration  signal  present 
in  all  channels  will  be  downconver'ed  and  applied  to  the  statistical  measureme..t  circuits’  IF 
feedforward  ports,  but  only  ''hannel  one’s  circuits  will  be  used  in  this  step  of  the  calibration 
procedure.  Therefore,  the  FF  IF  signal  for  channel  one  is  the  reference  to  which  the  feedback 
signals  will  be  calibrated. 
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2.  Enable  the  real  component  of  the  complex  weight  for  channel  one  only.  To  accomplish 
this,  set  Wj  =  1  +  jO,  all  other  wk  =  0  +  jO. 

3.  Measure  the  amplitude  and  phase/delay  of  the  feedback  signal  (from  channel  one 
only)  relative  to  the  reference  FF  IF  signal  using  channel  one’s  IF  correlator.  The  outputs  of 
the  correlator  circuits  are  in  rectangular  coordinates;  convert  this  data  to  polar  coordinates. 
This  calibration  data  is  RjA-Oj. 

4.  Repeat  step  numbers  2  and  3  for  each  of  the  channels,  yielding  all  the  cal  data  RkA.0k. 

5.  Once  the  summation  occurs,  the  only  way  each  term  in  the  summation  can  be 
corrected  for  (without  added  circuitry)  is  by  applying  the  correction  factors  to  the  weights 
rather  than  to  the  data.  Since  channel  one  is  selected  as  a  reference,  a  correction  factor  similar 
to  the  one  calculated  in  the  baseband  data  calibration  procedure  is  required  to  calibrate  the 
other  channels  to  channel  one’s  amplitude  and  phase  reference.  As  in  the  baseband  case,  these 
calibration  factors  may  be  written  as  ClkA.6lk  =  Rj/Rk  A.(0j  -  0k). 

6.  The  weight  estimates  are  calculated  by  the  adaptation  algorithm;  the  correction 
factors  described  above  multiply  the  weight  estimates  before  application  to  the  RF  weighting 
elements.  This  assures  that  y  *=  Ewkxk,  without  distributed  amplitude  and  phase/delay  errors 
corrupting  each  term  of  the  summation. 

The  second  step  of  the  IF  correlation  calibration  involves  calibrating  the  feedfoward  paths  in 
all  the  channel’s  relative  to  one  another.  These  paths  are  highlighted  in  figure  18.  This 
effectively  provides  calibration  of  all  the  inputs  xk.  This  procedure  is  almost  identical  to  the 
baseband  data  collection  calibration  procedure  except:  1)  The  feedback  signal  used  as  a 
reference  is  not  the  FB  oscillator,  but  is  the  signal  which  originates  in  the  RF  cal  generator  and 
is  allowed  to  pass  through  one  complex  weight.  2)  Some  of  the  statistical  measurement  circuits 
used  for  IF  correlation  measurements  are  different  than  those  used  for  baseband  data 
generation  and  manipulation.  3)  Only  one  term  needs  corrected  (xk)  rather  than  two  terms  of 
a  product  (xkx,).  Following  the  same  procedure  as  outlined  in  those  paragraphs:  calibration 
measurements  are  made,  correction  factors  to  calibrate  channels  2  thru  4  to  channel  1  are 
calculated,  real  data  measurements  are  made,  then  the  correction  factors  are  applied  to  the 
measured  data.  Again,  since  the  corrections  are  deterministic,  they  can  be  applied  after 
integration. 

Combining  the  two  steps  described  above  completes  the  IF  correlation  calibration  process.  The 
correction  factors  for  the  weights  are  determined  per  step  1.  The  input  data  correction  factors 
are  determined  per  step  2.  The  calibration  data  collection  cycle  is  then  complete,  external  signal 
reception  is  enabled,  and  adaptation  processing  begins.  The  IF  cross-correlation  measurements 
are  made,  block  by  block.  The  correction  factor  for  the  input  path  variations  are  applied  to 
these  cross-correlation  measurements.  The  adaptation  algorithm  calculates  weight  estimates 
based  on  this  corrected  data,  then  the  correction  factors  from  step  1  (which  corrects  for  input 
to  output  variations)  are  applied  to  the  weights  to  determine  the  desired  weights  to  be  applied 
to  the  CPMs. 
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FIGURE  18.  CALIBRATION  OF  IF  CORRELATION  DATA,  STEP  2 


5  fc| 

&  A 


o  o 


o  o  o 


Advanced 

CounterMeasure 

Systems 


84 


4. 3. 3.8  COMPLETE  BASELINE  AAP  RECEIVE  ARCHITECTURE 

The  baseline  AAP  receive  architecture  has  been  developed  through  out  section  4.3.3.  Options 
for  critical  components  and  subsystems  have  been  presented,  and  choices  for  the  baseline  were 
selected  and  justified.  The  resulting  AAP  receive  architecture  is  shown  as  a  fairly  detailed 
block  diagram  in  figure  19. 

Although  all  the  features  of  this  architecture  were  discussed  in  detail  in  previous  sections,  the 
followirg  comments  are  provided  as  concluding  remarks  which  summarize  the  AAP  receive 
architecture. 

o  The  antenna  array,  shown  in  the  block  diagram  to  include  four  elements,  receives  the 
RF  signals  present  and  provides  the  time  delay  between  sensors  which  is  exploited  by  the  AAP 
to  form  nulls  in  the  direction  of  jammers  and  a  beam  in  the  direction  of  the  desired  signal. 

o  Several  components  are  included  at  the  inputs  of  the  AAP  which  are  used  exclusively 
for  internal  calibration.  These  include  the  input  switches,  the  CAL  Generator,  the  input  four 
way  power  splitter,  and  the  input  couplers. 

o  A  front  end  network  is  provided  for  each  channel  which  is  shared  by  the  odd  and 
even  processors.  These  networks  determine  the  RF  operating  band  and  establish  the  AAP  unit 
noise  figure. 

o  For  anticipative  processing,  most  of  the  adaptive  processing  components  must  be 
repeated.  Since  four  channels  are  included  in  the  baseline,  eight  identical  processor  modules 
are  required.  These  modules  include  downconversions  of  the  input  signal,  the  IF  chains  which 
condition  the  feedforward  signals  for  correlation  measurements,  power  monitors  of  the  input 
signals,  statistical  measurement  circuits  (IF  correlator  /  baseband  downconverter  with 
associated  analog  or  digital  processing),  and  the  complex  weighting  devices. 

o  The  LO  signals  must  frequency  hop  with  the  host  radio.  Since  the  odd  processor  is 
operating  in  advance  of  the  even  processor,  two  sets  of  LO  signals  are  required,  hence  two 
synthesizer  modules  are  shown. 

o  The  feedfoward  paths  include  filtering  to  limit  correlation  estimates  to  the  message 
bandwidth,  thereby  maximizing  the  SIR.  Also  included  are  a  delay  adjustment,  step  attenuator 
for  amplitude  adjustment,  and  gain  to  overcome  losses  in  other  components. 

o  The  statistical  measurement  circuits  shown  include  the  calibration  circuits  required 
to  assure  the  accuracy  of  the  IF  correlation  and  baseband  data. 

o  The  complex  weighting  circuits  include  the  error  correction  circuits  which  map  the 
calculated  weight  values  to  the  transfer  functions  of  the  CPMs  used  to  implement  them. 
Temperature  dependent  variations  are  also  accounted  for  and  correction  is  provided. 

o  A  feedback  path  is  provided  which  is  very  similar  to  the  feedforward  path.  This  chain 
includes  downconversion,  filtering  to  the  message  bandwidth,  a  delay  adjustment,  a  step 
attenuator  for  amplitude  adjustment,  and  a  power  monitor.  The  feedback  path  also  includes 
an  IF  switch  which  determines  the  input  to  the  statistical  measurement  circuits.  One  of  three 
inputs  may  be  selected:  1)  The  feedback  signal  from  the  output,  used  for  cross-correlation 
vector  measurements.  2)  The  IF  oscillator,  which  is  used  in  the  collection  of  off  time  baseband 
data.  3)  A  50  ohm  termination,  used  in  calibrating  the  statistical  measurement  circuits. 
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4. 3.3.8  COMPLETE  BASELINE  AAP  RECEIVE  ARCHITECTURE  (continued) 

o  The  even  and  odd  processors  include  identical  circuits  and  work  in  parallel.  The 
output  is  formed  by  alternatively  selecting  the  processor  whose  output  includes  the  signal  of 
interest  (the  processor  working  with  on  time  data).  This  selection  is  done  by  the  RF  output 
switch. 


oThc  digital  controller  interfaces  with  external  devices  (such  as  the  host  radio), samples 
data  from  internal  subassemblies,  controls  the  internal  subassemblies,  and  performs  the 
calculations  required  for  adaptation.  One  digital  controller  can  serve  both  the  even  and  odd 
processors. 


4. 3. 3. 9  TRANSMIT  ARCHITECTURE 

The  concept  of  using  the  information  acquired  via  receive  adaptation  to  direct  the  transmit 
energy  of  the  host  radio  has  been  demonstrated  in  this  and  other  studies.  The  simulated  system 
performance  improvements  are  encouraging  and  deserve  additional  investigation.  While  the 
simulations  (which  assume  ideal  components)  are  encouraging,  the  directive  and  retrodirective 
transmission  techniques  will  be  difficult  to  implement  with  real  components.  A  top  level  block 
diagram  of  a  directive/retrodirective  AAP  which  illustrates  some  of  the  issues  requiring 
further  investigation  is  shown  in  figure  20.  Examples  of  these  issues  are  listed  below. 

o  Even  though  the  host  radio  offers  a  range  of  transmit  power  levels,  the  lowest  level 
is  too  high  -  the  transmit  weighting  elements  will  have  limited  dynamic  range.  The  transmit 
power  from  the  host  radio  needs  to  be  decreased  either  within  the  host  radio  (undesirable  since 
this  requires  modifying  an  established  design  or  customizing  some  units  for  AAP  operation) 
or  within  the  AAP.  This  may  cause  thermal  problems  in  the  AAP  or  degradation  in  the  quality 
of  the  transmitted  signal. 

o  A  calibration  procedure  which  maps  the  transmit  weights  calculated  by  the  digital 
controller  to  the  actual  weights  provided  by  the  RF  devices  used.  This  problem  is  similar  to  the 
receive  AAP’s;  the  same  type  of  correction  circuits  associated  with  the  CPMs  may  be  used  for 
the  transmit  path  as  well  as  the  receive  path. 

o  A  method  of  calibrating  the  transmit  weights  to  the  receive  weights  is  required. 

o  The  host  radio’s  transmit  signal  is  power  divided  into  four  signals;  each  of  these 
signals  are  processed  in  parallel.  The  output  processing  includes  the  power  divider,  an 
attenuator  to  decrease  the  levels  into  the  weighting  devices,  the  weighting  devices  themselves, 
RF  amplifiers  and  filters  (probably,  although  not  shown)  and  the  high  power  amplifier.  All  of 
these  transmit  paths  must  be  calibrated  relative  to  one  another.  This  is  similar  to  the  amplitude 
and  phase/delay  calibration  procedures  identified  in  the  receive  AAP  discussions. 

o  One  high  power  amplifier  is  required  for  each  channel  of  the  AAP.  Although  the 
transmit  power  required  of  each  amplifier  will  be  less  than  the  power  available  from  the  host 
amplifier  (since  the  AAP  outputs  are  summed  together),  the  inefficiencies  of  the  amps  will 
result  in  considerably  higher  DC  power  dissipations  and  may  cause  thermal  problems. 

An  effort  was  made  by  ACM  Systems  to  address  these  issues  and  develop  an  architecture 
capable  of  supporting  directive  and  retrodirecUve  transmissions,  but  time  and  money 
limitations  required  a  premature  conclusion  of  this  effort. 
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4.3.4  COMPUTER  MODELING  AND  SIMULATIONS 

Note:  The  notation  for  matrices  and  vectors  in  the  following  section  differs  slightly  from  that 
of  the  other  sections  in  this  report.  For  instance  M  represents  a  matrix,  and  v  represents  a 
vector. 

4.3.4. 1  INTRODUCTION 

Extensive  simulations  were  performed  to  demonstrate  the  effectiveness  of  blind  adaptive 
antenna  array  algorithms  and  architectures  on  frequency  hopped  signals  in  a  corrupt 
environment.  This  signal  environment  consists  of  directional  jammers,  and  a  directional, 
frequency  hopped  (FH)  signal  of  interest  (SOI).  The  data  is  modeled  as  the  output  of  an 
adaptive  array  processor  (AAP)  front  end  with  receiver  noise,  and  with  adjustable  gain,  phase, 
and  delay  variations  between  individual  elements.  The  vector  signals  generated  are  essentially 
down  converted,  de-hopped,  and  sampled. 

Two  algorithms,  Dominant  Mode  Maximin  (Block  Update)  and  Maximin  [20,  21,  22],  were 
implemented.  Each  algorithm  exploits  the  FH  characteristic  of  the  SOI  to  reject  the  jammers 
and  beam  up  on  the  SOI.  The  Dominant  Mode  is  a  Maximin-type  algorithm  which  takes  an 
eigen-equation  approach,  whereas  Maximin  is  a  gradient  method.  This  study  examines  the 
application  of  the  Dominant  Mode  algorithm  to  Army  tactical  radio  systems  for  the  first  time. 

In  order  to  evaluate  the  simulations,  a  variety  of  outputs  were  produced.  These  include  weight 
vectors,  aperture  vector  (steering  vector)  estimates,  and  output  signal  to  interference  and  noise 
ratios  (SINR).  The  weight  vectors  were  used  to  generate  receive  antenna  patterns  which  give 
a  visual  measure  of  the  effectiveness  of  the  adaptive  array  processor.  In  a  similar  fashion,  the 
aperture  vector  estimates  were  used  to  generate  transmit  patterns  and  power  footprints. 

4. 3. 4. 2  SIMULATION  DESCRIPTION 
Environment  and  AAP  Front  End  - 

The  data  used  in  the  simulations  was  generated  with  the  EPLRS  system  in  mind.  Due  to  the 
rapid  convergence  of  the  Dominant  Mode  algorithm,  the  amount  of  data  needed  to  perform  the 
simulations  could  be  reduced.  We  chose  256  chips  per  hop  for  processing,  whereas  the  EPLRS 
system  has  about  2400  chips  per  hop.  This  is  equivalent  to  saying  that  adaptation  was 
discontinued  after  256  chip  periods.  We  have  also  ignored  for  now  the  guard  time  which 
contains  no  SOI.  That  is,  processing  occurred  inside  a  window  during  which  an  SOI  should  be 
present  (rather  than  including  the  period  during  which  it  is  known  no  SOI  will  be  present). 

The  hop  increment  of  the  EPLRS  system  is  3  MHz;  the  data  is  spread  over  frequency  using 
direct  sequence  spread  spectrum  (DSSS)  techniques  with  a  3  Mbit/sec  chip  rate.  The  EPLRS 
signal  is  frequency  hopped  over  eight  channels,  or  bins.  In  the  simulations,  the  center 
frequencies  of  these  bins  are  (426  +  3n)  MHz,  where  n=0,l,2,...,7  is  the  bin  number  (see  Table 
1).  The  27  MHz  RF  bandwidth  represents  a  percent  bandwidth  of  6.2  percent. 


Table  1  Bin  center  frequencies  and  their  corresponding  bin  numbers. 


Bin  frequency 

426 

429 

432 

435 

438 

441 

444 

447 

Bin  number 

0 

1 

2 

3 

4 

5 

6 

7 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 

For  generality,  the  SOI  generated  for  each  simulation  was  a  pseudo-random  (PN)  constant 
modulus  (CM)  sequence  (s(t)=e'-*®^)  with  uniformly  distributed  phase  (6(t)=u[-ir,ir)).  A  timing 
error  (SOI  offset  advance)  was  provided  as  an  option,  but  not  investigated  as  part  of  this  study. 

The  jammer  signals  were  modeled  optionally  as  either  PN  CM  noise  (vCM(t)=e''ie^^)  with 
uniformly  distributed  phase,  or  complex  white  gaussian  (CWG)  noise 
(vCWG(t)=CWGN(p=O,0=l))  with  zero  mean  and  unit  variance.  The  receiver  noise  was  strictly 
CWG  (u(t)=CWGN(|i=0,o=l)). 

It  should  be  noted  that  these  sequences  are  independent  of  each  other,  independent  from  hop 
to  hop,  and  independent  between  on  and  off  times.  The  receiver  noise  is  independent  from 
element  to  element.  And,  fresh  sequences  were  generated  for  each  simulation. 

Two  sets  of  data  were  generated  for  each  hop  (timeslot):  off  time  data  and  on  time  data.  The 
off  time  data  contains  jammers  and  receiver  noise  only,  as  where  the  on  time  data  contains 
jammers,  receiver  noise,  and  the  SOI.  In  these  simulations,  the  off  time  data  was  modeled  as 
though  it  was  collected  from  the  frequency  bin  of  the  on  time  data  one  timeslot  prior  to  the 
actual  hop  to  that  frequency  bin.  That  is,  anticipative  off  time  data  collection  is  being 
simulated.  Thus,  the  off  time  data  will  have  the  same  noise  and  directional  jammers  present 
as  the  on  time  data.  Remember,  though,  the  jammers  and  noise  are  independent  between  on  and 
off  times,  but  have  the  same  statistics. 

The  data  provided  to  the  algorithm  is  essentially  down  converted,  de-hopped,  and  sampled.  It 
is  modeled  as  vector  output  of  an  antenna  array  front  end.  This  modeling  is  accomplished  by 
creating  aperture  vectors  for  the  SOI  and  for  each  jammer  present  in  the  current  f  requency  bin. 
The  aperture  vectors  provide  for  the  option  of  gain,  phase,  and  delay  variations  between 
antenna  elements.  In  this  study,  however,  time  did  not  permit  investigation  of  these  variations. 
The  dynamic  range  of  the  system  is  also  taken  into  account  when  adding  the  receiver  noise  and 
applying  the  automatic  gain  control  (AGC). 

Table  2  lists  the  SOI  parameters;  Table  3  lists  the  various  selectable  system  parameters  with 
their  default  values;  and  Table  4  lists  the  jammer  parameters. 


TABLE  2  SOI  Parameters 


Parameter 

Description 

Defaults 

o. 

SOI  direction  of  arrival  (DOA) 

0° 

SOI  timing  advance  (  0  s  t,  s  Nh  ) 

0 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 


TABLE  3  System  Parameters 


Parameter 

Description 

Defaults 

M 

number  of  sensors 

4 

(SU 

sensor  locations  (units  of  center 
frequency  wavelength,  Xc) 

At  vertices  of  a 
unit  square. 

DR 

dynamic  range  (dB) 

40 

voltage  gain,  delay,  and  phase 
offset,  respectively,  of  the  mth 
sensor’s  RF  chain,  assumed  to  be 
independent  of  frequency.  Delay  is 
normalized  to  the  period  of  the 
center  frequency. 

{1,0,0} 

fc 

center  of  frequency  band  (MHz) 

436.5 

Mh 

number  of  hop  bins 

8 

Nh 

number  of  chips  per  hop 

256 

TABLE  4  Jammer  Parameters 


Parameter 

Description 

Defaults 

lsnoi 

total  number  of  jammers 

2 

{6.} 

jammer  DOA’s 

O 

o 

O 

o 

{^k|} 

jammer  frequency  range  (bin  #’s) 

40*  -  {  3,  7  } 

110*  -{0,4  } 

{CM,} 

constant  modulus  (CM)  flag 

T  ->  CM  jammer 

F  «>  CWG  jammer 

40*  -F 

110*  -  F 

{S  JR,} 

in-bin  signal  to  jammer  ratios  (dB) 

40*  -  -20 

110*  -  -20 
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4.3. 4.2  SIMULATION  DESCRIPTION  (continued) 

The  following  comprehensively  describes  the  generation  of  the  simulated  antenna  array  front 
end  output.  First,  the  antenna  array  apertures  and  attenuators  are  initialized  from  the 
parameters  specified  in  the  tables  above.  The  jammer  levels  are  incorporated  into  their 
respective  aperture  vectors  for  simplicity.  The  apertures  and  attenuators  are  then  rescaled  by 
taking  into  account  the  dynamic  range  of  the  receiver  array  and  applying  the  AGC. 

(In  the  following,  k  =  frequency  bin  number,  m  =  sensor  number.) 


1.  Set  initial  (unnormalized)  receive  antenna  gains: 


RCV  /<2jt()T 

“fan  =  8me 


2.  Set  initial  SOI  receive  apertures: 


SOI 


“fan 


-/2«(7)d>e,> 

e 


a(6,)  » 


'  cosOj' 

,s“e., 


3.  Set  SNOI  (signal  not  of  interest  =>  jammer)  variables: 

Lk  *  No.  SNOI’s  at  frequency  bin  k 
0k!  -  SNOI  DOA’s  of  jammer  1  at  bin  k 
CMkl  =  SNOI  CM/CWG  flag 

4.  Set  initial  SNOI  attenuations  and  receive  apertures: 

af 0/  =  10-m‘,2° 


SNOI 

“ton 


SNOI 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 

5.  Set  the  rescaling  gain  factors  (AGC  normalization): 


6.  Rescale  (normalize)  all  apertures  and  receiver  gains: 


RCV 

.  RCV 

RCV 

“ton 

-  “  m 

X 

“ton 

SOI 

RCV 

SOI 

“ton 

-  “ton 

X 

“ ton 

SNOI 

RCV 

SNOI 

“tom 

-  “ton 

X 

“tom 

7.  Calculate  maximum-attainable  SINR  for  each  frequency  bin: 

SINR(k)  = 


where  the  sensor  subscript,  m,  has  been  incorporated  into  the  vector  notation  of  the 
apertures.  And,  the  composite  interference  and  noise  is  £(n). 


Once  the  apertures  and  attenuators  are  computed,  the  various  signals  are  generated  and 
combined  accordingly: 

Over  the  i‘h  hop: 

1.  Determine  k,  the  ith  hop  frequency  bin  (pseudo  random  number  generator). 

2.  Set  the  composite  interference  field,  g(n),  for  n=-Nh  to  Nh-1  (n=-Nh  to  -1  is  the 
anticipatory  portion)  by  combining  jammers,  Vj(n),  and  receiver  noise,  j£(n): 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 

h 

a(n)  =  u(H)*  S<f>) 

/-I 


Uin)  ~  LID.  CWGN(QJ) 


v/n)  ~ 


c  f 

zf.n)  , 


6~I.I.D.  u[-n,n),  CM{  =  T 
Z,~IJ.D.  CWGN( 0,1),  CM,  =  F 


where  I.I.D.  -  Independent,  Identically  Distributed. 


3.  Add  SOI: 


m  « 


fl^0/s(n+ 1 ,)+«(/!)  , 


-zsinz-xs+Nh-l 


{  fl(«)  ,  -Nhzn<- xt  ,  -x  g+Nh-l<nzNh-l 

5(n)  ~  ,  0"/J.D.  w[-n,n) 

(i(n)  is  referred  to  as  the  augmented  hop  data). 


Algorithms  - 

Processing  can  be  done  with  the  off  time  data  (no  SOI  present)  as  well  as  the  on  time  data.  Two 
algorithms  were  incorporated  into  the  processor  model,  but  adequate  simulations  of  their 
effectiveness  was  not  possible  due  to  time  limitations  on  the  SBIR.  The  two  algorithms  are: 
steepest  descent  null  steering  and  inverse  power  method  null  steering.  Both  techniques 
minimize  the  output  power  in  the  absence  of  the  SOI;  only  SNOIs  are  present  -  consequently, 
minimizing  output  power  during  the  off  time  steers  nulls  without  regard  for  SOI  DOAs.  The 
inverse  power  method  null  steering  algorithm  accomplishes  this  by  adapting  the  weights  to  the 
minimum  eigenvector  of  the  off  time  auto-correlation  matrix.  These  preprocessing  techniques 
could  be  explored  in  subsequent  studies  since  they  should  decrease  initial  adaptation  times. 

The  weights  are  initialized  isotropic  at  the  start  of  simulation  (prior  to  off-time  processing). 

The  two  on  time  algorithms  explored  in  this  study  are  Dominant  Mode  Block  Update  and 
Maximin.  The  computer  simulations  of  the  Dominant  Mode  Block  Update  algorithm  is 
described  in  the  following  text.  The  simulation  of  the  Maximin  algorithm  was  attempted,  but 
the  results  obtained  were  inconsistent;  to  avoid  wasting  time  on  established  results,  this  effort 
was  abandoned.  The  reader  is  referred  to  (20,  21,  22]  for  descriptions  of  the  Maximin 
implementation  and  its  simulations. 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 

Dominant  Mode  Block  Update  - 

The  flow  chart  in  Figuie  21(a-c)  depicts  in  detail  the  implementation  of  the  Dominant  Mode 
Block  Update  algorithm.  Table  5  lists  the  selectable  parameters  used  in  the  flow  chart.  Table 
6  gives  a  brief  description  of  the  variables  appearing  in  the  flow  chart. 

The  Dominant  Mode  Block  Update  algorithm  relies  on  the  power  method  to  adapt  the  weights 
to  the  maximum  eigenvector  as  shown  in  the  flew  chart.  In  the  implementation  of  the 
algorithm,  however,  the  on  time  autocor-elation  matrices  are  not  exactly  the  same  from  c'oek 
to  block.  This  method  is  still  valid,  thoug..,  since  the  matrices  will  exhibit  ine  same  statistics. 
As  the  collect  times  approach  infinity,  the  block  to  block  variation  is  eliminated.  This  is 
actually  a  ’stochastic  power  method’  being  implemented.  The  optimal  set  ot  weights  from  block 
to  block  (after  convergence)  will,  for  the  most  part,  be  unchanging  except  for  an  overall  gain 
and  phase  factor.  The  adapted  weights  are  therefore  saved,  according  to  the  current  frequency 
bin,  from  hop  to  hop  in  order  to  significantly  reduce  adaptation  time  in  subsequent  visits  to 
ine  same  frequency  bin. 

There  are  a  couple  of  architecture  dep-  ident  variations  which  were  investigated  as  part  of  this 
study.  The  first  uses  R.xyi  instead  of  RjW  in  the  adaptat  on.  This  was  simulated  oy 
bypassing  the  smoothing  applied  to  from  block  to  block.  That  is,  since  RXY1  =  £ki(Ci)w  » 
then  E-xyi  =  &i—  by  setting  R,i  •  Ri(m)  each  hop.  This  variation  is  motivated  by  t'  :  f-’ct 
that  computation  and  stor-ge  of  RXyi  much  less  intensive  than  computation  and  storage  ot 
Rt .  This  variation  did  not  significantly  alter  the  performance  of  the  processor,  and  has 
consequently  been  considered  a  viable  hardware  reducing  option. 

The  second  variation  attempted  a  similar  operation  to  construct  the  off  time  statistics,  Eo 
The  simulation  consisted  of  building  ^  one  column  (vector)  at  a  tim~  by  processing  the  off 
time  data  in  blocks.  Notice,  the  first  column  of  (denoted  [R_n]n  )  is  equal  to  F--Y0 
w.«=£0  =  [l  0  0  ...  0]T  ,  since  Rxyo  =  ■  ln  general,  then,  [R^  *=  Ro£j  •  The  variation 

considered  consists  of  obtaining  the  first  column  from  the  first  off  time  data  block  by  setting 
w  «  ;  the  second  column  from  the  second  block  w  =  £j  ;  and  so  on.  Thus,  R*,  would  not  have 

to  be  computed  directly,  reducing  hardware  components  (and  thus,  simplifying  the 
architecture). 

This  approach  was  not  successful,  though.  The  resulting  matrix  had  lost  some  of  its  most 
important  properties:  hermitian  and  positive  definite.  An  attempt  to  restore  these  propeit>',s 
(specifically  the  hermitian  property)  did  not  redeem  this  approach.  These  results  c^n  oe 
explained,  though,  by  realizing  that  Rq  is  not  the  same  from  block  to  block.  The  statistics 
f rom  blo'k  to  block  are  similar,  but,  only  if  a  large  enough  number  of  samples  were  used  would 
the  matrices  begin  to  be  equal.  The  relation  between  this  convergence  and  the  number  of 
samples  used  per  block  is  approximately  (l/v'n).  Thus,  the  various  columns  will  not  combine 
to  give  a  functional  estimate  of  the  off  time  statistics,  Rq  due  to  our  finite  collect  interval. 
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4. 3. 4. 2  SIMULATION  DESCRIPTION  (continued) 


TABLE  5  User  Specified  Processing  Parameters 


Parameter 

Description 

Defaults 

N0 

off  time  collect  interval 

256 

N, 

on  time  block  collect  interval 

32 

factor  for  accelerated  convergence 

0 

TABLE  6  Algorithm  Variables 


Variable 

Description 

k 

the  frequency  bin  of  the  current  hop  /  timeslot 

Mx 

number  of  on  time  blocks  to  process  (  Nh  /  Nj ) 

x(n) 

the  augmented  hop  data  generated  above 

Xo(n) 

off  time  data 

2Li(n) 

on  time  data 

Ro 

estimated  off  time  auto  correlation  matrix 

Ri(m) 

estimated  on  time  auto  correlation  matrix  over  block  m 

S.1 

smoothed  E.i(i)  over  i*0,...,m 

stored  weight  vector  for  kth  frequency  bin 

Ok 

stored  aperture  estimates  for  kth  bin 

Yk 

stored  SINR  estimates  for  kth  bin 

y(n) 

output 

hk 

number  of  visits  to  bin  k 
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Figure  21(a)  Dominant  Mode  Block  Update  Implementation  Flow  Chart 
For  an  arbitrary  hop  /  timeslot: 
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4.3.4. 3  SIMULATION  RESULTS 

In  order  to  evaluate  the  simulations,  a  variety  of  outputs  were  produced.  These  include  weight 
vectors,  aperture  vector  (steering  vector)  estimates,  and  output  signal  to  interference  plus  noise 
ratios  (SINR).  The  weight  vectors  were  used  to  generate  receive  antenna  patterns  which  give 
a  visual  measure  of  the  effectiveness  of  the  adaptive  array  processor.  The  equation  used  to 
generate  these  patterns  is: 


*V«>  = 


M'fjf'Ce) 


where  0g  is  the  direction  of  arrival  cf  the  SOI.  In  a  similar  fashion,  the  aperture  vector 
estimates  were  used  to  generate  transmit  patterns  and  radiated  power  footprints.  Recall  the 
calculation  of  the  maximum  attainable  SINR: 

SINR^k)  =  taf')s£8‘  .  af 


—4  4 


The  output  (true)  SINR  referred  to  in  the  following  discussion  is  estimated  as: 


SINRom(k)  = 


|<of 


i/  )2 


The  results  of  several  simulated  jammer  and  array  scenarios  follows.  All  of  the  simulations 
used  the  Dominant  Mode  Block  Update  algorithm. 

Receive  Adaptation  Baseline  - 

The  baseline  test  cases  simulate  a  7  element  antenna  array,  with  one  element  in  the  center  and 
the  remaining  equally  spaced  around  a  circle  of  diameter  X.CV2  (see  Figure  28,  7  elements  - 
centered).  The  scenario  for  the  5  jammers  present  in  the  environment  is  shown  in  Figure  22. 
Parameters  not  specified  are  assigned  the  defaults  listed  in  the  tables  in  section  4.3.4.2. 

Figure  23  shows  the  SINR  convergence  within  individual  time  slots  for  the  first  12  hops  of  a 
pseudo  random  hop  sequence.  The  sequence  of  frequency  bins  appears  at  the  bottom  of  the 
figure.  The  corresponding  SINR  data  (for  first  visits  only)  appears  in  Table  7.  One  can  see  the 
algorithm  adapts  within  the  first  3  data  blocks  for  each  initial  visit  to  a  bin.  The  jagged  line 
running  across  the  top  is  the  maximum  attainable  SINR.  The  maximum  attainable  SINR  varies 
slightly  from  bin  to  bin  due  to  the  jammer  scenario  of  the  various  bins.  Upon  revisiting  a  bin. 
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Figure  22  Baseline  Jammer  Scenario 


Jammer  Modulations: 


Constant  Modulus  White  Noise 


Complex  White  Gaussian  Noise 
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Table  7 


Output  8 I NR  for  initial  visit  to  each  frequency  bin. 


Bin  0 

Blk 

SINR 

(dB) 

Bin  l 


Bin  2 


Bin  3 


30.4 

16.5 

17.2 

17.3 
17.3 
17.3 
17.3 
17.3 


SINK 

(dB) 


30.8 

12.7 

15.1 

15.1 

15.1 

15.1 

15.1 

15.1 


SINR 

(dB) 


30.8 

11.4 

15.3 

15.3 

15.3 

15.3 

15.3 

15.3 


SINR 

(dB) 


30.8 

14.8 

15.8 
15.8 
15.8 
15.8 

15.8 

15.9 


Bin  4 

Blk 

SINR 

(dB) 

30.8 

10.8 

15.9 

15.9 

15.9 

15.9 

15.9 

15.9 


Bin  5 

Blk 

SINR 

(dB) 

31.1 

11.0 

15.3 

15.3 

15.4 

15.4 

15.4 

15.4 


Bin  6 

- 1 

Blk 

SINR 

(dB) 

-30.8 

9.5 

15.4 

15.4 

15.4 

15.4 

15.4 

15.4 


15.5 


Bin  7 

Blk 

SINR 

(dB) 

-30.4 

15.6 

17.1 

17.2 

17.2 

17.2 

17.2 

17.2 


17.3 


Table  8 


Output  SINR  for  various  jammer  levels. 


SJR 

=  -10 

Blk 

i 

SINR 

(dB) 

.0 
.5 
.5 
37.5 
37.5 
37.5 
37.5 
37.5 


SJR 

=  -20 

Blk 

SINR 

(dB) 

.0 
.2 
.5 
27.5 
27.5 
27.5 
27.5 
27.5 


SJR 

=  -30 

Blk 

SINR 

(dB) 

-30.0 

13.9 

17.5 

17.5 

17.5 

17.5 

17.5 

17.5 


SJR 

-  -40 

Blk 

SINR 

(dB) 

-40.0 

-0.1 

6.9 

7.5 

7.5 

7.5 

7.5 

7.5 


37.5 


27.6 


17.6 


4.3.4. 3  SIMULATION  RESULTS  (continued) 
Receive  Adaptation  Baseline  (continued)  - 


103 


the  weights  are  already  adapted  to  give  the  maximum  attainable  SINR  (refer  to  time  slots  8, 
10,  1 1  and  12).  This  is  so  because  the  weights  are  stored  and  recalled  for  each  frequency  bin. 
The  convergence  to  the  optimal  weight  vector  (ie.,  the  maximizing  of  the  SINR)  is  extremely 
rapid  —  within  96  chips  (32  chips  per  block,  3  blocks). 

The  antenna  patterns  for  the  adapted  weight  vectors  from  selected  frequency  bins  appear  in 
Figure  24.  The  grid  lines  are  spaced  lOdB  apart.  Notice,  in  each  case  the  nulls  are  quite  sharp 
and  are  at  least  50dB  below  the  beam  in  the  direction  of  the  SOI. 

In  order  to  simulate  a  communications  network  such  as  the  EPLRS  TDMA,  simulations  were 
run  where  the  SOI  moved  on  subsequent  visits  to  various  frequency  bins.  Figure  25  shows  the 
antenna  patterns  for  the  first  and  second  visits  to  frequency  bins  1  and  3,  with  the  SOI  moving 
in  each  case  to  one  of  the  inadvertent  nulls  of  the  first  visit.  In  each  case,  the  weights  adapted 
to  produce  a  beam  in  the  direction  of  the  SOI,  while  maintaining  the  nulls  in  the  direction  of 
the  jammers.  The  weights  re-converged  in  less  than  3  blocks  in  all  cases  tested. 

Another  situation  encountered  in  network  communications  is  the  absence  of  an  SOI. 
Investigation  of  the  Dominant  Mode  algorithm  reveals  that  in  the  absence  of  an  SOI,  the  off 
time  and  on  time  statistics  should  be  the  same  (R0  -  R,).  This  implies  that  the  weights  should 
not  change  since  w(k+l)  *  R^Rjwfk)  -  Iw(k)  «  w(k).  Thus,  if  the  weights  were  already 
adapted  to  the  optimal  for  a  particular  SOI,  and  then  the  SOI  disappeared  in  a  subsequent  visit 
to  the  bin,  the  weights  should  not  change.  In  reality,  however,  the  weights  will  drift  because 
the  off  time  and  on  time  matrices  are  not  identical  due  to  the  finite  collect  time.  Figure  26 
shows  the  drift  in  the  antenna  patterns  during  the  absence  of  an  SOI  in  5  consecutive  visits  to 
frequency  bin  3.  me  SINR"  shown  in  each  case  is  the  value  obtained  from  the  equation  for 
SINROUT  above  with  the  same  SOI  aperture  vector  as  if  an  SOI  were  present.  As  shown  in 
Figure  27,  this  "SINR"  also  drifts  from  the  maximum  attainable  SINR  over  time. In  actual 
operation,  this  drift  could  be  stifled  by  the  processor  once  the  loss  of  the  SOI  is  detected.  A 
way  of  detecting  this  is  by  a  significant  drop  in  the  value  of  the  dominant  mode  eigenvalue 
v  hich  is  calculated  by  the  Dominant  Mode  algorithm. 

Architecture  Dependent  Simulations  - 

All  of  the  simulations  model  the  off  time  data  as  one  long  block  of  anticipative  data  integrated 
over  the  entire  timeslot.  The  Dominant  Mode  Block  Update  algorithm  utilizes  the  power  method 
with  matrices  based  on  smaller  blocks  of  on  time  data  which  change  from  block  to  block  but 
are  statistically  similar.  This  results  in  a  stochastic  power  method  implementation.  All  of  the 
simulations  demonstrate  the  algorithm  to  be  robust,  supporting  the  validity  of  this  application. 
Notice,  as  the  collect  time  (within  a  block)  gets  larger  (more  chips  collected)  the  on  time 
statistical  matrices  become  more  constant  from  block  to  block.  The  Dominant  Mode  algorithm 
increases  the  collect  time  in  subsequent  blocks  within  a  timeslot  by  weighted  averaging  of  the 
previous  blocks.  This  is  the  smoothing  of  the  on  time  auto  correlation  matrix,  ^  ,  in  the 
algorithm  flow  chart  (Figure  21). 

The  architecture  actually  utilizes  Eocyi  instead  of  RjW  for  adaptation.  This  was  simulated 
by  suppressing  the  smoothing  mentioned  above.  This  subtle  variation  did  not  significantly 
impact  the  algorithm’s  performance.  Thus,  the  architecture  only  needs  to  compute  the  cross 
correlation  between  the  input  and  output  for  the  on  time  statistics. 
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Receive  Antenna 


A.  441  MHZ  (bin  5) 
Jl:  30*,  CWG,  30dB 
J3 :  180*,  CWG,  20dB 
J4 :  235*,  CWG,  20dB 
J5 :  270*,  CM,  20dB 


C. 
Jl: 
J4 : 
J5 : 


444  MHz  (bin  6) 
30*,  CWG,  30dB 
235*,  CWG,  20dB 


270* 


CM,  20dB 


9 


for  Baseline  Simulation 


J2 : 

110* , 

CM, 

3  0dB 

J3 : 

180*  , 

CWG, 

2  0dB 

J4 : 

235*  , 

CWG, 

20dB 

D.  435  MHz  (bin  3) 


Jl: 

30  *  , 

CWG, 

30dB 

J3 : 

180*  , 

CWo  , 

20dB 

J4 : 

235*  , 

CWG, 

20dB 
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Figure  25  Receive  Antenna  Patterns  for  Changing  SOI  DOA 


(SOI  DOA  =0*) 

429  MHZ 
J2 :  110 
J3 :  180 
J4:  235 


(SOI  DOA  =36*) 

(bin  1) 

,  CM,  30dB 
,  CWG,  20dB 
,  CWG,  20dB 


C.  First 

visit 
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4. 3. 4. 3  SIMULATION  RESULTS  (continued) 
Architecture  Dependent  Simulations  - 
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This  technique  is  not  directly  applicable  to  the  off  time  statistics,  since  the  entire  off  time  auto 
correlation  matrix  is  needed  for  the  Dominant  Mode  algorithm.  A  variation  of  this  method  was 
attempted,  though.  The  simulation  consisted  of  building  Rg  one  column  at  a  time  by 
processing  the  off  time  data  in  blocks.  Notice,  the  first  column  of  *o  (denoted  [R^g )  is  equal 
to  R.xyo  if  w  =  £<,  =  [1  0  0  ...  0]T  ,  since  RyVn  =  R^w  .  In  general,  then,  [Rgjj  =  R^j  .  The 
variation  considered  consists  of  obtaining  the  first  column  from  the  first  off  time  data  block 
by  setting  w  =  ^  ;  the  second  column  from  the  second  block  w  =  ^  ;  and  so  on.  Thus,  Rg 
would  not  have  to  be  computed  directly,  reducing  hardware  components  (and  thus,  simplif ying 
the  architecture).  This  approach  proved  invalid  over  the  finite  collect  time. 

In  trying  to  achieve  success  with  this  innovative  approach,  we  noticed  the  resulting  matrix  had 
lost  some  of  its  most  important  properties:  hermitian  and  positive  definite.  An  attempt  to 
restore  these  properties  (specifically  the  hermitian  property)  did  not  redeem  this  approach. 
These  results  can  be  explained  by  realizing  that  Jig  is  not  the  same  from  block  to  block.  The 
statistics  from  block  to  block  are  similar,  but,  only  if  a  large  enough  number  of  samples  are 
used  would  the  matrices  begin  to  be  constant.  The  relation  between  this  convergence  and  the 
number  of  samples  used  per  block  is  approximately  (1/v^n).  Thus,  with  a  finite  number  of 
samples  collected  the  various  columns  will  not  combine  to  give  a  f unctional  estimate  of  the  off 
time  statistics,  Rg  . 

An  important  result  is  the  overall  effectiveness  of  storing  and  recalling  weights  from  hop  to 
hop.  The  storage  table  consists  of  a  weight  vector  for  each  frequency  bin.  This  provides  for 
immediate  SOI  recovery  in  subsequent  visits  to  a  particular  frequency  bin,  as  long  as  the  RF 
environment  has  not  changed  significantly.  Additionally,  storing  and  recalling  weights  for 
separate  frequency  bins  provides  very  effective  frequency  compensation.  This  result  was 
illustrated  in  Figure  23,  where  the  SINR  of  a  revisited  bin  begins  with  nearly  the  maximum 
attainable  value  because  of  convergence  during  the  previous  visit. 

Antenna  Geometry  and  Spacing  - 

The  effects  of  the  antenna  array  geometry  and  spacing  were  investigated  in  this  simulation. 
This  includes:  varying  number  of  elements;  varying  geometry  of  arrays;  and  varying  the 
spacing  of  the  elements.  We  examined  the  performance  with  3,  4,  5,  and  7  elements.  For  the 
4  and  7  element  arrays,  the  effect  of  changing  the  geometry  was  studied.  Specifically,  one 
geometry  consisted  of  having  one  clement  in  the  center,  with  all  remaining  elements  on  the 
circumference  of  a  circle;  the  other  placed  all  elements  on  the  circumference.  In  this 
simulation,  all  of  the  array  diameters  are  y/2kc.  See  Figure  28  for  the  array  geometries.  The 
effect  of  varying  the  spacing  of  the  elements  was  performed  for  the  4  element  array,  with  all 
of  the  elements  on  the  circumference  of  a  circle.  The  diameters  compared  arc:  A.  ,  \/2 A.  , 
2AC ,  and  2>/2Ac  . 

The  jammer  scenario  for  each  of  these  simulations  is  the  same.  Each  simulation  was  forced  to 
take  place  in  frequency  bin  3.  Only  one  timeslot  was  visited  (8  blocks).  There  were  two 
jammers  present  in  frequency  bin  3  for  these  simulations:  40”  and  1 10*,  SJR  -  -20dB  for  both. 
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Figure  28 


Antenna  Array  Geometries 
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Antenna  Geometry  and  Spacing  (continued!  - 

Figure  29  shows  the  antenna  patterns  for  the  various  number  of  elements  and  geometries. 
Notice  the  improvement  of  the  antenna  patterns  as  the  number  of  elements  increases.  There 
is  not  a  significant  improvement,  though,  between  the  non-centered  4  element  case  and  the 
centered  5  element  case,  nor  between  the  3  element  case  and  the  centered  4  element  case.  The 
non-centered  4  element  case  outperformed  the  centered  4  element  case  in  this  scenario.  The 
centered  7  element  case,  on  the  other  hand,  appears  to  be  a  more  favorable  pattern  than  the 
non-centercd  7  element  case.  These  results  are  highly  dependent  on  angles  of  arrival  of  the  SOI 
and  SNOIs,  however,  and  are  not  conclusive.  The  model  exists  for  future  exploration. 


Figure  30  shows  the  antenna  patterns  for  the  various  array  diameters.  The  most  favorable 
pattern  appears  to  be  the  one  corresponding  to  a  diameter  of  Slk  .  Notice,  as  expected,  the 
number  of  inadvertent  nulls  increases  as  the  diameter  increases. 


Varying  Jammer  Levels  - 

The  effect  of  varying  the  jammer  levels  was  also  studied.  The  antenna  geometry  chosen  for 
WaS  the  4  element<  non-centered  array,  with  diameter  y/2k.  .  Jammer  levels 
(SJR)  of  -IOdB,  -20dB,  -30dB,  and  -40dB  were  studied  for  a  single  jammer  present  in  frequency 
civro  Wltu  3  direction  of  arrival  of  40°.  Only  one  timeslot  was  processed  (8  blocks).  Output 
cIktd  3t  the,start  °f  each  b,ock  is  listcd  for  cach  jamroer  level  in  Table  8.  Maximum  attainable 
SINR  is  included  at  the  end  of  each  list  for  comparison.  A  plot  showing  the  convergence  time 
ot  the  output  SINR’s,  normalized  to  the  corresponding  maximum  attainable  SINR  in  each  case, 
appears  in  Figure  31.  Notice  the  weaker  jammers  are  nulled  out  slightly  more  quickly  than  the 
higher  level  jammers.  In  each  case,  though,  the  weights  still  converge  within  3  block  updates. 
Actually,  the  weights  adapted  within  2  block  updates  for  the  weaker  jammer  levels  --  that  is 
all  but  SJR  *  -40dB  .  * 

One  test  performed  on  the  SJR  -  -40dB  case  provided  quicker  adaptation  -  within  2  blocks, 
this  simulation  set  the  factor  for  accelerated  convergence,  p,  equal  to  1.0. 
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ire  29  RECEIVE  ANTENNA  PATTERNS  vs.  NUMBER  OF  ELEMENTS 

(DIAMETER  OF  ANTENNA  ARRAYS  -  A.  *  J2) 
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Figure  30  RECEIVE  ANTENNA  PATTERNS  VS.  ARRAY  DIAMETER 

(4  ELEMENT  ANTENNA  ARRAY) 


DIAMETER  =  XQ 
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Table  8  output  SINR  for  various  jammer  levels. 
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4. 3. 4. 3  SIMULATION  RESULTS  (continued) 

Directive  Communications.  Antenna  Patterns  - 

Retrodirective  and  directive  antenna  patterns  were  generated  using  the  weights  and  aperture 
estimates.  For  the  4  element  array  which  is  non-centered  with  diameter  '/2kz,  the  receive 
antenna  pattern,  retrodirective  antenna  pattern,  and  retrodirective  radiated  power  footprint 
a>-e  shown  in  Figure  32.  In  this  figure,  the  antenna  pattern  grids  are  lOdB  increments.  The 
power  footprints  have  a  linear  scale,  and  are  normalized  to  the  transmit  direction  of  0e  (angle 
of  intended  receiver). 

The  antenna  patterns  and  radiated  power  footprints  in  Figure  33  are  for  directive 
transmissions  using  different  amounts  of  smoothing  applied  to  the  apertures.  Here  smoothing 
means  each  aperture  estimate  in  subsequent  visit  to  a  particular  frequency  bin  is  incorporated 
into  an  on-going  average  of  all  prior  aperture  estimates  in  that  bin  (see  Dominant  Mode  Block 
Update  algorithm  flow  chart).  The  patterns  clearly  improve  with  more  smoothing  (more 
timeslots  used  to  average).  The  ideal  patterns  were  generated  using  the  known  apertures 
calculated  in  the  Simulation  Description,  section  4.3.4.2,  and  are  included  for  comparison. 

For  each  of  the  above  cases,  two  jammers  were  present  (simulation  is  of  frequency  bin  3).  The 
arrival  directions  are  40°  and  110°,  and  the  SJR  is  -20dB  for  both. 

Directive  transmit  patterns  should  improve  when  more  elements  are  used.  Simulation  results 
shown  in  Figure  34  demonstrate  this  for  4,  7,  and  19  elements.  These  patterns  were  generated 
using  the  known  apertures.  See  Figure  28  for  the  geometry  of  the  19  element  array.  The 
intended  receiver  is  assumed  to  be  at  0°  in  all  of  the  transmit  cases  analyzed. 


Advanced 

CounterMeasure 

Systems 


Figure  32  4  ELEMENT  ANTENNA  RECEIVE  AND  RETRODIRECTIVE  TRANSMIT  PATTERNS 
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TRANSMIT  ANTENNA  PATTERN  TRANSMIT  POWER  FOOTPRINT 


Figure  33  4  ELEMENT  ANTENNA  ARRAY 

(SMOOTHING  APERTURE  ESTIMATES) 


19  ELEMENTS 
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4. 3. 4. 3  SIMULATION  RESULTS  (continued) 

Directive  Communications.  LPI  - 

The  application  of  directive  communications  to  an  LPI  system  is  described  in  sections  4. 3. 2. 7 
and  4. 3. 2.8.  The  results  of  the  corresponding  performance  simulations  are  presented  here. 

Figure  35  shows  a  set  of  intercept  footprints  for  four  antenna  arrays.  Each  set  includes  three 
cases:  1)  both  the  receive  and  transmit  arrays  having  isotropic  antenna  patterns  (this  is  the 
baseline  situation,  provided  for  reference),  2)  the  receive  array  having  an  adapted  pattern  but 
the  transmit  array  having  an  isotropic  pattern,  and  3)  the  receive  array  having  an  adapted 
pattern  and  the  transmit  array  providing  directive  transmissions.  The  dramatic  reduction  in 
intercept  footprint  areas  is  evident  in  the  shadings  which  represent  these  three  situations. 

The  first  example  is  a  four  clement  array  with  the  sensors  located  on  the  circumference  of  a 
circle.  In  this  case,  the  transmitter  -  receiver  link  is  positioned  relative  to  the  array  geometry 
such  that  virtually  the  largest  possible  intercept  footprint  results. 

The  second  example  shows  the  same  array  rotated  45  degrees.  The  resulting  intercept  footprint 
area  drops  significantly. 

The  third  example  also  simulates  a  four  element  array,  except  in  this  case  one  of  the  sensors 
is  located  in  the  center  of  the  array. 

The  fourth  example  is  a  seven  element  array.  The  expected  reduction  in  intercept  footprint 
area  as  the  number  of  elements  is  increased  is  demonstrated  by  comparison  of  this  example  to 
the  examples  of  the  four  element  arrays. 

The  dependence  of  this  area  on  the  number  of  elements  is  further  illustrated  in  figure  36, 
which  shows  the  footprint  area  reduction  versus  the  receiver  to  transmitter  link  direction  for 
some  of  the  examples. 

Directive  Communications.  A/J  - 

The  application  of  directive  communications  to  an  A/J  system  is  described  in  sections  4. 3.2.7 
and  4. 3. 2. 8.  The  results  of  the  corresponding  performance  simulations  arc  presented  here. 

Figure  37  shows  the  A/J  jammer  footprints  of  the  same  arrays  discussed  above.  The  same 
situations  are  represented:  1)  both  the  receive  and  transmit  arrays  having  isotropic  antenna 
patterns  (this  is  the  baseline  situation,  provided  for  reference),  2)  the  receive  array  having  an 
adapted  pattern  but  the  transmit  array  having  an  isotropic  pattern,  and  3)  the  receive  array 
having  an  adapted  pattern  and  the  transmit  array  providing  directive  transmissions.  Again,  the 
shadings  illustrate  the  dramatic  reduction  in  jammer  footprint  areas  for  these  three  situations. 

The  first  plot  uses  the  four  element  array,  all  sensors  on  the  circumference,  with  the  array 
rotated  45  degrees  f rom  the  original  case.  The  second  plot  uses  the  four  element  array  with  one 
sensor  at  its  center.  The  third  plot  represents  the  seven  element  array. 

The  areas  are  greatly  reduced,  and  are  difficult  to  see  without  changing  scale.  A  close-up  view 
of  the  critical  areas  for  each  of  the  previous  three  plots  are  shown  in  the  fourth  plot.  Again, 
the  dependence  upon  the  number  of  elements  in  the  array  is  demonstrated. 
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4.3.5  DOCUMENTATION 

This  Phase  I  SBIR  effort  has  been  well  documented  since  the  program’s  inception.  Monthly 
reports,  formally  delivered  to  the  Contracting  Officer  and  informally  to  other  technical 
representatives  from  the  Government,  have  provided  regular  status  reporting  of  technical 
accomplishments  and  cost  /  schedule  issues.  Five  monthly  reports  have  been  provided,  due  the 
20th  of  each  month  from  September  1990  to  January  1991.  The  last  deliverable  requirement  of 
this  SBIR  program  is  this  final  report. 

In  addition  to  these  reports,  a  proposal  for  a  follow-on  Phase  II  effort  was  submitted.  This 
proposal  is  entitled  "A  Multi-Purpose  Adaptive  Array  Processor  Hardware  Development  for 
Frequency  Hopped,  Direct  Sequence  Spread  Spectrum  Radios"  and  is  dated  January  7,  1991. 

4.4  ESTIMATES  OF  TECHNICAL  FEASIBILITY  /  PERFORMANCE  EVALUATIONS 

4.4.1  INTRODUCTION 

This  study  has  addressed  the  application  of  adaptive  array  processing  to  ground  based  Army 
tactical  radios.  It  was  organized  into  three  main  tasks:  algorithm  development,  architecture 
development,  and  computer  modeling  and  performance  simulation.  The  technical  feasibility 
of  the  algorithm  and  architecture  are  discussed,  and  the  results  of  the  performance  simulations 
are  evaluated  below. 

4.4.2  ALGORITHM  TECHNICAL  FEASIBILITY 

The  algorithm  developed  under  this  study  is  the  Dominant  Mode  Block  Update  algorithm.  It’s 
feasibility  and  application  to  FHSS  and  DSSS  Army  radios  has  been  proven  by  derivation,  a 
convergence  and  stability  analysis,  and  performance  analysis,  and  has  been  demonstrated  by 
computer  simulation. 

The  algorithm  developed  provides  flexibility  for  customization  to  a  specific  communication 
network.  The  method  for  data  collection  is  not  entirely  dictated  by  the  algorithm;  for  example, 
the  period  or  band  that  the  off  time  data  is  acquired  from  can  vary  depending  on  the 
application.  The  equation  which  generates  updates  of  the  weight  estimates  can  be  implemented 
using  on  time  data  in  vector  (R,,,)  or  matrix  (Rt)  form.  Monitoring  critical  points  within  the 
AAP  and  calculation  of  key  parameters  will  allow  modification  to  the  adaptive  process 
depending  on  the  jamming  scenario.  For  example,  if  on  time  data  suggests  that  no  SOI  is 
present,  the  algorithm  may  cease  adaptation  and  hold  the  weights  at  their  last  valid  state.  The 
signal  levels  internal  to  the  AAP  could  be  monitored  and  adjusted  for  optimal  adaptation 
performance.  The  algorithm  could  adjust  the  variables  which  affect  convergence  time 
depending  on  the  data  collected,  such  as  changing  the  value  of  u  in  the  power  method  matrix 
A  *  R0’Iri  *  uI-  If  a  specific  Army  radio  is  identified  for  operation  with  an  AAP  employing 
this  algorithm,  many  operational  features  of  the  algorithm  could  be  tailored  to  achieve  specific 
goals. 

The  application  of  AAP  to  directive  and  retrodirective  communications  was  addressed  as  well. 
Figures  of  merit  for  a  fully  adapted  (receive  and  transmit)  AAP  system  were  identified  for  LPI 
and  A/J  applications:  intercept  f ootprint  and  maximum  receiver  input  power,  respectively.  The 
analysis  shows  these  figures  of  merit  to  be  greatly  improved  over  an  isotropic  communication 
system  and  to  nearly  approximate  the  optimal  performance. 
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4.4.3  ARCHITECTURE  TECHNICAL  FEASIBILITY 

The  feasibility  of  the  architecture  developed  in  this  study  is  heavily  supported  in  the 
architecture  section  of  this  report.  That  section  includes  examples  of  how  every  element  of  the 
AAP  could  be  implemented  in  hardware;  in  many  cases  multiple  implementations  are  available 
and  the  merits  of  those  options  are  discussed.  Examples  of  AAP  elements  end  methods  of 
implementation  are  given  below: 

o  The  complex  weights  could  be  implemented  at  RF,  IF  or  baseband.  If  realized  at 
baseband,  digital  or  analog  circuits  could  be  used.  The  weights  could  be  set  using  magnitude 
and  phase  adjustments  or  in-phase  and  quadrature  adjustments.  The  suggested  implementation 
is  at  RF  (or  IF,  depending  on  application),  with  in-phase  and  quadrature  adjustments. 

o  The  statistical  measurement  circuits  could  be  based  on  RF,  IF  or  baseband 
multiplications;  the  sampled  baseband  data  could  be  fowarded  to  a  digital  calculator  which 
determines  the  signal  statistics  or  the  integrated  statistical  data  could  be  sampled  by  the  digital 
calculator.  The  architecture  suggested  allows  for  both  IF  multiplications  (for  cross-correlation 
vector  measurements  directly)  and  quadrature  downconversions  (for  generating  baseband  data 
used  in  performing  covariance  matrix  calculations).  This  provides  all  the  data  needed  for 
employing  a  number  of  adaptation  algorithms  with  very  little  added  complexity. 

o  It  is  suggested  that  the  weight  calculations  be  done  digitally;  this  is  really  the  only 
valid  option  if  the  architecture  is  to  support  more  than  one  radio  or  adaptation  algorithm. 
Furthermore,  a  microprocessor  based  circuit  is  suggested,  which  would  provide  additional 
flexibility  for  the  entire  AAP  system. 

o  All  the  major  subassemblies  are  identified  and  discussed  as  well,  but  they  are  not 
unique  to  AAP  applications.  For  example,  the  frequency  synthesizer  must  be  fast  hopping  -  a 
multi-loop  phase-locked  loop  synthesizer  is  suggested.  The  RF  front  ends,  downconverter 
circuits,  RF  and  IF  chains,  dc  power  supplies  and  the  antenna  array  are  also  discussed,  and 
descriptions  of  commercially  available  components  for  implementing  them  are  provided. 

o  The  key  components,  items  such  as  SAW  filters,  termination  insensitive  mixers, 
quadrature  multipliers,  etc.  are  suggested  for  use  and  an  explanation  of  their  importance  is 
given. 

The  performance  of  this  architecture  depends  on  the  accuracies  of  several  of  its  circuits:  the 
statistical  measurement  circuits,  the  weight  calculator,  the  complex  weighting  circuits,  and  the 
IF  and  RF  amplitude  and  phase  variations  through  various  paths.  These  accuracies  are  greatly 
improved  by  calibration;  components  and  procedures  for  implementing  these  calibrations  are 
described  in  detail,  and  help  assure  the  AAP  system  hardware  feasibility. 

A  major  advantage  of  this  architecture  over  others  is  that  it  allows  implementation  of  many 
algorithms,  a  flexibility  that  could  be  exploited  in  further  development  efforts. 


4.4.4  SIMULATED  PERFORMANCE  EVALUATION 

This  section  provides  an  evaluation  of  the  computer  simulation  results  presented  in  section 
4.3.4.  All  the  following  comments  refer  to  simulated  results,  and  underscore  the  significance 
of  the  results  in  the  context  of  an  operating  communication  system. 
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4.4.4  SIMULATED  PERFORMANCE  EVALUATION  (continued) 

o  The  fundamental  operation  of  the  AAP  was  verified  -  the  output  of  the  AAP 
approaches  the  maximum  attainable  SINR  in  all  cases  tested. 

o  The  AAP  converges  very  quickly  in  all  cases  tested.  Equating  the  results  to  the  EPLRS 
format,  convergence  is  achieved  in  less  than  256  chips,  or  85  usee  (one  chip  is  about  333  nsec). 
Although  a  timeslot  lasts  about  2  msec,  the  transmission  occurs  in  a  burst  of  about  800  usee. 
Thus,  the  maximum  attainable  SINR  is  achieved  in  about  ten  percent  of  the  transmission  time. 

o  With  this  fast  convergence  speed,  the  AAP  is  ideal  for  frequency  hopped  applications. 
The  simulations  provided  consistent  results  in  a  severe  jamming  environment,  with  a  different 
mix  of  high  level  jammers  at  each  frequency  bin. 

o  The  improved  performance  expected  by  storing  and  recalling  weights  unique  to  each 
frequency  bin  was  verified.  Once  the  bin  is  visited  and  the  weights  converge,  those  weights 
provide  a  very  good  starting  point  for  adaptation  during  the  next  visit  to  that  bin. 

o  The  antenna  patterns  illustrate  the  depth  of  the  nulls  in  the  directions  of  jammers  and 
the  high  gain  maintained  in  the  SOI  direction  (a  "beam"). 

o  The  AAP  supports  network  communications  as  well  as  LPI  communications.  The 
direction  of  the  SOI  is  not  consistent  in  the  network  case,  causing  potential  problems  for 
algorithms  attempting  to  maximize  SINR.  The  worst  case  condition  is  to  have  the  SOI  direction 
during  the  present  timeslot  be  the  same  as  an  inadvertant  null  from  the  previous  timeslot.  The 
impact:  the  SOI  is  partially  rejected  (30  dB  or  more)  at  the  beginning  of  adaptation.  This  case 
was  simulated,  and  the  AAP  achieved  convergence  at  rates  consistent  with  the  results  above. 

o  A  SOI  is  expected  to  be  present  during  on  time  processing  by  the  algorithm.  If 
conditions  are  such  that  the  level  of  the  SOI  is  too  low  for  detection,  the  previously  optimized 
weights  will  slowly  drift  away  from  optimal.  Although  the  amount  of  drift  was  demonstrated 
to  be  acceptable,  the  algorithm  could  easily  be  modified  to  cease  adaptation  under  these 
conditions. 

o  The  anticipative  method  of  collecting  off  time  data  was  proven  effective. 

o  The  use  of  cross-correlation  vector  measurements  R^  rather  than  auto-correlation 
matrix  measurements  R,^  for  on  time  data  collection  was  proven  valid.  This  will  speed 
adaptation  in  future  hardware  developments. 

o  Some  antenna  array  investigation  was  conducted  but  this  eff ort  was  not  sufficient  in 
scope  to  state  any  verifiable  conclusions. 

o  The  AAP  achieved  the  maximum  attainable  SINR  with  consistently  fast  convergence 
times  for  multiple  jammers  and  a  range  of  SJR  from  -10  to  -40  dB. 

o  The  use  of  adapted  receive  weights  to  calculate  the  transmit  weights,  for  application 
to  directive  and  retrodirective  communications  was  verified.  The  transmit  antenna  patterns 
achieved  approached  the  optimal. 

o  The  directive  communication  transmit  weights  require  averaging  of  the  receive 
aperature  estimates  over  a  number  of  timeslots.  This  is  easy  to  implement  and  very  affective. 
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4.4.4  SIMULATED  PERFORMANCE  EVALUATION  (continued) 

o  For  LPI  applications,  directive  communications  dramatically  reduces  the  geographical 
area  a  reconaissance  receiver  must  be  in  to  detect  the  transmission. 

o  For  A/J  applications,  directive  communications  dramatically  increases  the  SOI  power 
available  to  the  intended  receiver  given  a  fixed  transmit  power. 

4.4.5  SUMMARY  /  CONCLUSIONS 

This  study  has  resulted  in  a  number  of  unique  developments  and  provides  supporting  proofs, 
justifications,  and  simulations. 

o  A  flexible  algorithm  was  developed  specifically  for  FHSS  and  DSSS  systems  which 
has  been  demonstrated  and  would  allow  application  to  various  Army  communication  systems. 

o  An  extremely  flexible  architecture  concept  has  been  developed  which  would  provide 
compatibility  with  Army  tactical  radios  and  allow  implementation  of  a  variety  of  adaptation 
algorithms. 

o  The  basic  computer  model  and  simulation  programs  are  established  for  future  studies, 
and  were  used  in  this  phase  to  verify  the  suggested  algorithm  and  architecture.  The  models 
allow  many  parameters  to  be  variables,  thereby  providing  flexibility  to  simulate  environmental 
conditions,  hardware  limitations,  and  various  algorithms  and  architectures. 

The  conclusion  of  this  Phase  I  effort  allows  several  options  for  further  development: 

o  A  multi-purpose  AAP  hardware  development  which  allows  testing  and  verification 
of  a  variety  of  algorithms  under  real  world  conditions  could  be  pursued.  This  hardware  could 
be  designed  to  be  compatible  with  an  existing  Army  radio.  The  Phase  2  proposal  is  based  on 
development  of  this  AAP  hardware. 

o  Additional  study  of  the  Dominant  Mode  Block  Update  algorithm  could  be  undertaken. 

o  Further  development  of  the  architecture  suggested  in  this  study,  optimized  for  a 
specific  Army  radio,  could  be  undertaken.  Commitment  to  hardware  could  occur  at  a  later  date. 

o  Additional  study  of  the  application  of  adaptive  array  processing  to  directive 
communications  could  be  pursued.  The  primary  application  is  low  probability  of  intercept 
communication,  which  shows  good  promise. 

o  Other,  less  directly  related  applications  are  also  mentioned  in  the  Phase  2  proposal  for 
future  development,  including  direction  finding. 

4.5  LIST  OF  ALL  TECHNICAL  PUBLICATIONS 

Other  than  the  documents  generated  directly  for  this  SBIR  program  (ie.  monthly  reports  and 
this  final  report),  there  have  been  no  technical  documents  generated  as  part  of  this  contract. 
A  proof  showing  some  of  the  algorithms  investigated  to  yield  maximum  likelihood  estimates 
was  started,  and  is  expected  to  be  submitted  for  publication  at  some  later  date.  A  preliminary 
abstract  for  this  paper  is  provided  below;  when  this  paper  is  submitted  for  publication,  the 
support  of  this  study  will  be  recognized  and  credited.[31] 
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4.5  LIST  OF  ALL  TECHNICAL  PUBLICATIONS  (continued) 

Preliminary  Abstract  of  Maximum  Likelihood  Estimator  Proof  - 

A  unifying  framework  for  the  design  of  blind  adaptive  signal  extraction  algorithms,  based  on 
maximum  likelihood  (ML)  estimation  of  unknown  signals  with  known  modulation  properties 
is  prsented.  The  general  ML  estimator  of  a  single  signal  of  interest  (SOI)  with  known  properties 
is  derived,  under  the  assumption  that  the  SOI  is  received  by  a  narrowband  antenna  array  in  the 
presence  of  temporally-white  complex-Gaussian  interference  with  unknown  spatial  covariance, 
it  is  shown  that  the  ML  estimator  admits  a  simple  objective  function  if  the  SOI  can  be  modelled 
as  a  member  of  a  property  set  that  encompasses  all  siganls  with  the  properties  possessed  by  the 
transmitted  SOI  waveform  or  symbol  sequence.  This  objective  function  is  shown  to  have  a 
Rayleigh  Quotient  form  with  a  closed-form  optimal  solution  if  the  property  set  is  a  linear 
subspace  that  is  not  spanned  by  the  set  of  received  data,  for  example,  if  the  SOI  is  a  time- 
limited,  frequency-limited,  FHSS  or  DSSS  waveform.  The  objective  function  is  also  shown  to 
reduce  to  a  simple  normalized  least-squares  form  with  an  iteratively  achievable  optimal 
solution  for  a  number  of  useful  nonlinear  property  sets,  for  example,  if  the  SOI  waveform  or 
symbol  sequence  has  a  constant  modulus  or  a  known  modulus  variation.  Lastly  results  are 
extended  to  more  complicated  environments  where  multiple  SOIs  are  received  in ’temporally- 
white  complex-Gaussian  interference,  and  where  a  single  SOI  is  received  in  the  presence  of 
temporally-colored  complex-Gaussian  interference.  Using  these  theoretical  results  it  is 
demonstrated  that  many  of  the  blind  and  partially  blind  array  adaptation  algorithms  developed 
to  date  can  be  interpreted  as  maximum  likelihood  estimation  techniques. 


4.6  PERSONNEL  SHOWING  ADVANCED  DEGREES  EARNED  WHILE  ON  PROJECT 
None. 
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